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Abstract: Hysteresis is a phenomenon occurring naturally in several magnetic and 
electric materials in condensed matter physics. When applied to cosmology, aka cos¬ 
mological hysteresis, has interesting and vivid implications in the scenario of a cyclic 
bouncy universe. Most importantly, this physical prescription can be treated as an 
alternative proposal to inflationary paradigm. Cosmological hysteresis is caused by 
the asymmetry in the equation of state parameter during expansion and contraction 
phase of the universe, due to the presence of a single scalar held. This process is 
purely thermodynamical in nature, results in a non-vanishing hysteresis loop inte¬ 
gral pdV) in cosmology. When applied to variants of modihed gravity models -1) 
Dvali-Gabadadze-Porrati (DGP) brane world gravity, 2) Gosmological constant dom¬ 
inated Einstein gravity, 3) Loop Quantum Gravity (LQG), 4) Einstien-Gauss-Bonnet 
brane world gravity and 5) Randall Sundrum single brane world gravity (RSII), un¬ 
der certain circumstances, this phenomenon leads to the increase in amplitude of 
the consecutive cycles and to a universe with older and larger successive cycles, pro¬ 
vided we have physical mechanisms to make the universe bounce and turnaround. 
This inculcates an arrow of time in a dissipationless cosmology. Remarkably, this 
phenomenon appears to be widespread in several cosmological potentials in variants 
of modihed gravity background, which we explicitly study for- i) Hilltop, ii) Natu¬ 
ral and hi) Golemann-Weinberg potentials, in this paper. Semi-analytical analysis 
of these models, for diherent potentials with minimum/minima, show that the con¬ 
ditions which creates a universe with an ever increasing expansion, depend on the 
signature of the hysteresis loop integral (f pdV) as well as on the variants of model 
parameters. 
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1 Introduction 

Hysteresis is a phenomenon that arises in systems with a lag between its inpnt 
and ontpnt. When this lag is dynamic i.e it changes with time, we get hysteresis 
loops. It can be evalnated by pnrely thermodynamical expressions where the ontpnt 
depends on the cnrrent and past inpnts. Snch phenomenon natnrally occnrs in 
many laboratory systems like ferromagnetic and ferroelectric materials and is often 
incorporated artihcially in several electrical systems. 

In analogy with hysteresis, in cosmology, we have the phenomenon of cyclic nni- 
verse in which the nniverse re-borns repeatedly after each cycle. Jnst as in hysteresis, 
where the material nndergoes throngh the same process over and over again, in a 
cyclic nniverse, the nniverse starts from big bang and ends in big crnnch repeatedly. 
Several models of cyclic nniverse have been proposed in literatnre [1-3]. Snch models 
also arise natnrally as exact solntions of Einstein eqnations for a closed nniverse hlled 
with perfect flnid. However in most of these models, all the cycles are identical to 
one another. Also all these models does not provide any prescription for avoiding 
singnlarity. Hence these models remain nnsnccessfnl in solving some of the major 
problems of big bang model i.e the flatness and horizon problem, and avoidance of 
singnlarity. However Tolman [4] in his paper nsed a radically different approach 
by which one conld get an oscillating nniverse with increasing expansion maximnm 
after each cycle He postnlated the presence of a viscons flnid which gave rise 
to asymmetric, irreversible eqnation of motion. This created an ineqnality between 
the pressnres at the time of expansion and contraction phases which resnlted in the 
growth of both energy and entropy. Thns he showed a novel way of linking thermo¬ 
dynamical principles to the model of cyclic nniverse. This nnnsnal approach helped 
in solving the horizon and flatness problem. In later years, theory of inflation [7-9] 
was developed which addressed both the horizon and flatness problem. Bnt none of 
these models were able to avoid big bang singnlarity. Also the model proposed by 
Tolman led to an inevitable increase in entropy with each cycle. 

However, in [10, 11], the anthors have proposed a method of avoiding both the 
presence of singnlarity and increasing entropy by nsing a cosmological analog model 
of hysteresis. The basic idea of generating the cyclic nniverse with an increasing 

^The same phenomenon occurs naturally in other cyclic universe models as proposed in refs. [5, 6]. 
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maximum remains the same and in the words of Tolman is “if the pressure tends 
to be greater during a compression than during a previous expansion, as would be 
expected with a lag behind equilibrium conditions, an element of fluid can return 
to its original volume with increased energy..". The authors created the asymmetry 
in pressure using the scalar held dynamics generated during inhationary paradigm, 
thereby maintaining the symmetric nature of the equation of motion hence avoiding 
entropy production. Originally proposed in [10] and later extended in [11], the au¬ 
thors demonstrated that “a universe hlled with a scalar held possess the intriguing 
property of ‘hysteresis " The central idea was to show that the presence of a massive 
scalar held "under certain reasonable conditions at the bounce [12-25], gives rise to 
growing expansion cycles, the increase in expansion amplitude being related to the 
work done by/on the scalar held during the expansion/contraction of the universe." 
This leads to the production of hysteresis loop dehned as f pdV, during each os¬ 
cillatory cycle. The loop area is largest in case of inhationary potentials since they 
give rise to largest asymmetry between expansion and contraction pressures. But the 
phenomenon of hysteresis is generic i.e. independent of the nature of potential. Any 
potential with a proper minima which randomizes the phase of the scalar held as it 
oscillates around the minima during expansion, thereby making all possible values of 
0 probable at turnaround, is cable of causing the phenomenon of hysteresis. However 
potentials without any proper minimum will result in a unique value of (j) which will 
make 

Pexp Pconti (^•^) 

thereby making; 

^pdV = 0. (1.2) 

Such potentials are not suitable candidates for causing hysteresis. In order to avoid 
big bang and big crunch, the authors have made used of the presence of existing 
models like brane world scenario in the early universe, and the presence of negative 
density or phantom like density in the late universe. These models replaces the big 
bang singularity by bounce and the big crunch by re-collapse or turnaround. 

Our present paper is based on the analysis done by [11]. We have further in¬ 
vestigated the phenomenon of hysteresis in different models like the variants of cos¬ 
mological constant model including ACDM, higher dimensional models like Dvali- 
Gabadadze-Porrati (DGP) brane world gravity model. Loop Quantum Gravity model 
and Einstein Gauss-Bonnet brane world gravity model in brane world, and in mod¬ 
els where the dynamics of the scalar field gets modihed which can be achieved by 
making the cosmological constant held dependent. Our aim is to study not only the 
phenomenon of hysteresis in different models but also to constrain the parameters of 
the model using hysteresis. We have mostly studied the models which can give rise to 
both the phenomena of bounce in the early universe and turnaround in the late uni¬ 
verse. We have also investigated the equivalent conditions required to achieve such 
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bouncing and re-collapsing scenarios. We have shown that our analysis holds true 
for any general form of the potential of the scalar held with a proper minimum. We 
have also shown that the phenomenon of hysteresis or the asymmetry in pressure can 
be achieved irrespective of whether the slow roll conditions of inhation are satished 
or not. A notable feature of this analysis is that an increase in expansion maximum 
after each cycle now depends not only on the sign of f pdV but also on the parame¬ 
ters of the models that we have considered. Thus we see that using the remarkable 
cosmological effect of hysteresis as proposed by [10, 11], there are numerous methods 
and models in which a cyclic universe with an ever increasing amplitude maximum 
can be achieved. 

The plan of the paper is as follows: 

• In section 2, we have discussed the mathematical formulation leading to the 
phenomenon of hysteresis in cosmological scenario. 

• In section 3 we have tried to draw an analogy between the hysteresis in ferro¬ 
magnetic materials and in a cyclic universe. This simple analysis leads to the 
conclusion that the equation of state parameter w plays the role of magnetic 
held H in cosmology and the scale factor mimics the behavior of magnetization 

M. 

• In section 4, 5, 6, 7 ,9.1, we have studied the behavior of each of the above 
mentioned models for hat, closed and open universe and have discussed various 
unexplored important cosmological features. 

• In section 4, 6, 7 ,9.1, we have explicitly studied cosmological hysteresis in 
the context of higher dimensional theories, where we have shown the results 
for both space-like and time-like extra dimensions. For the sake of simplicity 
we restrict ourselves up to hve dimensions {D = 5). But one can extend the 
computation for dimensions, D > 5. 

In this paper, we have drawn various physical conclusions by explicitly solving the 
equations governing the dynamics of the system using the semi-analytical techniques. 
Though the analysis is perfectly true for any kind of cosmological potential with a 
proper minimum/minima, we have the studied the detailed features for three different 
potentials - hilltop potential, natural potential and Colemann-Weinberg potential. 
All these potentials have well dehned minimum/minima and have free parameters 
which can be adjusted to get the required results. Thus this analysis helps us to 
put stringent constraints on the characteristic parameters of these models in the 
bouncing scenario along with cosmological hysteresis. Though the analysis that we 
have performed holds good under certain physically acceptable approximations and 
limiting cases, but we can at least show mathematically if there are any limiting cases 
in which these potentials combined with the models can give rise to the phenomenon 
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of cosmological hysteresis i.e. make j> pdV non-zero. In this paper we have also 
explicitly derived the expression for work done in one complete cycle of expansion 
and contraction, and have shown it to be non zero. But the sign of the integral 
depends on how we have chosen the sign and magnitudes of the parameters of our 
models. The most interesting result of our analysis is that there are several models 
which can give rise to a cyclic universe with an increasing amplitude of expansion. 

2 Basics of cosmological hysteresis 

Before going to the technical details of the “Cosmological hysteresis” let us mention 
clearly the underlying assumptions. It is important to note that throughout the 
analysis in this paper we assume that during a specihed period in the time line of the 
universe, it is described by a single massive scalar held which is minimally interacting 
with the gravity sector. The presence of this scalar held is responsible to generate 
the required asymmetry in the pressure of the universe and this leads to an overall 
increase in the energy of the universe and hence an increase in its amplitude of the 
expansion rate. 

The action governing the dynamics of this scalar held 0 with potential V (0) 
within ehective held theory description, is given by 

^ = j =Seh + S^ (2.1) 

where the signature of the metric throughout our analysis is (-, +, +, +). The total 
action S can be written as the sum of the standard Einstein Hilbert action {Seh) 
and the action for the scalar held (5'^). The energy-momentum tensor for the scalar 
held can be computed from the matter part of the action 

S^ = J d^x^ ( 2 . 2 ) 

as: 

r* = j (2.3) 

and for a homogeneous and isotropic spatially hat {k = 0) FLRW cosmological 
background, the energy density and pressure for a scalar held can be computed from 
the energy momentum tensor as: 

p = ly + V {< P ), (2.4) 

. (2.5) 

-4 - 



and the resulting equation of state parameter w can be written as: 


p _ - 21 /((/>) 

P ~ (^2 + 21 /( 0 ) ’ 


(2.6) 


where we assume that the energy momentum tensor for scalar held can be approx¬ 
imated via perfect huid. Also in the spatially hat FLRW cosmological background 
the scalar held equation of motion is given by: 


□0 + 


dV 

dcj) 


4 > + 3i/0 -|- 


dV 

d(p 


0 , 


where □ is the d’Alembertian operator four dimension dehned as: 


(2.7) 




A 

df^ 


d 


(2.8) 


Here H is the Hubble parameter dehned as: 


1 da(t) 
a{t) dt 


(2.9) 


where a{t) is the scale factor. 

Origin of cosmic hysteresis: If we closely look at Eq (2.7), then we can conclude 
that when the universe expands i.e. H > 0 the second term 3H(j) mimics the role of 
friction and opposes the motion of the scalar held, thus serves the purpose of damping 
during its motion. This lowers the kinetic energy of the scalar held compared to its 
potential energy, giving rise to a soft equation of state (P = —p in case 0^/2 << H(0) 
i.e. slow roll regime). By contrast, in a contracting [H < 0) phase of the universe, 
the term 3P0 behaves like anti-friction and favors the motion of the scalar held 
and hence accelerates it. This makes the kinetic energy of the scalar held much 
larger than the potential energy, giving rise to stih equation of state (P = p in case 
0^/2 >> H(0)). As a result, from second law of thermodynamics, we can convey that 
a net asymmetry in the pressure (during expansion and contraction cycle) leads to a 
net non-zero work done by/on the scalar held. In addition, if we now postulate the 
presence of bouncing and recollapsing mechanisms during contraction and expansion 
respectively, one can expect that a non-zero work done during a given oscillatory 
cycle to be converted into expansion energy, resulting in the growth in the maximum 
amplitude and hence maximum volume of the universe of each successive expansion 
cycle. Thus producing older and larger cycles. In [11], using simple thermodynamic 
arguments, the authors have developed the equations which relate the change in 
maximum amplitude of the scale factor after successive cycles to the work done. The 
authors have shown that these equations have a universal form which is independent 
of the scalar held potential responsible for hysteresis. As has been discussed in 
[11], though the process is independent of the potential, “the presence of hysteresis 
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is closely linked to the ability of the held 0 to oscillate”. This urges the presence 
of potential minimum/minima. This is because only potentials with well dehned 
minimum can make the held oscillate. As a result Oscillations of the scalar during 
expansion makes its phase arbitrary, thereby making all values of 0 equally likely 
at turnaround. Thus assuring that the values of 0 and 0, when the universe turns 
around and contracts, are nearly uncorrelated with its phase space value when the 
held 0 began oscillating. As a result, the held almost always rolls up the potential 
along the diherent phase space trajectory compared to the one along which it had 
descended during expansion. This gives rise to unequal pressure during expansion 
and contraction hence to non-zero work done. 

In the present work, we have reconsidered the above phenomenon of hysteresis 
and applied it to models or physical mechanisms which can make the universe bounce 
and turnaround in the presence of potentials having well dehned minimum/minima. 
Further solving Eq. (2.5) and Eq. (2.7) simultaneously, along with the Friedmann 
equations derived from various cosmological background model, and applying the 
proper bounce and turnaround conditions, which will be discussed in the following 
sections, we get the explicit expressions for the scalar held and scale factor as a 
function of time. 

It has been hrst pointed out in ref. [10] that when we plot the equation of state 
given by tc = p/p vs the scale factor from a specihed cosmological model, we get a 
hysteresis loop whose area contributes to the work done by/on the scalar held during 
expansion and contraction of the Universe. The general expression for the work done 
by/on the scalar held during one cycle is given by 


pdV = / pdV + / pdV 


' cont 


’ exp 


( 2 . 10 ) 


where the contributions f ,pdV and f pdV represent the work done by/on the 
scalar held during the phase of contraction and expansion of the Universe respectively. 
The signature of the integral depends on the pressure during the phase of contraction 
and expansion i.e if the pressure corresponding to the contraction phase is greater 
than the pressure due to expansion phase i.e. 


Pcont ^ Pexp 

then the overall signature of the p — dV work done is negative i.e. 

pdV < 0. 


( 2 , 11 ) 


( 2 . 12 ) 


On the other hand, if the pressure corresponding to the contraction phase is smaller 
than the pressure due to expansion phase i.e. 


Pcont ^ Pexp 


(2.13) 
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then the overall signature of the p — dV work done is positive i.e. 


pdV > 0. 


(2.14) 


For a universe characterized by a scale factor a(t), the total volume at any given time 
is given by a^{t) (neglecting the overall constant factor). Hence using this input, the 
area of the cosmological hysteresis loop or equivalently the p — dV work done is given 
by; 
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/ pa^da + 

/ pa^da 

_J cont 

J exp 


Hence we follow the following algorithm; 


(2.15) 


• In the context of various cosmological frameworks i.e. DGP branewold gravity, 
Loop Quantum gravity (LQG), Einsitein-Gauss-Bonnet (EGB) gravity and in 
presence of pure and held dependent cosmological constant within Einstein 
gravity we explicitly derive an expression for this characteristic integral in terms 
of the scale factor and the scalar held degrees of freedom. 

• Then we solve the equations of motions for the scalar held under various limit¬ 
ing approximations and hence elaborately study the physical conditions under 
which the above integral gives non zero value. 

• Further we repeat this mentioned two step process for three diherent cosmo¬ 
logical potentials i.e. for Hilltop potentials. Natural potential and Goleman- 
Weinberg potential within the framework of ehective held theory prescription, 
which will be discussed in the next section in detail. 

• We also analyze the whole process graphically to see whether we get a net 
increase in the amplitude of the scale factor after one cycle of expansion and 
contraction. 


Though the phenomenon of hysteresis is commonly attached to magnetic and electric 
systems, its appearance for diherent cosmological models, almost naturally, makes 
us appreciate and acknowledge its importance in the held of cosmology. The various 
advantages which a cyclic universe along with the phenomenon of hysteresis, which 
is our main focus in this paper, have are; 

• The phenomenon of hysteresis is important due to the simplicity with which it 
can be generated. Only a thermodynamic interplay between the pressure and 
density, creating an asymmetry during expansion and contraction phase of the 
universe, succeeds in causing cosmological hysteresis. 
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• As the phenomena of cosmological hysteresis deals with the bonncing as well 
as the re-collapsing phase of the universe, one can avoid the appearance of Big 
Bang Singularity as well as the Big Crunch at early and late times. 

• Hysteresis can be generated by the presence of a single massive scalar field, 
which has already been studied for a wide variety of physical situations. The 
most exciting issue is, it do not require any other fields for its occurrence. 
Hence it is very easy to handle and its properties can be studied extensively in 
cosmological literature. 

• A cyclic universe having the conditions to cause hysteresis, can solve all the 
Big Bang puzzles, hence acts as an alternative proposal to inflation [27-31]. 

• In this scenario, we can always start with a closed or open universe and after 
allowing the universe go through a number of cycles, we get the present observ¬ 
able flat universe. In this paper, we will extensively deal with several models 
that can give rise to such cyclic universe. Through analytical calculations, we 
will show that irrespective of the nature of the universe, we can get a cyclic 
model with increasing amplitude of the scale factor for a wide variety of models. 

• It results in dissipative cosmology [10, 11, 26], which makes the whole process 
irreversible, which finally causes an arrow of time according to widespread 
belief. But recently in the ref. [35] the authors have explicitly shown that 
for cosmological hysteresis phenomena such an arrow of time can appear even 
if equations describing cosmological evolution are dissipationless, which makes 
the cumulative process reversible provided they possess cosmological attractors 
in the expanding and contracting phase of the universe. Additionaly, it is 
important to note that for flat cosmological potential the increament in the 
expansion cycles are routinely observed. But for the steep potentials two fold 
cyclic pattern with lesser expansion cycles is observed, which being nested in 
the larger expansion cylces. This phenomena is exactly analogous to the ‘beat’ 
formation in acoustics systems. In ref. [10, 11] the authors have studied the 
cosmological beat formation in the context of chaotic potential H(0) = 

• Cosmology in this scenario has not been explored in a very wide sense. Earlier 
it has been studied by the authors of refs. [10, 11, 35]. In this paper, we have 
further explored this phenomenon for several other models which have not been 
discussed earlier in refs. [10, 11, 35]. 

• In future we plan to connect these analyses with CMB observations, by rigorous 
study of the cosmological perturbation theory [32-34] in various orders of metric 
fluctuations and computation of two point correlations to get the expressions 
for scalar and tensor power spectrum in this context. Hence we extend the 
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study of this paper to compute the primordial non-Gaussianity in CMB from 
three and four point correlations [36-45]. We also plan to derive the explicit 
expression for various modihed consistency relations between the non-Gaussian 
as well as other cosmological parameters in the present context. 

• Also we plan to setup the formalism of magnetogenesis from the present setup 
which can be treated as the alternative version of the inflationary magnetoge¬ 
nesis [46-48]. 

• We also carry forward our analysis in the development of density inhomo¬ 
geneities, which is the prime component to form large scale structures at late 
times. Also the specihc role of cosmological hysteresis in the study of cosmolog¬ 
ical perturbations i.e. for interacting/decoupled dark matter and dark energy 
have not been explored at all earlier. We have some future plan to do some 
computations from this setup. 

• Further using the reconstruction techniques [49, 50, 52, 53] we want to study 
the generic features of scalar held potentials in the framework of cosmological 
hysteresis. 

• Last but not the least, we also plan to further study this outstanding cosmologi¬ 
cal phenomenon for different modihed gravity pictures i.e. for variants of f{R) 
gravity [54-57], two brane-world model in presence of the Einsitein-Hilbert 
term and the Einstein-Hilbert-Gauss-Bonnet gravity setup [58-64]. 

3 Analogy with magnetic hysteresis 

The phenomenon of lagging of magnetic induction B or magnetization M behind the 
magnetic held H when a specimen of a magnetic material (such as iron) is subjected 
to a cycle of magnetization is called “magnetic hysteresis”. The closed loop that is 
traced by the material in the B — H or M — if plane is known as the hysteresis 
loop. It is related to the change in the alignment of the magnetic dipoles, as one 
varies the magnetic held which leads to magnetization and demagnetization of the 
material. It is a beautiful way of depicting the ehect of varying H on the system. 
Fig. 1(a) is a representative schematic diagram in which we have explicitly shown 
one such hysteresis loop in ferromagnetic material. Since we want to draw an analogy 
between the phenomenon of hysteresis in magnetism and cosmology, we have also 
shown Fig. 1(b) to illustrate how an ideal hysteresis loop looks in the case of a cyclic 
universe. In a cyclic universe, the hysteresis loop is traced by the universe in the 
a — w plane where a is the scale factor of the universe and w = p/p is the equation 
of state parameter. Glosely following the two hgures we can say that just as in a 
magnetic material, we also vary the magnetic held and try to hnd the behavior of 
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the magnetization of the material, in this cyclic model, we vary the pressure or the 
equation of state in order to find the variation of the scale factor or expansion of 
the universe. This motivates us to consider that the parameter playing the role of H 
in cosmological hysteresis is w and that oi M/B is a. Just as in magnetic material, 
where the magnetization oscillates within a specihed maximum and minimum values, 
analogously in cyclic model, the scale factor goes through maximum and minimum 
values which has been clearly shown in [ 10 , 11 ] for few cosmological models and will 
again be shown in this paper for some other cosmological models. While the primary 
cause of magnetic hysteresis is the asymmetry in the behavior of the alignment of 
the magnetic dipoles with increase and decrease of iJ, the cause for hysteresis in 
cosmology is the asymmetry in pressure during expansion and contraction phases of 
the universe. In magnetic hysteresis loop, the minimum and maximum corresponds 
to the state when all the magnetic dipoles are aligned in reverse and along the 
direction of H respectively. In cosmological hysteresis, the maximum is reached 
when the scale factor reaches its maximum value and the density of the scalar held 
reached at its minimum i.e when the condition for re-collapse is generated, and the 
corresponding minimum is reached at bounce when scale factor becomes minimum 
and density of the scalar held reached at its maximum value. But unlike in magnetic 
hysteresis where the parameters can take all kind of values i.e positive, negative and 
zero, in cyclic universe, one of our primary goal is to avoid singularity, i.e. a 7 ^ 0 , so 
that singularity is replaced by bounce. While repeated cycles of magnetic hysteresis 
loop results in loss of hnite amount of energy, hence decrease in the area of the 
loop, in cosmological hysteresis, repeated cycles may be either larger or smaller in 
area compared to the previous one. Thus, we hnd that a closer look at these two 
phenomena draws lots of similarities in their behavior. In table. 1 we have depicted 
the analogy and the parallelism between the diherent features of these two kinds of 
hysteresis. 

In this section we have shown the complete analysis of different models and 
studied the conditions under which we get an increase in the amplitude of the scale 
factor after each cycle. We have studied open closed and flat universe. We have also 
explicitly calculated the work done in one cycle for three different potentials. 

4 Hysteresis from Dvali-Gabadadze-Porrati (DGP) brane world 
gravity model 

The DGP model is a model of modified gravity theory proposed by Gia Dvali, Gregory 
Gabadadze, and Massimo Porrati. The model consists of a 4D Minkowski brane 
embedded in a 5D Minkowski bulk like Randall Sundrum (RS) II model. But unlike 
in RSII model [65, 66 ] here the inhnitely large 5th extra dimension is flat. The 

^See also the details of Randall Sundrum (RS) I model studied in ref. [67] for completeness. 
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(a) An illustration of hysteresis in ferromagnetic materials. The loop has been plotted in the 
B — H plane. 



(b) An idealised illustration of cosmological hysteresis. The loop has been plotted in the w — a 
plane. 

Figure 1. Analogy between magnetic hysteresis in ferromagnetic materials and cosmo¬ 
logical hysteresis. 
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Cosmological Hysteresis 

Loop parameters 
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H 

w = p/p 
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§ w da {^ § p dV) 
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Cyclic Universe 

hysteresis loop 

ferromagnetic materials 

with inflationary conditions 

Upper and lower 

No limit on maximum 

'^max = H“ 1 

limits of the loop 

and minimum values of H 

'^min — 1 

Characteristic equation 

B = g{H + M) 

w = p/p 

for the loop 



Nature of 

Soft ferromagnetic materials 

Universe with softer equation 

the loop 

have smaller loop area 

of state have smaller loop area 


Table 1. Table showing the analogy between magnetic hysteresis and cosmological hys¬ 
teresis. 


Newton’s law can be recovered by adding a 4D Einstein-Hilbert action sourced by 
the brane curvature to the 5D action. While the DGP model recovers the standard 
4D gravity for small distances, the effect from the 5D gravity manifests itself for large 
distances. 

The DGP model is described by the following action [68]; 

/ d^x^R + I (4.1) 

where qab is the metric in the 5D bulk and 

= d^X^d,X^~gAB (4.2) 

is the induced metric on the brane with X^[x^) being the coordinates of an event 
on the brane labeled by x^. The hrst and second terms in Eq. (7.2) correspond to 
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Einstein-Hilbert actions in the 5D bulk and on the brane, respectively. Note that 
and are 5D and 4D gravitational constants, respectively, which are related 
with 5D and 4D Planck masses, M 5 and M 4 , via 



K 


2 


1 

wr 

1 _ 1 


SttG. 


(4.3) 

(4.4) 


The Lagrangian describes matter localized on the 3-brane. 

The equations of motion read 


_ n 

^AB — U 5 


(4.5) 


( 5 ) 

where is the 5D Einstein tensor. The Israel junction conditions on the brane, 
under which a Z 2 symmetry is imposed: 


G(^K ^ 1 / ) 

T r 


(4.6) 


where is the extrinsic curvature calculated on the brane, is the energy- 
momentum tensor of localized matter and Vc is the crossover length scale. 


Tr. = 


2 M|’ 


(4.7) 


because it sets the scale above which the effect of extra dimension becomes important. 
The modihed Friedmann equations in this model are given by [69]: 

2 





H + H^ = --{p + p) 
6 


iP 




4^2 


-1/2 


+ 


+ "TW + X- 

3 dr? 2r. 


^2p 


(4.8) 


(49) 


In the present analysis, we have set the extra dimension as time-like which is neces¬ 
sary for getting late time acceleration in DGP model, but in generalized prescription 
one can comsider both space and time like extra dimensions. Also it is important to 
note that, in Eq. (4.9) the spatial curvature k can take values 0 or ±1. Now following 
the proposal of [ 10 , 11 ], we know that in order to get a cyclic universe, the condition 
for bounce and turn around are given by: 


• Bounce: 


and 


H 


a 

a 


a = 


1 da{t) 
a(t) dt 

dt^ 


(4.10) 

(4.11) 
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• Turn around: 


(4.12) 


d 1 da{t) 
a a(t) dt 


and 



(4.13) 


In the following subsections we will discuss all of these possibilities in detail for DGP 
brane world gravity framework. 

4.1 Condition for bounce 

In the early universe, at high energy or equivalently in the high density regime of 
the brane world, pr^/M| 3> 1, one can expand Eq. (4.9) as [74] 




At first order approximation, Eq. 4.14 reduces to the following expression: 



(4.15) 


This is obviously a valid assumption because at the time of bounce, the density of 
the matter content of the Universe (in our case scalar held) is at its maximum, hence 
we can neglect the contribution from the other higher order terms in Eq. (4.14). 
Now applying the pr^/M| 3> 1 at bounce, and setting 77 = 0 in Eq. (2(c)), we get: 



(4.16) 


where pb and ab are the density and scale factor at bounce, respectively. Similarly 
the mass content at bounce (neglecting the constant factor) is given by. 




Hence applying the energy conservation in the present context, we get: 


(4.18) 


SM + SW = Q, 


where 5W is the work done during each expansion-contraction cycle which is given by 
f pdV which includes contribution from the area of the hysteresis loop. By setting. 



(4.19) 
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we get the expression for change in amplitude of the scale factor at each successive 
cycle as 


6a 


min 




f pdV 


CLb 


3 ( ^ 

at 



Qb — 2 


(4.20) 


Therefore we clearly observe that the change in amplitude of the scale factor after 
each cycle depends on the sign of the integral, the curvature parameter and the 
cross over length scale for DGP brane world. It is independent of the density of the 
matter content of the universe. During our computation we also observe that the 
if we fix the spatial curvature parameter k = —1, this makes the energy density of 
the scalar held imaginary, which is not at all physically possible. Hence, for k = —1 
the cosmological bounce from DGP brane world model is not at all possible. But if 
we neglect all the higher order terms, Eq. (4.14) reduces to the standard Friedmann 
equation for which the condition for bounce becomes 


Pb 


3kMi 

Ub ’ 


(4.21) 


which is possible for both k = ±1 but not for A; = 0. Hence we see that the bouncing 
condition depends largely on the order of terms which we are including because 
our analysis is possible only under approximations. But since our present observed 
universe is hat (/c = 0), we will be interested to study the hrst case in more detail. 
Hence in our further analysis for DGP model, we will show the results for Eq. (2(c)) 
only i.e exclude the k = —1 possibility. Here we have the following expression for 
6amin- 


f pdV 


^ \a,b 


J_ 

2 rr 


ttb 




for k = +1 


6am,in ^ 


(4.22) 


f pdV 

--- 7 ^ for A: = 0. 

i ^(* + 2) 

Let us now briehy mention the characteristic feature of the results for cosmological 
bounce for DGP brane world model in the following: 


1. For a closed universe, depending on whether the quantity appearing in the 
denominator for A: = +1 is positive or negative, an increase in the scale factor 
after each cycle is possible if: 


The hysteresis loop integral 


pdV < 0 


(4.23) 
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or equivalently 


Pexp ^ Pconti 


(4.24) 


• The hysteresis loop integral 


pdV > 0 


(4.25) 


or equivalently 


Pexp ^ Pcont- 


(4.26) 


2. On the other hand, for the case fc = 0, the quantity in the denominator is 
always positive and finally we get an increase in the scale factor only if the 
hysteresis loop integral 


pdV < 0. 


(4.27) 


4.2 Condition for acceleration 


From Eq. (4.9), at bounce or at high energy, contribution from l/vc term is small 
compared to the density of the scalar field and consequently one finally gets the 
standard cosmology from the condition for acceleration as: 

P<-^ (4.28) 

which clearly implies that the cosmological bounce can be obtained by violating the 
strong energy condition. 

Substituting the expression for ph for k = +1, 0, we get the conditions for accel¬ 
eration at bounce as: 




-Ml 




+1 


5pb < < 


(4.29) 


V 


Ml 

4^2 


for A: = 0. 


This implies that the pressure of the matter content of the Universe at the time 
of cosmological bounce is related to the scale factor and the cross over length scale 
Tc. Also we observe that the condition for acceleration depends on whether we have 
considered a possibility of a closed or flat universe. For flat universe, we see that the 
condition for acceleration is independent of the scale factor at cosmological bounce. 

Finally substituting the expressions for the energy density p and pressure p from 
Eq. (2.5) into Eq. (9.9), we get the condition for the expansion of the Universe in 
terms of the scalar field degrees of freedom as: 


0 " < V{ct>) 


(4.30) 
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which is same as that we get for the standard cosmological inflationary scenario. 
Therefore, we can conclnde that the conditions on the dynamics of scalar which is 
needed for cansing acceleration in standard case remains nnchanged for the DGP 
model at the time of cosmological bounce. This implies that the cosmological poten¬ 
tials in the standard cosmological scenario, capable of satisfying the above condition, 
will also be able to cause acceleration in a universe described by the DGP brane world 
model. 

In Fig. 2, we have shown the phenomena of bounce and acceleration in the DGP 
model. We can draw the following conclusions from the above figures: 


• In Fig. 2(a) and Fig. 2(b), we have plotted the r.h.s of Eq. (4.9) using the 
relation p = The graphs have been plotted only for the cases when 

k = 0, 1. This is because, Eq. (4.16) shows that density becomes imaginary 
for the case when k = —1, which is unphysical. 


• We have used tc = 0, since we require a soft equation of state for causing the 
acceleration and expansion. The case w = 1/3 made H only approximately 
zero. The value of has been chosen such that we get proper bounce and 
acceleration. Thus we find that we require a larger value of Tc for causing 
bounce in a closed universe as compared to a flat universe. 


• From Fig. 2(a), we get the bounce at p = = 0. Similarly from Fig. 2(b), we 

get p = pb = 7.0M^. 


• In Fig. 2(a) and Fig. 2(b), H going to negative values may be interpreted as 
the universe changing its direction of motion at bounce. As had been discussed 
in [10], the condition of bounce/turnaround can be imposed by either the con¬ 
dition of making the scale factor changing sign with other quantities remaining 
same, or, a going to negative values with other quantities remaining same. 


• Figs. 2(c) and 4.31 show the necessary condition of acceleration (a > 0) at 
the time of bounce. Here we have plotted the r.h.s of Eq. (4.9) and Eq. (4.9). 
These plots have also been obtained for the same parameter values as the earlier 
graphs. 


• Thus Fig. 2 shows graphically that the phenomenon of bounce is possible for 
dgp model having tc = 0. 
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H vs p for DGP model (k=0) H vs p for DGP model (k=1) 




(a) An illustration of the bouncing condition for (b) An illustration of the bouncing condition for 
a universe with an equation of state u; = 0, k = a universe with an equation of state w = 0, k = 
0, I r, |= 1. 1, I r, 1= 1.44. 



(c) An illustration of the acceleration condition (d) An illustration of the acceleration condition 
at the time of bounce for a universe with an equa- at the time of bounce for a universe with an equa¬ 
tion of state w = 0, k = 0, | rc |= 1. tion of state w = 0, k = 1, | Tc |= 1.44. 

Figure 2. Graphical representation of the phenomena of bounce and acceleration for 
DGP model. 
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4.3 Condition for turnaround 


To establish the condition for turnaround we can again rewrite Eq. (4.9) as [74]; 

k 




1 

Aip2 

c 

[V SMi 

1 

4^2 

V 3M|y 


+ 1 


+ 24 


'i + i^d + i 
3Mi ^ 


(4.31) 


In late time universe, pr^/M| << 1, hence keeping terms upto the hrst order we get: 




2p 

3M2 


(4.32) 


This is also a valid assumption because at the time of re-collapse, the density of the 
Universe is low or reaches at its minimum but not negligible. Hence we need to keep 
terms upto the hrst order in the binomial series expansion. Now by setting 77 = 0, 
we get: 


Pt = 


3Mi 

2 



(4.33) 


where pt and at are the density and scale factor at turnaround respectively. 

Similar to the Cosmological bounce case, here we also can equate the change 
in energy or mass content of the Universe to the work done after each expansion- 
contraction cycle and hnahy we get: 


6Mt = 


3M| 

2 



(4.34) 


Therefore the change in the scale factor after each successive cycle is given by the 
following expression: 


6a 


max 


2 

3M| 



(4.35) 


Therefore, just like in the case of bounce, where the change in the amplitude is de¬ 
pendent on the parameters of the cosmological model, in such a physical prescrption 
the change in the amplitude of the scale factor at turnaround also depends not only 
on the work done, but also on the cross over length scale Vc of the DGP brane world 
model. 
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f 


2 


^dV for k = +\ 


SMi 



SQ/m,in \ 


2k 




pdV 


for A: = 0 (4.36) 


3M2 1 + 


3al, 


pdV 


for k = —1. 


Let us now briefly mention the characteristic feature of the results for turnaround 
for DGP brane world model in the following: 


1. For A; = — 1 we clearly observe that for an open universe, an increase in ampli¬ 
tude of the scale factor after each successive cycle is possible if the hysteresis 
loop integral 


pdV > 0. 


(4.37) 


2. For A; = 0 we observe that the increase in magnitude of the scale factor depends 
not only on work done, but also on the cross over length scale. We get a 
positive change in the scale factor with each successive cycle if the hysteresis 
loop integral 

^ pdV > 0. (4.38) 

Hence, this is a case where increase in expansion maximum is possible only if 
the work done is positive. 


3. For A; = -|-1 if the denominator 



(^1 ^“7"^ > 0 

(4.39) 

or equivalently 

3r;2 

^^max ^ -| 

c 

(4.40) 

we get positive dttmax 

if the hysteresis loop integral 



A 

o 

(4.41) 


and vice versa. Thus depending on the relative magnitude of the scale factor 
and cross over length scale at turnaround, we can get an increase in expansion 
for both positive and negative signature of the work done. 
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4.4 Condition for deceleration 

To establish the condition for deceleration we hrst take the time derivative of Eq. (3(a)), 
and using the energy conservation or equivalently the continuity equation, 

p + 3i7(p + p) = 0, (4.42) 


we can write: 


{P + P) 

Ml 


1 

H—^ + 


2p 

3Mf 


From the above equation we get the condition for deceleration as: 


{pt + 3pi) > 


3M| 


(4.43) 


(4.44) 


Therefore, turnaround can be obtained without violating the energy condition. And 
just like for acceleration, here we see that the condition for deceleration depends on 
Tc, which was expected because in the late universe, the effect of will become more 
important. 

In place of Eq. (4.29), we get the conditions for deceleration at turnaround as: 


Ml 




for k = +1 


Pt= < 



for A: = 0 


(4.45) 


, (r? + 2o?) '■ 

Finally substituting the expressions for the density p and pressure p from Eq. (2.5) 
into Eq. (4.44), we get the following condition for expansion of the Universe in terms 
of the scalar field as: 

^ ( 4 - 46 ) 

This clearly implies that the lesser the contribution from the potential energy as 
compared to the standard cosmological infiationary case or equivalently for the ac¬ 
celeration case as mentioned earlier, easier will be to achieve contraction phase of the 
Universe. We also clearly observe that, even if the potential satisfies the deceleration 
condition, just like as for standard cosmological scenario, it is not necessary that it 
will also cause deceleration in DGP brane world model at turnaround as of now the 
deceleration is also the cross over length scale Tc dependent in the present context. 

In Fig. 3, we have shown the phenomena of turnaround and deceleration in dgp 
model. We can draw the following conclusions from the above figures: 
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H vs p for DGP model (k=0) H vs p for DGP model (k=1) 




(a) An illustration of the turnaround condition (b) An illustration of the turnaround condition 
for a universe with an equation of state w = for a universe with an equation of state w = 
1, k = 0, I r 1= 10. 1, fc = 1, I 1= 1.44. 



(c) An illustration of the deceleration condition (d) An illustration of the deceleration condition 
at turnaround for a universe with an equation of at turnaround for a universe with an equation of 
state u; = 1, A: = 1, | Tc |= 10. state w = 1, k = 1,\ rc \= 1.44. 

Figure 3. Graphical representation of the phenomena of turnaround and deceleration for 
DGP model. 
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• In Figs. 3(a) and 3(b), we have plotted the r.h.s of Eqns. (4.31) and (4.9) 

respectively, using the relation p = But for this case we have used 

w = 1, since we require a stiff equation of state for causing contraction and 
deceleration. 

• From Figs. 3(a) and 3(b), we get the turnaround at p = pt = 0 and at 
p = Pt = O.OSMp respectively. 

• Figs. 3(c) and 3(d) show the necessary condition of deceleration (a < 0) at the 
time of turnaround. Here we have plotted the r.h.s of Eqns. (4.9)and (4.31), 
(4.9) respectively. This plot has also been obtained for the same parameter 
values as the earlier graph. 

• Thus Fig. 3 shows graphically that the phenomenon of turnaround is possible 
for the DGP model having w = 1. 


4.5 Evaluation of the work done for one cycle 

Let us now compute the explicit contribution from the work done for aone complete 
cycle for DGP brane world cosmological setup. To serve this purpose we start with 
the following closed loop integral ; 


pdV = / pdV + / pdV 

J cont J exp 


mm 


-^max 


pdV 


min 


pdV 


(4.47) 


where i and i — 1 refer to the two successive cycle i.e. ith. and i — 1th cycle of expan¬ 
sion and contraction phase of the Universe. In the present context the work done 
corresponding to the expanding phase and the contracting phase can be expressed 
in terms of the work done between two successive cycle as: 



i-1 

min 


i-1 

max 


pdV, 


i-1 

min 


pdV 


(4.48) 

(4.49) 


where (aj^^c, are the maximum magnitude of the scale factor for i — 1th and 

ith cycle of the Universe. Similarly, (a^in) represents the the minimum magnitude 
of the scale factor for i — 1th cycle of the Universe for DGP brane world model. 
Now the volume of the Universe can be written as. 


V = a 


3 


(4.50) 
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and an infinitesimal change in the volume can be written as, 


dV = 3a^da = 3a^(t) dt = Sa^ddt, (4-51) 

which is frequently used for further computation of the work done. And we also 
know from Eq. (2.5), in the presence of scalar held, the pressure for the scalar held 
can be written as: p = 0^/2 — V{(j)). Substituting these in Eq. (4.47), we get 





a^ddt + 





a^ddt 





'min 




a^ddt 


(4.52) 


Further using the solution of scalar held 0 and the scale factor a, one can solve the 
above equation to get the estimation of this integral. Also it is important to mention 
here that the above integral can also be expressed in terms of scale factor only using 
Eq (4.9), Eq (4.9) and Eq (2.5). But since the Friedmann equations in case of DGP 
brane world model are highly complicated, we can use the late time and early time 
approximations of the Friedmann equations in order to get an physically relevant 
approximate analytical expression for the integral for the work done. Therefore, the 
work done can be decomposed into four parts as follows: 


pdV = 


pdV 


mm 


pdV + 


pdV 


pdV 


Contraction 


Expansion 


(4.53) 


where the hrst two terms corresponds to late and early times during the period of 
contraction respectively and the last two terms corresponds to early and late times 
during the period of expansion respectively. Here a' corresponds to the scale factor 
at the time of transition t' from early to late time or vice-versa, Qmax and Omm 
corresponds to the values of the scale factor at the time of turnaround tmax and 
bounce t„iin respectively. 

At early time, from Eq. (4.14), keeping terms upto hrst order and neglecting the 
contribution from 1/r^ terms in the acceleration equation [74], we get: 



d _ {p + 3p) 
a QMl 


(4.54) 

(4.55) 
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Hence using Eq. (4.54), Eq. (4.55) and Eq. (2.5) for the pressure content of the 
Universe p, we get: 


2 


-V(4,) 




(4.56) 


Let this equation be valid upto to a cut-off time scale t' where the value of the scale 
factor is a'. Therefore using Eq. (4.56) into Eq. (4.52), the second and third integral 
as appearing in Eq. (4.53) can be expressed as: 


pdV 



At late time, using Eq. (3(a)) and Eq. (4.44), we get: 


_ 3M| /ay k _ 1 

^ 2 \a/ rl 

a ^ jp + p) , ^ 
a Ml rl 3M| 


(4.59) 

(4.60) 


Further using Eq. (4.59), Eq. (4.60) and Eq. (2.5) for the pressure content of the 
Universe p, we get: 


2 


U(0) = Ml 



d 

a 


1 / a\^ k 

2 Vay ~ ^ 


(4.61) 


This equation is also valid upto to a cut-off time scale t' where the value of the scale 
factor is a'. Therefore using Eq. (4.61) into Eq. (4.52), the hrst and last integral as 
appearing in Eq. (4.53) can be expressed as: 




dt 


(4.62) 

(4.63) 


Therefore by substituting Eq. (4.58) and Eq. (4.63) into Eq. (4.53) the complete 
expression for the work done in a single expansion-contraction cycle is governed by 
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the following expression: 





(4.64) 


In the present context, we clearly visualize from our analysis that, by knowing the 
solution of the scale factor in the early Universe and as well as in the late Universe, 
we can get an idea of the nature of the hysteresis loop. Additionally it is important 
to note that while the evaluation of work done is independent of the parameters of 
this model at early times, parameter dependence enters through late time evaluation 
of the integral for work done. 


4.6 Semi-analytical analysis for cosmological potentials 

In this section, we try to hnd simple analytical expressions for the work done during 
one cycle of expansion and contraction from various cosmological models in order 
to get some idea of the behavior of the cosmological hysteresis loop. Though the 
analysis is independent of any particular functional form of potential, but in order 
to study the physical signihcance as well as the nature of the derived results in the 
previous section, we need to specify the functional form of the cosmological potential. 
While doing the analysis for DGP brane world model and for all the other subsequent 
cosmological models, we will consider three different potentials i.e. Hilltop, Natural 
and Coleman-Weinberg potential, which will be discussed in the next sections in 
detail. 

Since from Eq. (4.64), we observe that the cosmological work done can be eval¬ 
uated if we know the explicit form of the variation of the scale factor with time, our 
main motivation is to hnd an explicit expression for the scale factor in terms of time 
‘T. For this we need to solve Eq. (4.9), Eq. (2.5) and Eq. (2.7) consistently. There¬ 
fore, we will basically substitute the expression for the Hubble parameter from the 
Friedmann equation into Eq. (2.7) and replace p by Eq. (2.5) with a specihed form 
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of the cosmological potential. Then we can solve for the scalar held and substitute it 
back to the expression as appearing in the Friedmann equation to get an expression 
for scale factor in terms of ‘f'. For the sake of simplicity henceforth we will only 
concentrate for the case of A; = 0, which is also a valid and natural assumption in the 
present context, since we know that our present observations predict a nearly hat 
universe. In order to simplify the analysis further, we will consider the case when 
the contribution from the kinetic term is lesser than the potential energy i.e. 

02 « 1/(0) (4.65) 

during the expansion and similarly in the physical situation where the kinetic term 
is larger than the potential energy i.e. 

02 » 1/(0) (4.66) 


during the contraction phase of the Universe. Hence the approximate forms of 
Eq. (2.5) and Eq. (2.7) during expansion and contraction phase of the Universe 
will be given by: 


3770+ — ^0, 
0 + 3770 0, 


p ~ U(0) (during expansion) 

0 ^ . . 

p ^ — (during contraction) 


(4.67) 

(4.68) 


Next step is to specify the specihc form of the potential, in order to get further 
informations and the constraints from the equations. In the next subsections we 
have explicitly shown the analysis for three different potentials i.e. Hilltop, Natural 
and origanted Coleman-Weinberg potential. 


4.6.1 Case I: Hilltop potential 


In case of hilltop models the potential can be represented by the following functional 
form [47]: 

r / A \ p-\ 


^(0) = Co 


1 + 0 



(4.69) 


where Vq = is the tunable energy scale and 0 is the index which characterizes the 
feature of the potential. In principle 0 can be both positive and negative. Addition¬ 
ally it is important to note that, in the present context, Vq mimics the role of vacuum 
energy. Since our present job is to substitute the expression for the scale factor into 
Eq. (4.64), we need to hnd separate expressions for the scale factor for both early 
and late times during the expansion and contraction phases of the Universe. 


A. Expansion 
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i) Early time; 


When the scale factor is lying within the window, amin < a < a' or eqnivalently 
when pr'l/Ml » 1, we can use the approximated version of the Friedmann equation 
as given by Eq. (4.14), with fc = 0, where the energy density of the scalar field p 
is for now described by the hilltop potentials. But instead of considering only the 
zeroth order term (as we had done in order to find the condition for bounce), here 
we will consider upto the first order terms (as was done for finding an expression for 
cosmological work done). Then substituting the resulting expression for the Hubble 
parameter H into Eq. (4.67), we get the following integral equation in DGP brane 
world as: 


/ 


d 





(4.70) 


The exact solutions of the above integral equation is given in the Appendix, which 
we see has a very complicated form for Randall Sundrum (RSII) limiting situation. 
Hence to simplify the integrals, we use the following redefinition of the field variables: 

= (4.71) 

where now we will solve for A. In order to further simplify the expressions, we solve 
for two limiting cases: 


a) 0 /M 4 << 1 : 


For this case we can expand the exponentials upto linear order and then using the 
result the integral on the left hand side of Eq. (4.70) becomes: 

Ml) f^P 


1 + 


(3pX 


1 

2(i + /3); 


(l+A)[l-(p-l)A]ciA = 


Rn 


t+c 


which is true under the assumption that the quantity 

/3Xp 

<< 1 


3M| 

(4.72) 

(4.73) 


(1 + /3) 

is small. It is important to note that in the present context c is an arbitrary integra¬ 
tion constant. 

Hence solving the above integral on the left hand side, we get the expression for 
A, hence for (p/M^ = as: 


A = 


Tn? 


7 Vn A (1+/3)1/" I 1 

/9p 2rc 


(ti — t) Xi 


(4.74) 


\3MiJ 
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Here Aj is the value of A at the initial time scale t = ti. 

Substituting the expression for A back into the Friedmann equation, we get the 
expression for scale factor as 


a{t) = tti exp 


3Mi) 


+ /3p 


Vo^ 

■iMt 


1/2 


2 ( 1 +/?)V 2 



F 

3M2 

11 

r‘ (^: 

.1/2(i^/3)1/2 ^ 

/3p 2rc ^ 

ij 


3M| 


/ ( _fb 

Pp 

\ 


1/2 




2(1 + /3)V2 


Vo 


2 ( F) A ^ (i+/3)^/^ I_ i_ 

j /Jp ^ 2 r, 




where a* is the value of the scale factor at time scale t = ti which is in our case hxed 
at the bouncing time scale ti = tb- 


b) 0 /M 4 >> 1 : 


For this physical situation Eq. (4.70) simplifies to the following expression; 




Hn 


3 M|) 


1/2 


oP '^/2 _|_ 
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Hn 
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(4.76) 


where we use Eq. (4.71). Here we have also assumed that » 1, which is a 
valid assumption since we are working in the limit where 0 /M 4 >> 1 . 

Hence solving the above integral we get the expression for A as 


A = In 
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Here Aq is an arbitrary integration constant which can be expressed as: 

1/2 
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Here we have neglected l/2rc as the following condition holds good: 
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— << 
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Vo 
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(4.77) 


(4.78) 


(4.79) 


which is a valid assumption since we are studying the cosmological consequences in 
the context of early universe in this limiting case. 
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Further substituting back the expression for A into the Friedmann equation and 
integrating, we get the following expression for the scale factor as: 


a(t) = tti exp 


/3p(p-4) ((2- 2) (Ao - 

4 


P ■ 
4 —p 


(4.80) 


ii) Late time: 


When the scale factor is lying within the window, a' < a < a^ax or equivalently 
when pr‘l/M‘1 « 1, we can use the approximate Friedmann equation given by 
Eq. (3(a)) with A: = 0, where p is now characterized by the hilltop potential. Then 
substituting the resulting expression for the Hubble parameter H into Eq. (4.67), we 
get an integral equation of the following form; 


/ 


d 




P\ 1/2 



(4.81) 


The exact solution of the above integral is given in the Appendix for the Randall 
Sundrum (RSII) model. For the sake of simplicity, we follow the same procedure as 
we have done in the previous case i.e. first use the field redefinition 0 /M 4 = and 
then study two limiting cases. 


« 1 : 


Once again expanding the exponentials upto linear order we get 
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A =-- 
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t + Ai 


(4.82) 


where Ai is the arbitrary integration constant. The expression on the left hand side 
has been obtained under the underlying assumption that the following condition; 


2U0/3 

3 Mi 
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o I 2 Vb+Vb/3 I 1 
^ 1 3 —^ r ::2 


Mf 


A << 1 


(4.83) 


is satisfied. As we performing the analysis in late time in the present context, hence 
1 /r^ is a large quantity, but to strictly satisfy this condition, we need to choose Vq 
in such a way that the following condition: 


2 Vo 1 2 Vo/3 

3M! ^ ^ ^ 3 


» 1 


(4.84) 
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is always satisfied during the late time. 

Further solving the integral as stated in Eq. (4.82), we get the following simplihed 
expression for A as: 


(^1 3M|^) 

( 2 Vo+VoP , ' 

Ml + rl) 


(4.85) 


where Ai is given by: 
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Vo 

3M| 


if 


(4.86) 


where A/ is the value of A at turnaround corresponding to time scale t = tf. 

In order to get the expression of the scale factor, we substitute this expressions 
back into the Friedmann equation and integrate the equation to get the following 
result: 


a{t) = A!^ exp 


where 
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(4.87) 


A!^ = af exp 
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Here a/ is the value of the scale factor at the time of turnaround. 
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(4.88) 


B. Contraction 

From Eq. (4.68), we observe that for the contraction phase of the Universe, un¬ 
der the approximation which we have assumed, the analysis becomes independent of 
the specihc form of the cosmological potential. Hence, the results for the DGP brane 
world model hold for any form of the potential. 


i) Early time: 


When the scale factor is lying within the window, a^m < a < a', following the same 
analysis as we did for the case of expansion, with the new expression for density, the 
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integral equation is now expressed in the following form: 

f I ■ f 

—d(j) = — dt 




(4.89) 


^ \/671^4 27* c 

Solving the above integral equation we get the following expression for cj) as: 


2rr. 


exp 


d-M) 

2rc 


(4.90) 


where A 2 is the arbitrary integration constant given by: 

y/GMi 


A 2 = ti- 2rc In 




+ 1 


(4.91) 


Here 0* is the value of the derivative of the scalar held at time of bounce. 

Further integrating the above equation once again, we get the expression for the 
held 0 as: 


0(f) = 




In ( 1 — 2e" 


12 , 


Tr-t 


where X 2 is the arbitrary integration constant given by: 


K =ti-rc hr 


- 1 - exp 


0i^c 

y/GM4 


A ti ] — 


(4.92) 


(4.93) 


where 0* is the value of the held at the time of bounce t = ti. Hence substituting 
Eq. (4.92) in the Friedmann equation, we get the following solution for the scale 
factor as: 

(4.94) 


a{t) = A 2 ^e 2 ’'c 


^2 


where A 2 is the arbitrary integration constant given by: 


A'' = 


di 


t ^2 

g 2rc — g 2rc 


(4.95) 


ii) Late time: 


When the scale factor is lying within the window, a' < a < amaxi using the Fried¬ 
mann equation given by Eq. (3(a)) and the late time condition, we get the following 
solutions from DGP brane world model: 
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(4.96) 

(4.97) 
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where ^3 and A 4 are the arbitrary integration constants given by: 


r 

A 3 = - In 
^ 3 

1 _ 1 

+ 

Tr 

STc 


= 


0 j 


(4.98) 

(4.99) 


Here 0/ and 0/ are the values of the scalar field and its derivative at turnaround 
time scale t = tf. 

Further substituting the expression for 0 in the Friedmann equation, hnally we 
get the expression for the scale factor as: 


a(t) 


H 5 exp 




4 A 


(Ai-t) 


— tanh ^ 
Tc 


where is the arbitrary integration constant given by: 







(4.100) 
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(4.101) 


C. Expression for work done 

Using the solutions for the scale factor, one can further compute the expression 
for the integrals as appearing in Eq. (4.64), giving the work done in one cycle for the 
cases for which we have the hnal expressions for the scale factor. 

For the case 0 /M 4 << 1, using the solutions of the scale factor as computed in 
the earlier section for hilltop potential, we hnd that the expressions for work done 
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for spatially flat case /c = 0 given by the following expression; 
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(4.102) 
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(4.104) 

(4.105) 


and consequently the total work done in a one cycle can be expressed as: 

pdV = -3Mla'2 - 18ai(e“i‘™" - e“i‘') 

+ 902+ 6a3(t„2in - t') 

-3Mlai{Ei{[a5{2a'jtmin - ae)] - Erf[a5(2a7f' - oe)]) (4.106) 

where ai...a 7 are constants that depends on the model parameters present in the 
expressions for the scale factor whose explicit expressions have been given in the 
appendix. Here we have quoted the results corresponding to each integral of the 
work done as given in Eq. (4.64). Thus, we see that we get 

pdV ^ 0, (4.107) 


whose signature depends on the numerical values of the constants. Thus, one can 
conclude that the phenomenon of hysteresis is fruitfully achieved for small held hill¬ 
top potentials. 


4.6.2 Case II: Natural potential 


In case of natural models the potential can be represented by the following functional 
form [70]; 


V{cl>) 


Eo 


1 -|- cos 



(4.108) 
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where Vq = is the tunable energy scale and / plays the mass scale of the pot. It 
is important to note that, the potential is periodic or equivalently shift symmetric 
for the shift in the field coordinate: 


0 -)■ 0 + 27r/, 


(4.109) 


where 27r/ mimics the role of angular shift in the field coordinate. Since we need to 
substitute the expression for the scale factor into Eq. (4.64), we need to hnd separate 
expressions for the scale factor for both early and late times during expansion and 
contraction phases. 

A. Expansion 

i) Early time: 

When the scale factor is lying within the window, Omm < a < a' or equivalently 
when pr‘l/M‘1 » 1, we can use the approximate Friedmann equation given by 
Eq. (4.14) with spatially flat case fc = 0, where p is now given by the potential for 
natural potential. But instead of considering only the zeroth order term (as we had 
done in order to hnd the condition for bounce), here we will consider upto hrst order 
terms to compute the expression for work done. Then substituting the resulting 
expression f or Hubble parameter H into Eq. (4.67), we get an integral equation of 
the following form; 



(4.110) 


The exact solutions of the above integrals are given in the Appendix for Radall 


Sundrum single brane world model (RSII). Simplihed analytical expressions could 
be obtained only for the small held case {(j)/f « 1) which has been discussed 
below. 


For 0// << 1 case we take small argument approximations of the trigonometric 
functions after which we get the following solution for the sclar held 0 as: 



(4.111) 


where Aq is the arbitrary integration constant is given by: 



(4.112) 
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Here (f)i is the value of 0 at the initial time of bouncing time scale t = ti. 

Further substituting the expression for 0 back into the Friedmann equation, we 
get the following simplihed expression for scale factor as: 



(4.113) 


where A-/ is the arbitrary integration constant is given by: 



(4.114) 


ii) Late time: 

When the scale factor is lying within the window, a' < a < a^ax oi' equivalently 
when prl/Mj « 1, we can use the approximated form of the Friedmann equation 
given by Eq. (3(a)) with spatially flat case A; = 0. Hence substituting the resulting 
expression for the Hubble parameter H into Eq. (4.67), we get the following integral 
equation as given by: 



(4.115) 


The exact solution of the above integral is given in the Appendix for Radall Sundrum 
single brane world model (RSH). Therefore, in order to simplify the analysis, we 
again study small held limiting case as mentioned earlier. For 0// << 1 case we 
take small argument approximations of the trigonometric functions and hnally we 
get the following simplihed expression for the scalar held 0 as: 



(4.116) 


where Ag is the arbitrary integration constant is given by: 



(4.117) 


Here is the value of the scalar held at turnaround time scale t = tp- 


36 






















In order to get the expression of the scale factor, we substitute this expression 
back into the Friedmann equation and integrate the equation to get following sim¬ 
plified expression for the scale factor as: 


a{t) = Agexp 



where Ag is the arbitrary integration constant is given by: 


Ag = apexp 



(4.118) 


(4.119) 


where ap is the value of the scale factor at turnaround time scale t = tp. 


B. Contraction 


As has already been mentioned while performing the analysis for expansion, the 
conclusions for contraction phase is independent of any potential, hence the analysis 
remains same for natural potential also. 


C. Expression for work done 


Here we also study the work done for all the cases for which we have expressions 
for the scale factor for all the integrals appearing in Eq. (4.64). We consider the 
case 0 // << 1 for which the expression for work done for spatially flat case /c = 0 is 
given by: 
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and consequently the total work done in a one cycle can be expressed as: 


pdV = 0 — 3Mla'2 


.2e3ap™„ + - 18 ai(e“i* 


t 


9a2(e^“'h™" - e^“'h') + 6a3(tmin - i') - q 

(4.124) 


Here ai..b 2 are all constants dependent on the model parameters that appearing 
in the expressions for the scale factor. Their explicit forms are in the appendix. The 
results obtained in this section implies that, we get 


pdV 7 ^ 0 


(4.125) 


for small held limit 0 // << 1 and the signature of the integral is completely deter¬ 
mined by the numerical values of the model parameters and the arbitrary integration 
constants. Most importantly, the analysis shows that the phenomenon of hysteresis 
is fruitfully achieved for natural potential. 


4.6.3 Case III: Coleman-Weinberg potential 

Let us start with a theory of hve dimensional M = 2 bulk supergravity in which by 
compactifying the extra hfth dimension it is possible to derive an four dimensional 
effective theory described by TV = 1 supergravity theory in brane world. Within this 
prescription, the four dimensional one-loop effective Coleman-Wienberg potential 
embedded in the brane world can be expressed as [71-73]: 


VW) 


v„ 


1 + 





(4.126) 


where Vq sets the energy scale of supergravity theory. Additionally the model param¬ 
eter a signihes the tree level effect and the parameter (3 characterizes the effect of 
one-loop correction to the leading order result. Here M 4 represents the background 
mass-scale of theory. For sake of simplicity one can consider M 4 to be the UV cut-off 
i.e. the Planck scale of the gravity theory. 


A. Expansion 
i) Early time: 


When the scale factor is lying within the window, Omm < a < a' or equivalently 
when pr‘l/M‘1 » 1, we can use the approximated form of Friedmann equation given 
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by Eq. (4.14) with spatially flat case /c = 0, where in the present context p is given 
by the snpergravity motivated potential. Instead of considering only the zeroth order 
term, here we will consider upto first order terms. Then substituting the resulting 
expression for the Hubble parameter H into Eq. (4.67), we get an integral equation 
of the following form; 
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dt (4.127) 


For the sake of simplicity let us consider the following transformation or redefinition 
in the held; 


Ma 


= e 


(4.128) 


as we had done for the previous cases, where now we will solve for the redefined held 
A. In order to further simplify the expressions, we solve for two limiting cases: 


a) 0 /M 4 << 1 : 


For this case we can expand the exponentials upto linear order after which the inte¬ 
gral on the left hand side of Eq. (4.127) becomes 
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Next we compute the above integral equation on the left hand side under the as¬ 
sumption that the values of model parameters, a and (3 satisfy the constraint: 
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Consequently we get the expression for A, hence 0 /M 4 = as: 
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Here Aj is the value of A at the initial time f of bounce. 

Substituting the expression for A back into the Friedmann equation, we get the 
expression for scale factor as: 


a{t) = Aiiexp 




\3MiJ 

where An is the arbitrary integration constant given by: 


An = tticxp 
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1/2 
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Here a, is the value of the scale factor at the time of bounce. 


(4.134) 


(4.135) 


b) (jy/Mi » 1: 


For this case Eq. (4.127) simplifies to 
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(4.136) 


The above expression is obtained provided that the following constraints are satisfied: 
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We have also used the condition pr^/M| >> 1, which is true in the case of early 
universe. Here An is an arbitrary integration constant. 

Solving the above integral equation we get the following expression for A as: 
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A = 



where the explicit form of An is given by: 
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Here A* is the value at the time of bounce t = ti. 

Further substituting back the expression for A into the Friedmann equation, we 
get the following expression for the scale factor as: 
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a(t) = Ai 3 exp 



(4.141) 


where ^13 is the arbitrary integration constant given by: 



(4.142) 


ii) Late time: 


When the scale factor is lying within the window, a' < a < a^ax oi' eqnivalently 
when pr‘l/M‘1 « 1, we can use the approximated version of the Friedmann equation 
given by Eq. (3(a)) with spatially flat case fc = 0. Hence substituting the resulting 
expression for the Hubble parameter H into Eq. (4.67), we get the following integral 
equation given by: 
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For the sake of simplicity during the analysis, we follow the same procedure as 
mentioned before and use the held redehnition: 



and then study two limiting physical situations. 


(4.144) 
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a) (jy/M^ « 1: 


For this small field limit expanding the exponentials upto linear order we get the 
following integral equation; 
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where the integral in the left hand side has been obtained under assumption that: 
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3 + ct + 4 q;A + (dX) « 1 , 
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which is a valid assumption since we are in the late universe and in the limit of small 

A. 

Solving the above integral, we get the expression for A as 
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where A/ is the value at turnaround corresponding to time t = tf. 

In order to get the expression of the scale factor, we substitute this expression 
back into the Friedmann equation and integrate the equation to get: 
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where Ai^ is the arbitrary integration constant given by: 
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b) (j) » M 4 : 


For this case the integral in Eq. (4.143), writing 0 /M 4 = e^, simplifies to: 
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But in order to get an solution for A, we need to further simplify the integral, which 
is possible if we assume that the values of the parameters of the model satisfy the 
condition; 
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is satisfied. This is possible because though we are in the late time i.e. 
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we have both the above conditions being satisfied by 1/r^ simultaneously. Under 
this assumption, we get the solution for A as: 
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where Aiq is the arbitrary integration constant given by: 
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Substituting the expression for A back into the Friedmann equation, we get the 
expression for the scale factor as: 


a{t) = Ai 7 exp 
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where An is the arbitrary integration consant given by: 
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Here, aj is the value of the scale factor at turnaround t = tf. 

B. Contraction 

As has already been mentioned that while performing the analysis for expansion, 
the conclusions for contraction phase is independent of any potential, hence rest of 
the analysis remains same as mentioned earlier. 

C. Expression for work done 

The expression for work done in this case is same as we get for the case of natu¬ 
ral potential with the expressions for constants now given by the parameters of this 
supergravity motivated model. Hence we see that in this case also, applying certain 
approximations and physical limits, we can get an analytical expression for work 
done which is non zero, thus leading to the phenomenon of cosmological hysteresis 
in the present context. 


4.7 Graphical Analysis 
4.7.1 Case I: Hilltop potential 
Graphical Analysis: 


All the graphs in this section and in the following sections have been plotted in 
units of Mp = 1, Hq = c=l, where Mp is the Planck mass, Hq is the present 
value of the Hubble parameter and c is the speed of light. Throughout the analysis 
Tc will be expressed in units of hence only its magnitude will be written explicitly 

In Fig. 4 we have shown the evolution of the scale factor and the potential for 
early and late time expansion phase. From the above plots, we can draw the follow¬ 
ing conclusions; 

• Fig. 4(a), shows the plot of the scale factor in the small held limit for hilltop 

potential given by Eqn. (4.75) with the parameter values Vq = p = 

2, /3 = 0.001, Ai = 10-®, U = 0.1, tti = 1, = 1. 

• Higher values of the integration constants only results in an increase in the 
amplitude of the scale factor, not affecting the nature of the graph. 

• From Eqn. (4.75), we can conclude that solution will only be possible for 
(3 > 0. Changing the value of f3 mildly modihes the amplitude of expansion. 
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10 ® 

10 ® 


a vs t for DGP model (early time) 


V(i/i(t)) vs t(i/i«IWp) for DGP model 




tpnMp-') t(inMp-') 

(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during the early expansion phase tial during early expansion phase for (p << Mp 
for (p « Mp with Vq = p = 2, j3 = with Vq = 4.1xlO“^Mp, P = 0.72, p = 2, Xi = 

0.001, X, = 10“®, ti = 0.1, Tc = 1. 6.6, ti = = 1, Tc = 1 . 

a vs t for DGP model (late time) V(0(t)) vs t (0«iWp) for DGP model 




(c) An illustration of the behavior of the scale (d) An illustration of the behavior of the poten- 
factor with time during late time expansion phase tial during late time expansion phase for p << 
for (j) « Mp with Vq = 10“®Mp, p = 3, /3 = Mp with Vq = 3.7xlO“^A/p, p = 3, P = 
0.001, Ai = 10-®Mp, A[ = 10-®, rc = 1.0. 0.14, Ai = 20.8Mp, = 1.0 . 

Figure 4. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for DGP model. 

but the nature of the graph remains same. The plot is almost independent of 
any variation in p. Larger values of Vc increases the amplitude of the graph. 
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• Fig. 4(b) shows the plot of the behavior of the potential with time for small held 
hilltop potential. This graph has been obtained with the help of Eqn. (4.74) 
with parameter values Vb = 4.1xlO“^Mp, /3 = 0.72, p = 2, A* = 6.6, U = 
95Mp"^, Tc = 1. The evolution of the potential is not much affected by the 
value of p. Large values of Vc increases the height of the potential and makes 
the potential fall less steeply. 

• Fig. 4(c) shows the plot of the scale factor in the late time expansion phase 

for hilltop potential given by Eqn. (4.87). This plot has been obtained for 
Vo = p = 3, f3 = 0.001, All = 10-®Mp, A[ = 10“®, = 1.0. Detail 

graphical analysis have shown that the expansion of the universe is almost 
independent of the values that /3, p, Ai takes i.e in this case both negative and 
positive values of (3 can give rise to an expanding universe. We do not get the 
required expansion for a value of Tc < 0.04. 

• Fig. 4(d) shows the evolution of the potential in the late time expansion phase 
for hilltop potential. Fig. 4(d) has been obtained by with the help of Eqn. 
(4.85 with the parameter values Vq = 3.7xlO“^Mp, p = 3, /3 = 0.14, Ai = 
20.8Mp, Tc = 1.0. Detail graphical analysis show that the larger the value of (3, 
smaller is the the range of Vc for which we get potential giving rise to expansion. 
The graph is almost independent of the values that p, Ai takes. 

• If we compare Fig. 4(b) and 4(d), we hnd that the potential falls more steeply 
during the late time than in early time. This is expected because the potential is 
near its end of expansion phase, hence kinetic term starts dominating resulting 
in steeper fall of the potential. 

In Fig. 5 we have shown the evolution of the scale factor and the potential for early 
and late time contraction phase. From the above plots, we can draw the following 
conclusions: 

• Fig. 24(a), shows the plot of the scale factor during late time contraction phase 
for hilltop potential given by Eq. (4.100) with the parameter values ^4 = 
mM-\ Tc = 6 . 

• The nature of the graph (i.e the decrease in amplitude of the scale factor with 
time, which is expected for contraction phase) remains almost unchanged if 
we change the parameter values. Larger values of r decreases the amplitude of 
expansion (but the change is very small), but smaller values make the potential 
fall more steeply. 

• As we can see from Fig. 24(a), the change in amplitude of the scale factor is 
very small. This change becomes negligibly small for very small values of A 4 , 
and very large values of A 4 are also not allowed. 
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(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the po- 
factor with time during the late contraction phase tential during late contraction phase with Vq = 
with A 4 = 100M-\ rc = 6 . A 3 = 152M-\ A 4 = 730Mp, P = 

0.52, rc = 1.62, p=l . 

a vs t for DGP model (early time) V(0(t)) vs t for DGP model 




(c) An illustration of the behavior of the scale fac- (d) An illustration of the behavior of the po- 
tor with time during early time contraction phase tential during early time contraction phase with 
with A 2 = lOOM-i, Ag = 1, \rc 1= 2.0. Vq = p = 1 , p = 0.392, A^ = 63M-^, \ 

Tc 1= 1-4 . 

Figure 5. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for DGP model. 
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• Fig. 24(b) shows the plot of the behavior of the potential with time for late time 

contraction of hilltop potential. This graph has been obtained with the help 
of Eq. (4.97) with parameter values Vq = 6.1xlO“^Mp, A 3 = A 4 = 

730Mp, f3 = 0.52, = 1.62, p = 1. 

• The potential in this case rises with time,which is expected since we are in the 
contraction phase and kinetic energy is dominating initially. Detail graphical 
analysis show that for very small value of Tc, the correct nature of the potential 
is obtained only if we take smaller values of the other parameters. Larger values 
of the parameters make the potential rise more steeply and attained the nearly 
flat region faster. 

• Fig. 5(c) shows the plot of the scale factor in the early time contraction phase 

for hilltop potential given by Eq. (4.94). This plot has been obtained for 
A 2 = lOOMp"^, A!^ = 1, I 1= 2.0. Detail graphical analysis have shown that 
as we decrease the value of | |, the fall of the scale factor becomes more 

linear. Larger values of | | increases the amplitude of the scale factor but 

makes contraction possible only for smaller values of A 2 . 

• Fig. 5(d) shows the evolution of the potential in the early time contraction 

phase for hilltop potential. Fig. 5(d) has been obtained by with the help of 
Eq. (4.92) with the parameter values Vq = p = 1, /? = 0.392, A '2 = 

63Mp"^, I Tc 1= 1.4. Detail graphical analysis show that the nature of the graph 
do not depend on the parameter values. Larger values of the parameters only 
increases or decrease the height of the potential. But in order to get the correct 
nature of the potential which will result in contraction, we need /? > 0. 

• Though fig. 24(a) and Fig. 24(b) are for late time and Fig. 5(c) and Fig. 5(d) 
are for early time, in one complete cycle, the late time phase appears before 
the early time phase in contraction under the convention which we follow. But 
in this case, the time range has been chosen not according to the convention, 
but in order to get some suitable output due to the complicated nature of the 
equations. 

• If we compare Fig. 24(a) and Fig. 5(c), we find that the decrease in the ampli¬ 
tude of the scale factor is much less during the late time contraction phase, as 
compared to early time phase. 

• If we compare Fig. 5(d) with Fig. 4(a), we find that there is a net increase in 
the amplitude of the scale factor after one complete cycle. 
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a vs t for DGP model (early time) 



t(inMp-') 


V(0(t)) vs t (0«f) for DGP model 



(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during the early expansion phase tial during early expansion phase for (j) « f with 
for 4'« f with Vb = 10-®M4 rc = 1, A 7 = 1- Vo = Vc = 1.67, Aq = 32, / = 


6 

a VS t for DGP model (late time) 





.7Mp . 


V(^(t)) VS t(^«f) for DGP model 



t(inMp-') 


(c) An illustration of the behavior of the scale (d) An illustration of the behavior of the poten- 
factor with time during late time expansion phase tial during late time expansion phase for (j) << f 
for (/) << / with Vq = 10“®Mp, rc = 1.17, Ag = with Vq = 4.7xlO“^Mp, rc = 1.67, Ag = 1, / = 

6.7M„ . 


1 . 


Figure 6. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the DGP model. 
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4.7.2 Case II: Natural potential 


In Fig. 6 we have shown the evolution of the scale factor and the potential for early 
and late time expansion phase. From the above plots, we can draw the following 
conclusions: 

• Fig. 6(a), shows the plot of the scale factor in the small field limit for natural 

potential given by Eq. (4.113) with the parameter values Vq = 10“®Mp, = 

1, At = 1 during the early phase of expansion. 

• Higher values of the parameters only increase the amplitude of expansion, keep¬ 
ing the nature of the plot unchanged. 

• Fig. 6(b) shows the plot of the behavior of the potential with time for small 
field natural potential during the early phase of expansion. This graph has 
been obtained with the help of Eq. (4.111) with parameter values Vq = 
4.7xlO“^Mp, Tc = 1.67, Aq = 32, / = 6.7Mp. Large values of Tc and / 
increases the height of the potential. Only for certain range of values of Aq, we 
get the required evolution of the potential (such as expansion is possible if Aq 
lies within (1 to 3) or (6 to 9) or (13 to 15) etc.). 

• Fig. 6(c), shows the plot of the scale factor in the small field limit for natural 
potential given by Eq. (4.118) with the parameter values Vq = 10“®Mp, Vc = 
1.17, Hg = 1 during the late phase of expansion. 

• Higher values of the parameters only increase the amplitude of expansion, keep¬ 
ing the nature of the plot unchanged. 

• Fig. 6(d) shows the plot of the behavior of the potential with time for small field 
natural potential. This graph has been obtained with the help of Eq. (4.116) 
with parameter values Vq = 4.7xlO“^Mp, Vc = 1.67, Hg = 1, / = 6.7Mp. 
The conclusions regarding the variation of the nature of the potential with 
parameters remains same as for the early time expansion. 

• If we compare the amplitudes of the scale factor in Fig. 24(a) with Fig. 6(a), 
we find that after one cycle of expansion and contraction, we can get a net 
increase in amplitude of the scale factor provided the parameters are chosen 
properly. 

4.7.3 Case III: Coleman-Weinberg potential 

In Fig. 7 we have shown the evolution of the scale factor and the potential for early 
and late time small field expansion phase for supergravity potential. From the above 
plots, we can draw the following conclusions: 
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a vs t for DGP model (early time) 


V(^(t)) vs t {^«Mp] for DGP model 




tlinlWp-') t(inMp-') 

(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during the early expansion phase tial during early expansion phase for (j) << Mp 
for (j) « Mp with Vq = Vc = with Vq = 10~^Mp, rc = 0.86, Aiq = 5Mp, a = 

0.86, All = 1, a = 0.1. 0.02, /3 = 0.03 . 


a VS t for DGP model (late time) 



V(0(t)) vs t (0«Mp)for DGP model 





(c) An illustration of the behavior of the scale fac- (d) An illustration of the behavior of the poten- 
tor with time during late time expansion phase tial during late time expansion phase for (j) « 


for (j) « 
1.66, Ai4 = 10 


Mp with 


l/o = 6 x 10-4M4, r, 
Ai 5 = 1, a = 0.145. 


= Mp with 


Vo = 10-5M4, re = 1.32, Am = 


10-®Mp, a = 0.032, ^ = -1.9 . 


Figure 7. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the DGP model. 
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• Fig. 7(a), shows the plot of the scale factor in the small held limit for super¬ 
gravity potential given by Eqn. (4.135) with the parameter values 

Vo = Tc = 0.86, All = 1, a = 0.1 

during the early phase of expansion. 

• Higher values of the parameters only increase the amplitude of expansion, keep¬ 
ing the nature of the plot unchanged. The variation in the nature and amplitude 
of the plot is almost independent of the value of a parameter. 

• Fig. 7(b) shows the plot of the behavior of the potential with time for small 
held supergravity potential during the early phase of expansion. This graph 
has been obtained with the help of Eqn. (4.132) with parameter values Vq = 
10“®Mp, Tc = 0.86, Aio = 5Mp, a = 0.02, [3 = 0.03. Large values of the 
parameters causes an increase in the amplitude of the expansion, keeping the 
nature of the plot unchanged. Largest increase in amplitude occurs for any 
change in value of the parameter a. Detail graphical analysis have shown that 
we do not get the required nature of the plot for (3 < 0. Thus only positive 
values of (3 are allowed in this case. 

• Fig. 7(c), shows the plot of the scale factor in the small field limit for su¬ 
pergravity potential given by Eqn. (4.149) with the parameter values Vq = 
6 xlO“^Mp, Tc = 1.66, Ai 4 = 10“®Mp, A 15 = 1, a = 0.145 during the late 
phase of expansion. 

• For Vo > 6x10“*^Mp, expansion is possible only if the values of /3, a and Tc lie 
close to or lesser than unity. If we make the values of any of the parameters 
/?, a, Tc large, then expansion is possible only for the case when both A 14 and 
Vo is 0 ( 10 “^) or less. 

• Large and positive values of f3 make the expansion more linear. Expansion is 
possible for /3 < 0 , provided Vq takes large values and a and A 14 are smaller 
than unity. 

• Fig. 7(d) shows the plot of the behavior of the potential with time for small 

held supergravity potential during late time expansion. This graph has been ob¬ 
tained with the help of Eqn. (4.147) with parameter values Vq = 10“®Mp, Tc = 
1.32, Ai 4 = a = 0.032, j3 = —1.9. Expansion is possible only li j3 <2. 

Again for expansion to happen for positive value of /3, all the other parame¬ 
ters must be much less than unity. Larger and positive values of (3 makes the 
potential fall more linearly. For expansion to occur, A 14 must be << 1. 
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• If we compare the amplitudes of the scale factor in Fig. 24(a) with Fig. 7(a), 
we find that after one cycle of expansion and contraction, we can get a net 
increase in amplitude of the scale factor provided the parameters are chosen 
properly. 

For large field case {4> » Mp) during expansion, the conclusions and the nature 
of the graphs are almost similar, hence have not been shown explicitly. Thus, a 
increase in the amplitude of the scale factor after one complete cycle of expansion 
and contraction, is possible for this case also. 

5 Hysteresis from cosmological constant dominated Einstein 
gravity model 

In this section we will consider three different cases by taking different functional 
forms of the cosmological constant like term in the effective action. One will be the 
case of standard model of cosmology i.e when the cosmological constant is a constant. 
Then we will further extend our analysis for the case when the cosmological constant 
is field dependent. This dependence will be included either in the form of a power 
series form of the scalar field or by including a dilaton field. 

Case 1: Cosmological constant of ACDM model 


Let us start our discussion with ACDM model, which is the standard model of Big 
Bang cosmology, where the Einstein’s field equations are modified by the additional 
cosmological constant term denoted as A. The modified Friedmann equations in this 
model can be expressed as: 



(5.1) 



(5.2) 


where we are denoting the Planck mass by M. 

Case 2: Scalar field dependent cosmological constant 

Let us consider a situation where we make the cosmological constant term scalar 
field dependent i.e we replace A by A(0), where 


4 



(5.3) 


i=0 
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then the representative four dimensional effective action for this case is given by: 



where Aq is the general cosmological constant (denoted by A in case 1) and 
Ai, A 2 , A 3 &: A 4 are the new constants. 

The hrst action S/<^eh is the normal action for ACDM model, hence the Fried¬ 
mann equations for this part are same as Eqns. (5.1) and (5.2). The scalar held 
dynamics get modihed due to the action Sa^. The equation of motion for the scalar 
held is given by: 


5d^^(j) 5(j) 


(5.5) 


where Caip is the Lagrangian density for the scalar held and for FLRW metric \/—g = 
a^. Solving the above equation we get the equation of motion for (j) as: 


4 

0 + 3i70 - ^ + 1/^ = 0 (5.6) 

i=l 


where the hrst two terms come from the hrst term in Eq. (5.5) and the rest of the 
terms come from the last term in Eq. (5.5). 

The energy momentum tensor for scalar held is given by 


= di,(j)d^(j) - g^,yCA4> (5.7) 

= d^(f)dy(i) - g^y (i)d„(i) - 1/(0) ^ Ai0* 

For perfect huid the expressions for the energy density p and pressure p can be 
expressed as: 

1 . ^ . 

P = 2*^^ + y{(l>) - 

i=l 

1 , 4 

i=l 

Solving Eq. (5.6), Eq. (5.9) and Eq. (5.1) simultaneously, we get the solutions for the 
scalar held 0 and scale factor a. Thus we see that both the density, pressure as well 
as the equation of motion of the scalar held gets modihed from the standard case. 


(5.8) 

(5.9) 
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which was expected. 


Case 3: Cosmological constant along with a dilaton field 


In this case we add a dilaton field along with the general cosmological constant 
i.e we replace A by A(0), where 


A(0) = Ao + Ae-^/^A 

Then the representative action for this case is then given by: 

S = j (fx^g + Ao + - C(0)^ 

= j d‘^x^/^ {R + Ao) + j 00^.0 


(5.10) 


(5.11) 


^ 0 ) , 


Saeh 


Sa4i 


where Aq is the general cosmological constant and A is another constant. 

Repeating the analysis as we did for the previous case 2, we get the following 
equation of motion for the scalar field 0 as: 

0 + 3id0 - + 1/^ = 0 (5.12) 

IVlp 

Similarly the expressions for density and pressure for a perfect fluid in this case 
can be expressed as: 

p = ^02 + 1/(0) - (5.13) 

p = ^02 - 1/(0) + (5.14) 

Solving Eq. (5.12), Eq. (5.14) and Eq. (5.1) simultaneously, we get the solutions for 
the scalar field 0 and scalae factor a. 


5.1 Condition for bounce 
For case 1, 2, 3 


In the analysis mentioned below, we have first shown the results for open and closed 
universe only. The analysis for fiat universe has been shown separately. 

Let us start the discussion for the case 1 where at bounce, setting the Hubble 
parameter id = 0 in Eq. (5.1) we get: 


pb = 3M^ 




(5.15) 
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where pi, and at are the representative density and scale factor at bonnce. 

The mass content at bonnce (neglecting the overall constant factor) is given by: 


Mb = PbCt-b = “iM I kttb — 


Aal 


(5.16) 


(5.17) 


(5.18) 


Therefore the inhnitesimal change in the mass content can be expressed as: 

SM = 3M^ {k — Aal) Sab- 
Now using energy conservation one can write: 

6 M + 6W = 0, 

where SW is the work done during each expansion-contraction cycle. 

Further setting 

6 M =-6W = - j) pdV, (5.19) 

we get the expression for change in amplitude of the scale factor at each successive 
cycle in terms of the work done as: 


S a^irt 


3M2 {k - Aa^ 


pdV 


(5.20) 


Thus we clearly observe that just like for DGP model as discussed in the earlier sec¬ 
tion, here also the increase in amplitude of the scale factor depends on the parameter 
of the model i.e the variants of cosmological constant as mentioned in case 1, case 2 
and case 3. 

( 

pdV for A: = -1-1 


3M1 (1 - A<.„) 


^^min A 


(5.21) 


pdV for k = —1. 


t 3M2 (1 + AalJ 

Let us now briefly mention the characteristic feature of the results for cosmological 
bounce for the previously mentioned variants of cosmological constant model in the 
following: 


1. For a closed universe, with A; = -|-1: 
• The hysteresis loop integral 


pdV < 0 


(5.22) 


56 



or equivalently 


(5.23) 


Pexp ^ Pcont 

if (1 — Aa^j„) > 0 or < 1 and 5amin > 0, 


The hysteresis loop integral 


or equivalently 


pdV > 0 


Pexp ^ Pcont 


if (1 — Aa^j„) < 0 or Aa^j„ > 1 and Samin > 0. 

2. On the other hand, for the open universe, with k = —1: 


(6,24) 


(5.25) 


• The hysteresis loop integral 


pdV < 0 


or equivalently 

Pexp ^ Pcont 

if (1 + < 0 or Aa^^^ < -1 and > 0, 

• The hysteresis loop integral 


pdV > 0 

or equivalently 

Pexp ^ Pcont 

if (1 + > 0 or Aa^j„ > -1 and Samin > 0. 


(5.26) 

(5.27) 


(5.28) 

(5.29) 


For the other two cases i.e case 2 and case 3, since the Friedmann equations re¬ 
main the same, the condition for bounce in both the cases will be given by Eq. (5.20). 
Hence the above analysis holds true for case 2 and case 3 also. 


5.2 Condition for acceleration 
For case 1 

From Eq. (5.2), at bounce the condition for acceleration is given by 

Pb + 3pb < 2AM2 (5.30) 
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Substituting the expression for pi, for different values of k we get the conditions 
for acceleration at bounce as: 


AM^-^ for k = +1 

at 


Pb= < 


AM^ + 



for k = —1. 


(5.31) 


Thus we see that the condition for acceleration at bounce now depends not only on 
the scale factor but also on the cosmological constant. 

Substituting the expressions of p and p from Eqns. (2.5) into Eq. (5.30), we get 
the condition for acceleration in terms of the scalar held as 


02 < E(0) + AM2 (5.32) 

Thus we see that even if the contribution from the potential is lesser than the standard 
canonical kinetic term i.e. even if 0 ^ > \/(^00 the above condition will still be satished 
if 0^ is lesser than the sum of E(0) + AM^. Thus we can use potentials of smaller 
values than standard case and yet produce bounce. Thus the presence of cosmological 
constant modihes the effective potential present at bounce. 

In Figs. 8 (a) and 8 (b), we have shown the phenomena of bounce and acceler¬ 
ation in this cosmological constant dominated model. We can draw the following 
conclusions from the above hgures: 

• In Fig. 8 (a), we have plotted the r.h.s of Eqn. (5.1) using the relation p = 
a-3(i+"'). The graph has been plotted only for the case when k = 1. This is 
because, Eq. (5.15) shows that density becomes negative for the cases when 
A; = 0 , — 1 , which is unphysical. 

• We have used tc = 1/3, since we require a soft equation of state for causing the 
acceleration and expansion. The case tc = 0 made H only approximately zero. 

• From Fig. 8 (a), using the value of A from Planck 2015 data [75], we get the 
bounce at p = pb = OASM^i. This value increases if we take higher values of A. 

• Fig. 8 (b) shows the necessary condition of acceleration (a > 0) at the time of 
bounce. Here we have plotted the r.h.s of Eqns. (5.1) and (5.2). This plot has 
also been obtained for the same parameter values as the earlier graph. 

• Thus Fig. 8 shows graphically that the phenomenon of bounce is possible for 
a closed universe dominated by cosmological constant like term and having 
w = 1/3. 
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H vs p for cosmological 
constant dominated Einstein gravity 



a VS p for cosmological 


constant dominated Einstein gravity 



(a) An illustration of the bouncing condition for (b) An illustration of the acceleration condition 
a universe with an equation of state w=l/3 and at the time of bounce for a universe with an equa- 
k=l. tion of state w=l/3, k=l. 


Figure 8 . Graphical representation of the phenomena of bounce and acceleration for 
cosmological constant dominated Einstein gravity. 


• Since the condition of bounce is true for all the three cases being studied in this 
section, we can say that bouncing universe is possible for all the three cases. 

For case 2 


Condition for acceleration in terms of the pressure and density is again given by 
Eq. (5.30). Then substituting the expressions for p and p from Eq. (5.9), we get the 
condition for acceleration in terms of the scalar held as: 

4 

02 < 2AoM2 + 1/(0) - (5.33) 

i=l 

Thus depending on the signs of constants (Ai, A 2 , A 3 , A 4 ), we observe that less or 
more contribution from the potential than the canonical kinetic term as appearing 
in the standard case may be required. To realize the essence of this statement let us 
consider an example, in which all the new constants are positive or the resultant sign 
of the summation is positive, then this decreases the right hand side from standard 
case given by Eq. (5.32) and this condition for acceleration can be achieved only 
if the contribution from the potential is more than compared to the standard case. 
Again if all the constants are negative or the resultant sign of the summation terms is 
negative, then we need lesser contribution from the potential than the standard case 
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to attain the condition of acceleration. Thus we see that apart from the general cos¬ 
mological constant term, the presence of different powers of scalar fields along with 
different constants increases or decreases the effective potential at the time of bounce. 

For case 3 

Condition for acceleration in terms of the pressure and density is again given by 
Eq. (5.30). Then substituting the expressions for p and p from Eqns. (5.14), we get 
the condiition for acceleration in terms of the scalar held as: 


y < 2A(,iVn + V{<p) - Ae*/'^r 


(5.34) 


Thus if A is positive, then this decreases the right hand side from standard case 
and this condition for acceleration can be achieved only if more potential than the 
standard case is present. Again if A is negative, even potential lesser than the 
standard case can satisfy the condition of acceleration. 

Therefore from the above analysis we can conclude that, whether the standard 
potential will be able to cause acceleration at bounce now depends on the extra terms 
present in the action. 

5.3 Condition for turnaround 
For case 1, 2 ,3 

The condition for turnaround is exactly same as the condition for bounce appearing 
in case 1 i.e. one can write: 



(5.35) 


Rest all the conclusions remain same as what we had got for the bounce case. It is 
important to mention here that the conclusions for case 2 and case 3 also remains 
the same as appearing in the case 1. 

5.4 Condition for deceleration 
For case 1 

From Eq. (5.2), at turnaround the condition for deceleration is given by 


Pt -|- 3pf > 2AAf^ 


(5.36) 


i.e turnaround can be obtained without violating the strong energy condition. 






In place of Eq. (5.31), we get the conditions for deceleration at tnrnaround as: 




M2 

—^ for k = +1 


Pt> < 


(5.37) 


AM^ -|-— for k = —1. 

at 

Substituting the expressions of p and p from Eqns. (2.5) into Eq. (5.36), we get the 
condition for expansion in terms of the scalar held as 

0^ > E(0) + AM2 (5.38) 

Thus we see that a potential lesser than the standard case is required to satisfy this 
condition because even if the condition (j) > V ( 0 ) is satished, we will get acceleration 
only if 0 > V{(j)) + AM^. Thus, now the cosmological constant increases the effective 
potential relative to the standard case. 


H vs p for cosmological 
constant dominated Einstein gravity 



(a) An illustration of the turnaround condition 
for a universe with an equation of state w=l and 
k=l. 


a vs p for cosmological 
constant dominated Einstein gravity 



(b) An illustration of the deceleration condition 
at turnaround for a universe with an equation of 
state w=l, k=l. 


Figure 9. Graphical representation of the phenomena of turnaround and deceleration for 
cosmological constant dominated Einstein gravity. 


In Figs. 9(a) and 9(b), we have shown the phenomena of turnaround and decel¬ 
eration in this cosmological constant dominated model. We can draw the following 
conclusions from the above hgures: 
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• In Fig. 9(a), we have once again plotted the r.h.s of Eq. (5.1) using the relation 

p = for a closed universe. But for this case we have used w = 1, since 

we require a stiff equation of state for causing contraction and deceleration. 

• In Fig. 9(a), H going to negative values may be interpreted as the universe 
changing its direction of motion at turnaround. As had been discussed in [10], 
the condition of bounce/turnaround can be imposed by either the condition of 
making the scale factor changing sign with other quantities remaining same, 
or, a going to negative values with other quantities remaining same. 

• From Fig. 9(a), using the value of A from Planck 2015 data [75], we get the 
turnaround sX p = pt = O.OldM^;. This value increases if we take higher values 
of A. 

• Fig. 8(b) shows the necessary condition of acceleration (a < 0) at the time of 
turnaround. Here we have plotted the r.h.s of Eqns. (5.1)and (5.2). This plot 
has also been obtained for the same parameter values as the earlier graph. 

• Thus Fig. 10 shows graphically that the phenomenon of turnaround is possible 
for a closed universe dominated by cosmological constant like term and having 
w = 1. 

• Since the condition of turnaround is true for all the three cases being studied 
in this section, we can say that bouncing universe is possible for all the three 
cases. 

For case 2 


Condition for deceleration in terms of the pressure and density is given by Eq. (5.36) 
with A replaced by Aq in this specific case. Then substituting the expressions for den¬ 
sity p and pressure p from Eq. (5.9), we get the following condition for deceleration 
in terms of the scalar field as: 

4 

02 > 2AoM2 + 1/(0) - (5-39) 

i=l 

Thus depending on the signs of constants (Ai, A 2 , A 3 , A 4 ), we see that less or more 
contribution from the potential than the standard case may be required. To justify 
the validity of this statement let us consider an example, where all the new constants 
are independently positive or the cumulative effect appearing through their summa¬ 
tion is positive, then this decreases the right hand side of Eq. (5.39) from standard 
case and the condition for deceleration can be achieved even if more contribution 
from the potential plays crucial role compared to the standard case. On the other 
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hand, if all the constants are independently negative or the cumulative effect ap¬ 
pearing through their summation is negative, deceleration can be achieved even if 
the contribution from the potential is lesser than the standard case. 

For case 3 


Condition for deceleration in terms of the pressure and density is again given by 
Eq. (5.36) with A replaced by Aq. Then substituting the expressions for density p 
and pressure p from Eq. (5.14), we get the condition for deceleration in terms of the 
scalar field as: 

02 > 2AoM2 + 1/(0) - (5.40) 

Thus if the signature of A is positive, then this decreases the right hand side from 
the standard case and this condition for acceleration can be achieved even if more 
contribution from the potential compared to the standard case is present. On the 
other hand, if A is negative, we need lesser potential than the standard case to attain 
the condition of deceleration. 

For k = 0 case 


The results below have been shown only for case 1. But they perfectly hold good for 
case 2 and case 3 also. In the present context the Friedmann equations for this case 
are given by: 


H + H‘^ 



a {p + 3p) 


6M2 




(5.41) 

(5.42) 


From Eq. (5.42), we can infer that the same Friedmann equation cannot give rise to 
both bounce and turnaround condition for k = 0 case. This is because A is constant 
in the present setup, hence it is not possible for it to be equal to both the high 
density at the time of bounce and low density at the time of turnaround. 

For causing bounce/turnaround the following condition holds good: 

pt/t = -AM2 (5.43) 

Therefore in place of Eq. (5.20), in the present context we get: 

^iflmin/max) AAf^ f 


Thus, the amplitude of the scale factor increases: 
• if ^ pdV > 0 for A > 0. 


(5.44) 







if j> pdV < 0 for A < 0. 


Similarly the condition for acceleration/deceleration for fc = 0 case can be written 
as: 


Pb < AM^ (for acceleration at bounce) (5.45) 

Pb > AM^ (for deceleration at turnaround) (5.46) 

Thus we observe that the pressure at bounce/turnaround required to cause accelera¬ 
tion/deceleration becomes independent of the scale factor at the time of bounce/turnaround. 
Also the results shows that it depends only on the cosmological constant. 

5.5 Evaluation of work done in one cycle 
For case 1 


Following the same procedure as we have done for DGP model, here also the general 
expression for the hysteresis loop is given by Eq. (4.52). In order to express the loop 
in terms of the scale factor only, using Eq. (5.1) and Eq. (5.2), we get: 


P = y + ^(0) = 3M2 


k A 
“I H—^ ~ w 


+ H = - = 
a 


jp + 3p) A 
6M2 3 


Further using the above equations and Eq. (2.5) for pressure p we get: 


2 


1 /( 0 ) = 




(5.47) 

(5.48) 


(5.49) 


Therefore the expression for work done is given by: 



3M" 


(Aa^a 


2 dda 


kd) dt. 


(5.50) 


It is important to note that the total work done in one cycle now depends on the 
cosmological constant also. 


For case 2 


For case 2 using the expression of p given by Eq. (5.9) into Eq. (4.52), we get the 
expression for the total work done in one expansion-contraction cycle in terms of the 
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0 and a as 




(5.51) 


(5.52) 


Since the Friedmann equations are same as that appearing in the earlier section, 
the expression for the above equation in terms of the scale factor is again given 
by Eq. (5.50) in which A is replaced by a new constant Aq. Thus from Eq. (5.52) 
we find that the work done now depends on the value of all the new constants 
(Ai, A 2 , A 3 , A 4 ) of the model. Thus by allowing sufficient amount of tuning in 
these constants we can make the signature of the representative integral positive or 
negative. 


For case 3 


For case 3 using the expression for pressure p as given by Eq. (5.14), into Eq. (4.52), 
we get the expression for the total work done in one expansion-contraction cycle in 
terms of the scalar field (j) and scale factor a as: 



' cont 


3 ( -0^ — V{(j))+ o?adt 


’ exp 


3 ( -0^ — V{(j))+ a^ddt 


(5.53) 


3 ( - 


-*'max 


02 - 1 /( 0 ) + a^a 


ddt 




f^max / 1 


a^ddt 


(5.54) 


In this case also, the expression for the above equation in terms of the scale factor 
is again given by Eq. (5.50) where A is replaced by the constant Aq. Thus adjusting 
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the values and signatures of the constants (Aq, A) we can change the signature of the 
representative integral. 

5.6 Semi-analytical analysis for cosmological potentials 

Here also we will perform the analysis for the case when /c = 0. Though, the exact and 
complete treatment would have been for the case k = 1, simple analytical solutions 
were obtained only for flat universe which has been discussed below. Also, since we 
know that for /c = 0, this model can cause either bounce or turnaround, hence to 
get an expression for the complete work done in one cycle,, we will consider that 
cosmological constant is present at late times (i.e causing turnaround) whereas the 
early universe is governed by brane world cosmology model for which the analysis 
has been shown in the Appendix in the context of Randall-Sundrum single brane 
world (RSII). Now, if we compare Eq. (9.2) with Eq. (6.30), we see that they are 
both same. Hence the results for early universe which will be calculated for loop 
quantum model, can be used here. Below, we show the late time analysis for ACDM. 
In this context, we have denoted the Planck mass by Mp. 

5.6.1 Case I: Hilltop potential 
For Case 1: 


A. Expansion 


At late times within the window a' < a < Qmax using Eq. (5.1) and Eq. (4.67) with 
the expression for density p now given by the hilltop potential, we get an integral 
equation of the following form; 


/ 


d 






(5.55) 


The exact solution of this integral is given in the Appendix for RSII model. To 
compute the left hand side of the above integral equation we again follow the same 
procedure as we done for DGP model i.e. we redehne the held as: 



and hence using this new dehnition we study two limiting cases. 


(5.56) 


a)4>/Mp « 1: 


66 










In this limit, we can expand the exponentials npto linear order and npon further 
simplification we get the following expression for the redefined field A as: 


A 




I3p 


(3^(1 + + ty 

where Bq is the arbitrary integration 

r“- 


(5.57) 


constant given by: 


5, 


>0 




3M2 

Here A/ is the value at the time of turnaround i. 

The above solution has been obtained under 


(5.58) 


the time of turnaround i.e aX t = tf. 

has been obtained under the assumption that the quantity 


Vq 

7- — -v/SpA « 1, 

(^(1 + /^) + t) 

which is possible if we choose the value of Vq and A accordingly. 

Substituting this expression back into the Friedmann equal 
pression for the scale factor as 


(5.59) 


a{t) = Bi exp 



3^(1 + / 3 ) + xI + 




equation, we get the 
_Yo_[j 2 2 


ex- 


7-'-^ 


where 

Hi 


(3]%) t ? 

4(jj(l+/3) + f) 2 


(5.60) 


: af exp 



^( 1 +«+!)+ 


_hL/52 2 \ 

3My P 

—7-:v-°o tf 

2 (3^(1+ /^) + 3) / 

/ \2 _ „ 1 


4 f_hL 


{fkypv 


.n j_ _L a' 


i 


Here Of is the 


ii) <p/Mp » 1: 
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In this limit, the solution for the redefined held A is given by: 


A = 


2 

P 



where B 2 is the arbitrary integration constant given by: 


Bo = 




Vat 


ot/ 


pi/2p2 


3M2 


(5.62) 


(5.63) 


The above solution has been obtained under the assumption that the conditions: 


» 1 , 

Vol3eP^ A 

3M2 3 ’ 


(5.64) 

(5.65) 


are satished. These conditions can be satished since we are in the large held limit 
and we can choose the values of the parameters of this model {/3,p, Vq) accordingly. 

Substituting the above expression into the Friedmann equation, we get the ex¬ 
pression for the scale factor as: 



where B^ is the arbitrary integration constant given by: 


S3 = ttfCXp 


In 


V3M2; 




1/2 


+ 


/ \ 1/2 


/3l/2p2 


,2 / _Y0 


1/2 N 


( 5 . 67 ) 



B. Contraction 


This phase is independent of any choice of the potential provided the condition 
0^ >> ld( 0 ) holds good in this phase. Using Eq. (4.68) and Eq. (5.1), we get the 
following solutions: 


A54 f + A 
_ 1 I 3 


where is the arbitrary integration constant given by: 


^ A\ 3 /- i' 


3j A ye^A + lJ ■ 

From the above expressions, we get the solution of the scale factor as: 


a{t) = 






oV^t 


1 / AB 4 


where B^ is the arbitrary integration constant given by: 


B5 — af 




-If 




1 / A-B4 


(5.68) 


(5.69) 


(5.70) 


(5.71) 


C. Expression for work done 

As will be discussed in the next section, that the expression for integral for work 
done in the early universe for loop quantum model contains a large no. of terms, 
hence will not be shown here. But the results of the integral corresponding to late 
universe have been discussed below. 

i) 0/Mp << 1: 


pdV = bi (erf [-62 + b^tmax] - erf [-62 + b^t']) + 3Mp{4:5smh{VAtmax) 

- 9 smh{2y/Atmax) + smh{3y/Atmax) - 45 sinh(A/Af^) 

+ 9 sinh(2yAf') - sinh(3yAT) - 30tmax + 30t') (5.72) 
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ii) cj)/Mp » 1: 


/ 


pdV 


+ A)(Erfi(A/36^(66 + bytmax)) - Erfi(A/36^(66 + Mmm)) 

-Oq-;- 

^ 267 

+ 3Mp{4:5smh{^/At„^,ax) 

— 9 sinh( 2 V^tmax) + sinh(3A/Atmaa;) — 45 sinh(A/At') 

+ 9 sinh(2-\/At') — sinh(3V^t') — SOtmax + 30t') (5.73) 


Here bi-.-bj are constants that depend on the parameters present in the expression 
for the scale factor. Their explicit forms have been given in the appendix. 

Including the contributions from the other integrals also we see that we get a 
non-zero work done for both large and small field cases for hilltop potential when we 
have pure cosmological constant in the background. 


For Case 2 


A. Expansion 


Within the window Omm < a < ttmax at late times using Eq. (5.5) and Eq. (5.9) 
along with energy density of scalar field, p, now given by the hilltop potential, we 
get an integral equation of the form (neglecting 0): 

/r^o(i+/3(4,)")-EtiAhMW(4;h a; 

/ 3M2 3 

For the sake of clarity, we again follow the same procedure as we have done for DGP 
model i.e. redefine the field variable, as: 


3M„ 


dt. 


(5.74) 



A 

Mp 


e . 


(5.75) 


Simplified analytical expressions can only be obtained for the small field limit given 
by 0/Mp << 1, so that we will concentrate in this specific physical situation for 
this case. In this limit, we can expand the exponentials upto linear order and upon 
further simplification we get the following expression for A as: 


(a+ 

^/K 


(5.76) 
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where the arbitrary integration constants J, K and -B 14 is given by: 


Ai A 2 AgMp A 4 M, 


3M„ ^ 3 ^ 


j ^0 Aq 

3M2 ^ 3 ^ 3M2 ■ 

^ + 2A2Mp + 3A3M2 + 4A4M3, 


+ 


5i 4 = A,(v^/i^) - 


3M„ 


(5.77) 

(5.78) 

(5.79) 


Here A/ is the value at the time of turnaround i.e. at the time scale t = tf. 

The above solution has been obtained under the assumption that the conditions; 


JA 

M << 
N\ 


(5.80) 

(5.81) 


are satished. This is possible because we are in the small A limit and we can choose 
the values of the other constants accordingly. Here we introduce two new constants 
M and K dehned as: 


M = 


VqPP 

3M2 

Hn 


Ai 2A2 3A3MP 4A4M2 


3Mr 




K = -^/3p{p - 1) + 2A2Mp + 9AsM^ + 


(5.82) 


(5.83) 


Further substituting this expression back into the Friedmann equations, we get 
the following expression for the scale factor as: 


a(t) = i?i 5 exp 




6MpJ3/2 J 


2 \QMp.PA 


(5.84) 


where H 15 is the arbitrary integration constant given by: 


Hi 5 = a/ exp 


/ BuKM\ t) f KM \ 

^ 6MpJ3/2 j + Y \mp.PA) 


(5.85) 


Here a/ is the value of scale factor at turnaround time scale t = tf. 


B. Contraction 


The conclusion remains same as for ACDM model, because in this case all the extra 
terms containing 0 , will be neglected at leading order approximation. 

C. Expression for work done 
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The solution of the scale factor (i.e its time dependence) is same as that for hill¬ 
top potential for constant A case. Hence the expression for work done is same except 
that the constant will now depend on the parameters of this model. Thus the phe¬ 
nomenon of hysteresis also holds true for this model. 

For Case 3 


A. Expansion 


At late times within the window a' < a < ttmax, using Eq. (5.12) and Eq. (5.14) 
with scalr field density p now given by the hilltop potential, we get an integral equa¬ 
tion of the following form (neglecting 0): 


/ 


d 





3M2 


+ 







(5.86) 


In order to simplify the above integral for computational purpose, we follow the same 
procedure as we have already done for DGP model i.e. redefine the held by using: 


A 

Mp 


= e 


(5.87) 


Simplihed solutions were possible only for small held case i.e. for 0/Mp << 1, which 
has been discussed below. 

In this limit, we can expand the exponentials upto linear order and upon further 
simplihcation we get the expression for A as 


A = 


3M, 


+ B 


18 




+ 


Vq I Vo^^ 
3M2 3M2 


Ae I Ao 
3M2 ■T- 3 


where Big is the arbitrary integration constant given by: 


.Bis — 


X , / Vq I ^OP _ Ae I Ao 

^T\/ “T n^,f2 QA/f 2 “ Q 




Ae 


3M2 


3M2 


3Ar2 


Ao 


' —VqPp I Ae 5 

' Mp ' Mp> 


3M„ 


(5.88) 


(5.89) 


Here Aj is the value at the time of turnaround i.e at the time scale t = tf. 

The above solutions have been obtained under the assumption that the condi¬ 
tions: 


{Vo/3p-Ae)\ 

(Co + I/o0-Ae + AoM2) 
((Ae - Vo(3p{p - 1)) /Mp)X 
{i-VoPp + Ae)/Mp) 


«1, 

(5.90) 

« 1 

(5.91) 
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are satisfied, which is possible if we choose the values of the constants accordingly. 

Further substituting this expression back into the Friedmann equation, we get 
the following expression for the scale factor as: 


a(t) = Biq exp 


Va((i + 


V V" ' 20'3/2 J ^ J 

where Big is the arbitrary integration constant given by; 


Big = ttfCXp 




P’tj 


VV 20'3/2y ^ 


in which we introduce two new symbols given by: 

^0 VoP Ae 

3M2 ^ 3M2 3M2 ^ 

, ^ l' Vo/3p-Ae \ / Vol3p 
~ ^ 3M2 )\ Mp 

Here a/ is the value of scale factor at turnaround. 


3 ’ 


(5.92) 


(5.93) 


(5.94) 

(5.95) 


B. Contraction 

Solutions remain same as obtained before. So the conclusion also remain unchanged 
for contraction phase of the Universe. 

C. Expression for work done 

The solution of the scale factor (i.e its time dependence) is same as that obtained for 
hilltop potential. Hence the expression for work done is exactly same except that the 
constant will now depend on the parameters of this model. Thus the phenomenon 
of hysteresis also holds good for this model. 


5.6.2 Case II: Natural potential 
For Case 1: 


A. Expansion 
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At late times within the window a' < a < amaxi using Eq. (4.67) and Eq. (5.1), 
we get the following integral equation of the form: 


/ 


d 




(5.96) 


The exact solution of the left hand side of the above equation is given in the Ap¬ 
pendix for RSII model. For the sake of simplicity, here we study solutions for two 
limiting physical situations. 


i) 0// « 1: 


Using the small argument approximations for trigonometric functions, we get the 
expression for scalar held (j) as: 


/ 


Bq exp 


3/2 



where Bq is the arbitrary integration constant given by: 


(5.97) 


Be = —exp 


Vetp 


3/2 1 + 


Vo 

3mJ 


1/2 


Vo , A 
3M| 3 


(_L A 
I 3M2 


1/2 


(5.98) 


Here (f)p is the value of the scalar held at turnaround corresponding to t = tp. 

Further substituting the expression for 0 into the Friedmann equation, we get 
the following expression for the scale factor as: 

1/2 


a(t) = By exp 


1 + 


Vo 

3M2 


1/2 


(^ + -] t 
3% + ty V3M2 3 I 


(5.99) 


where By is the arbitrary integration constant given by: 

X 1/2 


By = Oi? exp 


- 1 + 


3M2 


/- 


Un A 


Vo A / V 3/1+2 

3M2 + 3 / V'J 


1/2 


+ 3 I iF 


(5.100) 


Here ap is the value of the scale factor at turnaround. 
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ii) 0// » 1: 


Since we know that for large argument (p/f » 1, value of cosine function is very 
small, hence using this concept and the following constraint condition; 


cosec — >> 


3M2 cos ^ / 


fJ + f) (I 


we get the solution for the scalar held p as: 


Ht) 

f 


= 2 tan 


-1 


( Vpt 

V3M2 


+ B, 


Vo I A 
3M2 T 3 


where Bs is the arbitrary integration constant given by: 


(5.101) 


(5.102) 


Ss = tan ( ^ 


Vo A ^ Votp 
^ 3" ~ sM 


(5.103) 


Here pp is the value of the scalar held at turnaround. 

Further using the expression for p in the Friedmann equation, we get the follow¬ 
ing expression for the scale factor as: 


a(t) = Bgexp 


Hn A 


3M2 + 3 ) 


1/2 


where G{t) is dehned as: 

/ 


G{t) = 


Vo 

3mI 


t 


R„ AqF 1 
2 ( ^0.. + A 1 V ) 


Vo 

^ T72 tan-1 
Vo , A ■ 


f-hn) 

\ZMl) 


V 


and i?8 is the arbitrary integration constant given by; 


Bg = apexp 


Hn A 


1/2 


(5.104) 


Bs+ 


Vot 


\ 


Vo , A 


1/2 


(5.105) 


3M2 + 3 ) 


(5.106) 


Here ap is the value of the scale factor at turnaround. 
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B. Contraction 


This phase is independent of any choice of the potential provided the condition 
0^ >> V{(j)) holds good in this phase. This implies that the hnal conclusion remains 
same as obtained for previous hilltop potential in the background of ACDM model. 

C. Expression for work done 

As will be discussed in the next section, that the expression for representative in¬ 
tegral for work done in the early universe for loop quantum model contains a large 
number of terms, hence will not be explicitly shown in the present context. But the 
results of the integral corresponding to the late universe have been discussed below. 

i) 0// « 1: 


In this limit, the total work done is given by: 

y pdV = - 1 - 3Mp(45sinh(VAtmax) - 9 sinh(2VAf^aa:) 

-t- sinh(3V^tmax) — 45 sinh(V^t') -|- 9sinh(2V^0) — sinh(3V^0) 

-30t^,, + 30T).(5.107) 

Here 64...69 are constants that depend on the parameters present in the expression 
for the scale factor whose explicit expressions have been given in the appendix. 

Including the contributions from the other integrals also we see that we get a 
non-zero work done for small held case for hilltop potential when we have cosmolog¬ 
ical constant in the background. 

ii) 0// » 1: 


In this limit, the total work done is given by: 



X 


3 &io (—1 + &11&12 — &12&13) 

[-36^4 max —bii+bi3){bis+bi2t 

max ) _ g36io(t min fell+fel3)(6l5+fel2t™n)^ 

-j- ^4 ^g36io(i max —bii+bi3){bi5+bi2t 

max ) _ g3feio(t min fell+fel3)(fel5+6l2t™„)) 
max — bll+bl3){bis+bl2tmax) g^lO(imin—^l 1+^>13)(^15+^12*771171) ^ 

-h 3Mp{4:5smh{y/Atmax) - 9 sinh(2yAfmax) + sinh(3yAfmax) 

— 45 sinh(v/At') -|- 9 sinh(2v/At') — sinh(3V^tO — 30tmax + 300 

(5.108) 
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Here 64...614 are constants that depend on the parameters present in the expres¬ 
sion for the scale factor. Thus we get a non-zero work done for large held case for 
hilltop potential when we have cosmological constant. 

For Case 2 


A. Expansion 


Analytical solutions in this case can also be obtained for small held limit i.e. when 
0// << 1. Hence taking the small angle approximations of the trigonometric func¬ 
tions and redehning the held using following transformation equation; 


/ ’ 

we get the following solution for the transformed held A as: 

where Bi^ is the arbitrary integration constant given by; 

_ yiiT tp 

16 F 2^ 

where we introduce two new constants K' and M' given by: 

A' = y + Ai + 2A2Mp + 3 A 3 M 2 + 4A4Mp^ 


. 2Vq Aq 

M = —- H- - 

3M2 ^ 3 


Ai ^ A2 ^ AgMp ^ 


3Mp 3 


(5.109) 

(5.110) 

(5.111) 

(5.112) 

(5.113) 


Here AiT’ is the value at turnaround. 

The above expressions have been obtained under the assumption that the fol¬ 
lowing constraint conditions are satished: 


L'\ 

AT 

N'X 


«1, 

« 1 . 


(5.114) 

(5.115) 


Since we are in the small held limit, these conditions can be satished if we choose 
the numerical values of the model parameters accordingly. Here we additionally 
introduce two new constants L' and N' given by: 


Ai 2 A 2 3A.^A4p AA.4^A4p 

3A^ ^ IT ^ 3 ^ 3 


(5.116) 
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N' = j + 2A2Mp + 9A3M2 + I2A4M, 


( 5 . 117 ) 

The expression for the scale factor that follows from the above solution is given by: 


a{t) = Bn exp 


n _ L'K'e \ 

VV 2M'3/2 ) VlM'^Bf) 


where Byj is the arbitrary integration constant given by: 
Byj = Of exp 


L'K't\ 


2M'3/2 


Here is the value of the scale factor at turnaround. 


(5.118) 


(5.119) 


B. Contraction 


The conclusion remains same as for ACDM model, because in this case all the extra 
terms containing 0, will be neglected at leading order approximation. 

C. Expression for work done 

The solution of the scale factor (i.e its time dependence) is same as that for hill¬ 
top potential. Hence the expression for work done is same except that the constant 
will now depend on the parameters of this model. Thus the phenomenon of hysteresis 
also holds good for this model. 


For Case 3 


At late times within the interval a' < a < amax, using Eq. (5.12) and Eq. (5.14), 
we get the integral equation of the form similar to Eq. (5.86) with H(0) given by the 
natural potential. Simple analytical expressions for the scalar held and scale factor 
could be obtained only for the case when 0// << 1 which has been discussed below. 
In this limiting situation using the small argument approximations for trigonometric 
functions, we get the following expression for the scalar held 0 as: 


Ht) 

f 



A 


f j 2Vh A I Aq ’ 
J Y 3Ar2 3M2 3 


(5.120) 


where B 20 is the arbitrary integration constant given by: 


B 20 — 


/ 2^0 _ A I ^ 

3M2 3 


A 


tp 


(5.121) 
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Here (j)F is the value of the scalar held at turnaround corresponding to t = tp. 

Finally, substituting the expression for 0 into the Friedmann equation, we get 
the following expression for the scale factor as: 


a{t) = B 21 exp 


fp, ( q 

tv 2Q«l^fM^ 


t — 




where B 21 is the arbitrary integration constant given by: 


B 21 = Of exp 


Ill- 

VV 


tp — 


AHl 




Here we introduce a new constant Q' dehned as: 

, ^ 2Vo _ A ^ 

^ 3M2 3M2 3 

Here ap is the value of the scale factor at turnaround. 


(5.122) 


(5.123) 


(5.124) 


B. Contraction 

Solutions remain same as obtained before. So the conclusion also remain unchanged 
for contraction phase of the Universe. 

C. Expression for work done 

The solution of the scale factor (i.e its time dependence) is same as that for hill¬ 
top potential. Hence the expression for work done is same except that the constant 
will now depend on the parameters of this model. Thus the phenomenon of hysteresis 
also holds good for this model. 


5.6.3 Case III: Coleman-Weinberg potential 
For Case 1: 


A. Expansion 


At late times within the interval a' < a < amax, using Eq. (4.67) and Eq. (5.1), 
we get the following integral equation of the form: 


d ^ 



Un 


3 r 


^ 4a + /3-F4/31nU^ 


3M2 


dt. (5.125) 
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To compute the left hand side of the above integral equation, we once again use the 
field redefinition via the following transformation; 

(5. 

and study the two limiting physical situations as mentioned below, 
i) 0/Mp << 1: 




In this limit, we expand the exponentials upto linear order and get the resulting 
solution for the transformed held A as: 


A = 


Vh 


3M2 


t + B 


10 


(4q; + /3) 




1/2 ’ 


(5.127) 


where Biq is the arbitrary integration constant given by; 


-Bio — 


( j/l _l Yob. _l \ 

1^3Af2 3^2 3 y Vq 


4q; T /3 




(5.128) 


The above expressions have been obtained using the assumptions that the constraint 
conditions: 


4/3X 

(4q; + /3) 

^{4a + /3)X 

Vq I Vpa I A 
3M2 3M^ ^ 3 


(5.129) 

(5.130) 


are satished. Since we are already in the limit in which the value of A is very small, 
hence by choosing the other parameters of the model properly, we can easily satisfy 
the above constraint conditions. 

Further substituting the above expression into the Friedmann equation, we get 
the following expression for the scale factor as: 


a(t) = Bii exp 


y ^ p p ^ ^ 3M2 + 3 y ) 

iz. . . n 


I/o^(4a + A) 

-^-r-(4a + A)A 

12 I F) I Yob. a ) A/f 2 

\3M2 ^ 3M2 ^ 3 J 
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where -Bn is the arbitrary integration constant given by; 
Bii = a/exp 


' - 71 4-0 ~r I ~r 


V3M2 3 M2 3 


^(4a +13) 


9 f ^0 I Vb« I A^i 
^ l^3M2 T 3^2 3 y 


-Bio(4q; + /3) tf 


Vo^{4a + /3) 


1 2 I ^0 -L Yq^ 4- Y] M2 
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.(5.132) 


Here a/ is the value of the scale factor at turnaround. 


ii) 0/Mp >> 1: 


In this case we get an integral of the form 


3^ + t + 3^(“ + 


4A 


Vo 


/ d\ ^ / dt. (5.133) 

J (4a + /3 + 4/3\) 3>Mp j 

To compute the left hand side of the above integral equation and hence to get an 
analytical expression for A, we assume the following three constraint conditions i.e. 

K) A\ _ Ho 


A 

3M2 ^ i 


<< 


3Af| 


(a T /3\)€ 


4A 


4q; T 13 4/3A, 

a << /3A. 


(5.134) 

(5.135) 
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Since we are in large A limit, by choosing the values of the parameters of the model, 
the above conditions can be easily satished. Under these assumptions, we get the 
expression for A as: 


A = 


Un 


3M2 


f+ -B 


12 


4/9 

3M2 


where B 12 is the arbitrary integration constant given by; 
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(5.137) 


(5.138) 


Here A/ is the value at turnaround. 

Further substituting the expression for A, hence 0 into the Friedmann equa¬ 
tion and applying the above mentioned constraint conditions, we get the following 
expression for the scale factor as; 


a(t) = -B 13 exp 
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(5.139) 
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where i?i 3 is the arbitrary integration constant given by: 


Biz = a/ exp —2e V 




(5.140) 


Here af is the value of the scale factor at turnaround. 


B. Contraction 


This phase is independent of any choice of the potential provided the condition 
(jp' » V{(j)) holds good in this phase. This implies that the hnal conclusion remains 


same as obtained for previous hilltop and natural potential in the background of 
ACDM model. 

C. Expression for work done 

The solution of the scale factor (i.e its time dependence) is same as that for hill¬ 
top and natural potential. Hence the expression for work done is same except that 
the constant will now depend on the parameters of this model. Thus the phenomenon 
of hysteresis also holds perfectly for this model. 

For Case 2 

A. Expansion 

At late times within the window a' < a < amax, using Eq. (5.5) and Eq. (5.6), we 
get an integral equation similar to Eq. (5.74), with the potential V{(j)) now given by 
the Supregravity motivated Coleman-Weinberg potential. For the case 2 the math¬ 
ematical form of the governing integral equation is exacly same as that obtained in 
case 2, only in the present case the constant A is replaced by the held dependent 
A(0), which is dehned earlier. To compute the left hand side of the master integral 
equation we will once again use the previously mentioned held redehnition: 



One can hnd the analytical solutions only for the case when the small held limiting 
approximation is valid i.e. p/Mp « 1 is satished. 

In this limit, we expand the exponentials upto linear order and get the resulting 
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solution for the redefined field A as: 


A 




where Bi^ is the arbitrary integration constant given by: 


(5.142) 


Bie — 
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tf 

3Mp' 


(5.143) 


Here we introduce two new constants O and P given by the following expressions: 


O = -(4a + /3) + Ai + 2A2Mp + + 4A4Mp^ (5.144) 


_ Vo I ^0 _ / Ai ^2 A^Mp A^M‘^\ 

~ 3m| ^ y “ ^ y ^ 3 ^ 3 ) ■ 

Here A/ is the value at turnaround when t = tf. 

The above expressions have been obtained using the assumptions that the con¬ 
ditions, 
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(5.146) 

(5.147) 


are satisfied. Since we are already in the limit in which the value of A is very small, 
hence by choosing the other parameters of the model properly, we can easily satisfy 
the above conditions. Here we introduce two new constants Q and L given by the 
following expressions: 


M ^ + SAgMp 4A4M^y 

3M2 ^ 3M2 \3Mp ^ 3 ^ 3 ^ 3 / ’ 


(5.148) 


L = - {7/3 + 12a) + 2A2Mp + GAgM^ + UA^M^. (5.149) 

Further substituting the above expression into the Friedmann equation, we get the 
following expression for the scale factor as: 


a{t) = Bij exp 


((i + iH?) ‘ + i2pmM, 


OQ ^2 


where B^ is the arbitrary integration constant given by: 
Bn = a/ exp 


/p / A OQBie X OQ . 

^ Vv 2P3/2 ) ^ + 12P3/2]Vfp^/ 


(5.150) 


(5.151) 
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Here af is the value of the scale factor at turnaround. 


B. Contraction 

The conclusion remains same as for ACDM model, because in this case all the extra 
terms containing 0, will be neglected at leading order approximation. 

C. Expression for work done 

The solution of the scale factor (i.e its time dependence) is same as that for hill¬ 
top and natural potential. Hence the expression for work done is same except that 
the constant will now depend on the parameters of this model. Thus the phenomenon 
of hysteresis also holds true for this model. 

For Case 3 

At late times within the window a' < a < a^ax, using Eq. (5.12) and Eq. (5.14), 
we get the integral equation of the form similar to Eq. (5.86) with the potential 
now given by the supergravity motivated Coleman-Weinberg potential. Simplihed 
analytical solutions can be obtained if we use the following field redehnition: 



which we have used earlier and take small held limiting approximation 0/Mp << 1 
in the present context. 

In this limit, we expand the exponentials upto linear order and get the resulting 
solution for the redehned held variable A as: 



(5.153) 


where B 22 is the arbitrary integration constant given by: 



(5.154) 


Here for the sake of clarity we introduce four model dependent constants R', S',U', V 
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given by; 
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The above expressions have been obtained using the assumptions that the conditions, 

(5.159) 

(5.160) 
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are satished. Since we are already in the limit in which the value of A is very small, 
hence by choosing the other parameters of the model properly, we can easily satisfy 
the above conditions. 

Further substituting the above expression into the Friedmann equation, we get 
the following expression for the scale factor as: 


a{t) = B 23 exp 


where 


B 23 = a/exp 
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(5.162) 


Here a/ is the value of the scale factor at turnaround. 

B. Contraction 

Solutions remain same as obtained before. So the conclusion also remain unchanged 
for contraction phase of the Universe. 


C. Expression for work done 
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The solution of the scale factor (i.e. its time dependence) is same as that for hilltop 
potential. Hence the expression for work done is same except that the constant will 
now depend on the parameters of this model. Thus the phenomenon of hysteresis 
also holds true for this model. 


5.7 Graphical Analysis 
5.7.1 Case I: Hilltop potential 

All the graphs in this section and in the following sections have been plotted in units 
of Mp =1, Hq = 1, c=l, where Mp is the Planck mass, Hq is the present value 
of the Hubble parameter and c is the speed of light. While performing the analysis, 
for the value of the standard cosmological constant, we have used Planck 2015 data 
[75]. The analysis for all the three potentials has been done for late times since, in 
the early times the solutions are those that we get for Randall-Sundrum single brane 
world (RSH). 

For Case 1: 

Figs. 10 and 11 show the evolution of the scale factor and kinetic term of the 
potential during contraction phase. We can draw the following conclusions from the 
plots: 


• Fig. 10(a), shows the plot of the scale factor in the small field limit for hilltop 
potential given by Eqn. (5.60). 

• While obtaining the above plot, we kept the value of the integration constants 
Bq and Bi as 1. Higher values of the integration constants only results in an 
increase in the amplitude of the scale factor, not affecting the nature of the 
graph. 

• In Figs. 10(a) and 10(b), the plot has been shown for two particular cases which 
gives interesting results. The red curve has been obtained for the case when 
Vo = 10“®Mp, p = 2, (3 = 1, Bo = 1, Bi = 1. While the blue curve has been 
obtained for the case when Vq = 10“^Mp, p = 3, [3 = —8, Bo = 1, Bi = 1. 

• While analyzing the results graphically, we have found that when the value 

of Vq is low (for example ~ 10 “®Mp), the universe undergoes expansion for 
any values of the parameters p, (3, etc. But if we make Vq > then 

expansion is possible for any values of the parameters provided Bo takes larger 


a vs t (^«Mp) for cosmological a vs t (0«Mp) for cosmological 

constant dominated Einstein gravity constant dominated Einstein gravity 




(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the scale 
factor with time during expansion phase for (j) << factor with time during expansion phase for <j) « 
Mp with Vo = p = 2, /3 = 1, Bo = Mp with Vo = p = 3, /? = -8, Bq = 

Mp, Bi = 1 . Mp, Bi = 1 . 


V(0(t)) vs t {(p«Mp) for cosmological 
constant dominated Einstein gravity 



(c) An illustration of the behavior of the potential during expansion 
phase for (j) « Mp with Vo = 10 “'*(3 = 2, p = 2, Bo = Mp . 


Figure 10. Graphical representation of the evolution of the scale factor and the poten¬ 
tial during the expansion and contraction phase for the cosmological constant dominated 
Einstein gravity for case 1. 
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a vs t (contraction) for cosmological 


vs t for cosmological 



(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the kinetic 
factor with time during contraction phase with term of the potential with time during contrac- 
B 4 = i ?5 = 10 ^. tion phase with = 5AM^. 


Figure 11. Graphical representation of the evolution of the scale factor and the kinetic 
term of the potential during the contraction phase for the cosmological constant dominated 
Einstein gravity for case 1. 


value. One such example has been shown in Fig. 10(b), where the universe 
undergoes starts contracting instead of expanding at later times. Thus the 
values of parameters for which Fig. 10(b) has been obtained, are not favourable 
for the present scenario. 

• Fig. 10(c) shows the plot of the behavior of the potential with time for small 
field hilltop potential. This graph has been obtained with the help of Eqn. 
(5.57) with parameter values Vq = lO'^M^, /3 = 2, p = 2, Bq = 1. Negative 
values of f3 and very large values of Bq does not give the required variation of 
the potential for causing expansion. Thus, such values of the parameters are 
not favorable for our model. 

• In Fig. 11(a), we have plotted Eqn. (5.70), which results in the contraction 
phase of the universe provided we choose the value of the constants accordingly. 
This plot has been obtained for B^ = 10“^, B^ = lO'^. Taking higher values 
of Bi decreases the amplitude of the scale factor, whereas higher values of B^ 
increases the amplitude of the scale factor. If we compare Fig. 8(a) with Fig. 
11 (a), we find that there occurs a net increase in the amplitude of the scale 
factor after one expansion-contraction cycle. 






a vs t for cosmological v(^(t)) vs t (^»Mp^)for cosmological 

constant dominated Einstein gravity constant dominated Einstein gravity 



(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during expansion phase for (j) » tial with time during expansion phase for cj) » 
Mp with Vb = 3xlO“^Mp, p = 3, (3 = 1, B 2 = Mp with Vq = 3xlO“^Mp, p = 3, (3=1, B 2 = 

Mp, B 3 = 1 . Mp,. 

Figure 12. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the cosmological constant dominated Einstein gravity for 
case 1. 

• In Fig. 11(b), we have plotted the kinetic term of the potential during the 
contraction phase with the help of Eqn. (5.68) using 84 ^ = 5.4Mp. Negative 
values of B 4 do not give the correct nature of the kinetic term required for 
causing contraction. The graph shows a decreasing kinetic term with time 
because this will result in contraction phase. 

• Fig. 12(a) shows the plot of the scale factor in the large field limit for hill¬ 
top potential given by Eqn. (5.66). This plot has been obtained for Vq = 
3xlO“^Mp, /I = 1, p = 3, 82 = Mp, B 3 = 1. From Eqn. (5.66), we can con¬ 
clude that negative values of (3 are not allowed in this case, since this will make 
the expression imaginary. Larger and positive values of (3 results in a decrease 
in amplitude of the scale factor Also, if we make p > 3, then the amplitude of 
the scale factor decreases to a very low value. If we make p < 2, then the evo¬ 
lution of the scale factor becomes very uneven, i.e.it remains constant for most 
part of its evolution, with a sudden increase in amplitude at the end. Thus for 
smooth and large amplitude expansion for the values of the other parameters 
that we have chosen here, the allowed range of p is 2 < p > 3. 

• Fig. 12(b) shows the plot of the behavior of the potential with time for large 
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field hilltop potential. This graph has been obtained with the help of Eqn. 
(5.62) with parameter values Vq = 3xlO“^Mp, /3 = 1, p = 3, B 2 = M^. 

For Case 2: 

In Fig. 13, we have plotted the evolution of the scale factor and the potential 
during expansion for case 2. From the figures we can draw the following conclusions: 

• Fig. 13(a) shows the variation of the scale factor with time for small field 
hilltop potential for two cases of parameter values for a scalar field dependent 
cosmological constant. 

• Fig. 13(a) has been obtained from Fq. (5.84). Detail graphical analysis show 
that an expanding universe in the present scenario is possible only if the values 
of Ai, A 2 , A 3 , A 4 lies between —10“^ to 10“^. They can take any sign provided 
their magnitudes remain within the specified range. The expansion in this case 
is almost independent of the values of (3, p and Vq. 

• In Fig. 13(a), we have shown two such cases which give rise to an expanding 
universe. 

• Fig. 13(b) and Fig. 13(c) show the variation of the potential with time during 
expansion for case 2. Detail graphical analysis have shown that keeping the 
magnitude of the parameters within the above mentioned range, if we make Ai 
negative, then expansion is possible only if A 4 is also negative or only slightly 
positive. This also requires a higher value of Vq. The other two parameters 
can be either positive or negative. Also, if the values of Ai, A 3 and A 4 are very 
close to the upper bound, then we won’t get the correct nature of the potential 
required for evolution. 

• From Figs. 13(b) and 13(c), we can conclude that higher and positive parameter 
values increases the height of the potential. But smaller and negative values of 
the constants make the potential flatter. Also, from Fig. 13(c), only a negative 
value of A 4 can produce the correct evolution of the potential if all the other 
constants are positive and small. 

For Case 3: 

In Fig. 14, we have plotted the evolution of the scale factor and the potential 
during expansion for case 3. From the figures we can draw the following conclusions: 

• Fig. 14(a) shows the variation of the scale factor with time for small field 
hilltop potential for Vq = 10“®Mp, Bi^ = M“^, A = 10“^Mp, Sig = 1. 
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- Ai =A2=A3=A4=10-®(in units of ,Mp^ ,Mp“) V((^(t)) vs t for cosmological 

-A,=A2=A3=A4=-10-^(in units of Mp^/Wp^Mp^Mp“) constant dominated Einstein gravity 



(a) An illustration of the behavior of the scale factor (b) An illustration of the behavior of the potential 
with time during expansion phase for 0 << Mp with with time during expansion phase for 0 << Mp with 
Vo = 10-^M*, p = 2, /3 = 1, Bi4 = M-2, Si 5 = Vo = p = 3, /S = 1, Bn = Mp, Ai = 

1, Ai = A2 = 10 -®M 2 , As = IQ-^Mp, A4 = A2 = A3 = A4 = 

10“® for red curve and Vb = 10“®Mp, p = 3, fi = —10“^. 

2, Bn = , Bn = 1, Ai = —A 2 = 

A 3 = A 4 = - 10-2 for blue curve 


V(0(t)) vs t for cosmological 
constant dominated Einstein gravity 



(c) An illustration of the behavior of the potential with time during expansion phase 
for (j> « Mp with Vb = 4xlO-®Afp, p = 3, /3 = 1, -B 14 = Mp, Ai = 10~^Mp, A 2 = 
10-3m|, As = 4xlO-2Mp, A 4 = -3x10-3. 


Figure 13. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the cosmological constant dominated Einstein gravity for 
case 2. 
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a vs t for DGP model 


V(^(t)) vs t for cosmological 
constant dominated Einstein gravity 




(a) An illustration of the behavior of the poten- (b) An illustration of the behavior of the poten¬ 
tial with time during expansion phase for (j) « tial with time during expansion phase for cj) « 
Mp with Vo = Bis = Mp"^, A = Mp with Vq = p = 3, /3 = 1, Bis = 

Big = 1. A/p-2, A = -10-3M4. 

Figure 14. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the cosmological constant dominated Einstein gravity for 
case 3. 


• Fig. 14(a) has been obtained from Eqn. (5.92). Detail graphical analysis show 
that an expanding universe in the present scenario is possible only if the value 
of A lies within the range —10“^ to 10“^. But as we go nearer to the upper 
bound, expansion occurs only for very small values (very close to unity) of Big. 
The expansion in this case is independent of the values of /3, p . Higher values 
of Vq increases the amplitude of the expansion. But values of Vq greater than 
10 ~^Mp do not give rise to an expanding universe. 


• Fig. 14(b) shows the variation of the hilltop potential for an expanding uni¬ 
verse for case 3 with 0 << Mp with Vq = 2x10“'^p = 3, 0 = 1, Big = 
Mp'^, A = — 10“^Mp. Graphical analysis show that within the allowed range 
of the parameter values, as discussed before, only a negative A can give the 
correct nature of the evolution of the potential when Vq is very small (around 
10“®Mp. For this case, (3 and p can take any values. But for larger values of 
Vq, we can increase the range of A and make it positive by making (3 positive. 
For example, if Vq = 10“^Mp and A = 8xlO“^Mp, then the minimum value of 
(3 required to give the correct evolution of the potential is 5. 
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5.7.2 Case II: Natural potential 
For Case 1: 


Fig. 15 show the evolution of the scale factor and the potential for the case 
of natural potential for case 1 during expansion. The following conclusions can be 
drawn from the graphs: 

• Fig. 15(a), shows the plot of the scale factor in the small held limit for natural 
potential given by Eqn. (5.99). 

• While obtaining the above plot, we kept the value of the integration constants 
Bq and i?7 as 1. Higher values of the integration constants only results in an 
increase in the amplitude of the scale factor, not affecting the nature of the 
graph. Higher values of Vq also only results in an increase in amplitude of the 
scale factor. 

• Fig. 15(b) shows the plot of the behavior of the potential with time for small 
held natural potential. This graph has been obtained with the help of Eqn. 
(5.97) with parameter values Vo = 3.6xlO“^Mp, Bq = 1, f = 8Mp. Larger 
values of the parameters gives rise to larger amplitude of expansion, but the 
nature of the graph remains the same. If we compare Fig. 15(a) with Fig. 
11 (a), we hnd that there occurs a net increase in the amplitude of the scale 
factor after one expansion-contraction cycle. 

• Fig. 15(c) shows the plot of the scale factor in the large held limit for nat¬ 
ural potential given by Eqn. (5.104). This plot has been obtained for Vq = 
10“^Mp, Bg = Mp, Bg = 1. Expansion is obtained only if the value of Vq is 
lesser than 2Mp. 

• Fig. 15(d) shows the plot of the behavior of the potential with time for large 
held natural potential. This graph has been obtained with the help of Eqn. 
(5.102) with parameter values Vq = 10“^Mp, Bg = Mp. 

For Case 2: 

In Fig. 16, we have plotted the evolution of the scale factor and the potential 
during expansion for case 2. From the hgures we can draw the following conclusions: 

• Fig. 16(a) shows the variation of the scale factor with time for small held 
natural potential for two cases of parameter values for a scalar held dependent 
cosmological constant. 
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a vs t for cosmological 
constant dominated Elsntein gravity 


V(^(t)) vs t (^«f) for cosmological 
constant dominated Elsntein gravity 




(a) An illustration of the behavior of the scale 
factor with time during expansion phase for (j) « 
f with Vo = 10-®M4, Be = 1, = 1. 

a VS t for cosmological 
constant domlanted Einstein gravity 


(b) An illustration of the behavior of the poten¬ 
tial during expansion phase for (j) « f with 
Vo = 3.6x10-3m4, Be = l, / = 8Mp . 

V(0(t)) vs t (0»f) for cosmological 
constant domlanted Einstein gravity 




(c) An illustration of the behavior of the scale 
factor with time during expansion phase for (j) » 
f with K) = Bs = Mp, Bo = 1. 


(d) An illustration of the behavior of the poten¬ 
tial with time during expansion phase for (j) » f 
with Vo = Bs = Mp,. 


Figure 15. Graphical representation of the evolution of the scale factor and the poten¬ 
tial during the expansion and contraction phase for the cosmological constant dominated 
Einstein gravity for case 1. 
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-Ai=A 2 =A 3 =A 4 =- 10 '^(in units of Mp^Mp^Mp^Mp‘') 
— Ai=A2=A3=A4=3*10'^(in units of Mp^Mp^Mp^Mp‘’) 



V(^(t))vs t (^«f) for cosmological 
constant dominated Einstein gravity 





(a) An illustration of the behavior of the scale factor (b) An illustration of the behavior of the potential 
with time during expansion phase for (j) « f with with time during expansion phase for (p « f with 
Vo = 10-^M^, f = 2Mp, Bie = Mp^, Bir = 1, Ai = Vo = f = 4Mp, Big = Mp'^, Ai = 

A 2 = -10-^M^, As = -10-^Mp, A 4 = A 2 = 1.2xlO-^M^, A 3 = 2.4x10-3^4, A 4 = 

— 10-3 for red curve and Vq = 10 -*Mp, / = 6.5x10-3. 

6Mp, Big = M-3, Bi7 = 1, Ai = 3x10-3m|, A 2 = 

3xlO-3Af3, A 3 = 3x 10—3Mp, A 4 = 3x10—3 for blue 
dashed curve . 


Figure 16. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the cosmological constant dominated Einstein gravity for 
case 2. 


• Fig. 16(a) has been obtained from Eq. (5.118). Detail graphical analysis show 
that an expanding universe in the present scenario is possible only if the values 
of Ai, A2, A3, A4 lies between —10“^ to 10“^ similar to the case for hilltop 
potential. They can take any sign provided their magnitudes remain within 
the specihed range. But expansion is possible only if / > IMp and value of Biq 
is very close to unity. If we make the value of Biq too large, then expansion is 
possible eithe rif all the other constants are positive except A2 or if only A2 is 
positive and any two or all three other constants are negative. Increase in the 
value of Vq only increases the amplitude of scale factor. 

• In Fig. 16(a), we have shown two such cases which give rise to an expanding 
universe. 

• Fig. 16(b) shows the variation of the potential with time during expansion for 
case 2. Detail graphical analysis have shown that keeping the magnitude of 
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the parameters within the above mentioned range, expansion is possible only 
if As > SxlO-^Mp, A4 > 5x10-3 and Vq > 

• Higher and positive parameter values increases the height of the potential. 

For Case 3: 


In Fig. 17, we have plotted the evolution of the scale factor and the potential 


a vs t for cosmological 
constant dominated Einstein gravity 



V(0(t)) vs t (0«f)for cosmological 
constant dominated Einstein gravity 



(a) An illustration of the behavior of the po- (b) An illustration of the behavior of the poten- 
tential with time during expansion phase for tial with time during expansion phase for </><< / 


10 


« / with H, = 10-8M4, B 20 = M-2, 

iy±p , 


S21 = 1. 


A = with Vo = 4x10-3m4, / = 1 , 
- 10 - 


B 20 — M A — 




Figure 17. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the cosmological constant dominated Einstein gravity for 
case 3. 


during expansion for case 3. From the figures we can draw the following conclusions; 

• Fig. 17(a) shows the variation of the scale factor with time for small field 

natural potential for Vq = lO-^M^, B 20 = A = B 21 = 1. 

• Fig. 17(a) has been obtained from Eq. (5.122). Detail graphical analysis 
show that an expanding universe in the present scenario is possible only if the 
value of A lies within the range —10“^ to 10“^. But if / < IMp, expansion is 
possible only if B 20 is very close to unity and A > —2.8xlO-^Mp. For / > IMp, 
expansion is possible if A lies within the range — IQ-^Mp to IQ-^Mp 
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• Fig. 17(b) shows the variation of the natural potential for an expanding uni¬ 
verse for case 3 with 0 << / with Vq = 4xlO“^Mp, / = 1, i? 2 o = A = 

— Fig. 17(b) has been obtained with the help of Eqn. (5.120). Graph¬ 

ical analysis show that within the allowed range of the parameter values, as 
discussed before, larger values of Vq and A results in steeper fall of the potential 
and larger potential height. 

5.7.3 Case III: Coleman-Weinberg potential 
For Case 1: 


Fig. 18 show the evolution of the scale factor and the potential for the case of 
supergravity potential for case 1 during expansion. The following conclusions can be 
drawn from the graphs: 

• Fig. 18(a), shows the plot of the scale factor in the small held limit for su¬ 
pergravity potential given by Eqn. (5.131) with the parameter values Vq = 
2xlO-^M*, Bio = Mp, Bn = 1, a = 0.5, (3 = 1. 

• While obtaining the above plot, we kept the value of the integration constants 
Bio and Bn as 1. Higher values of the integration constants only results in 
an increase in the amplitude of the scale factor, not affecting the nature of the 
graph. Higher values of Vq also only results in an increase in amplitude of the 
scale factor. The evolution of the scale factor is almost independent of a. Large 
increase in amplitude of the scale factor occurs for large value of (3 irrespective 
of its sign. 

• Fig. 18(b) shows the plot of the behavior of the potential with time for small 
held supergravity potential. This graph has been obtained with the help of Eqn. 
(5.127) with parameter values Vq = 4.6xlO“^Mp, Bio = 0.5, a = 0.18, (3 = 2. 
Larger values of the parameters gives rise to larger amplitude of expansion, but 
the nature of the graph remains the same. If we compare Fig. 18(a) with Fig. 
11 (a), we hnd that there occurs a net increase in the amplitude of the scale 
factor after one expansion-contraction cycle. 

• Fig. 18(c) shows the plot of the scale factor in the large held limit for su¬ 
pergravity potential given by Eqn. (5.139). This plot has been obtained for 
Vo = 10“®Mp, Bi 2 = 10“®Mp, Bio = 1, (3 = 7. Proper expansion is obtained 
only if the value of Bn is << 1 and the value of (3 is large. Larger value of Vq 
results in nearly constant expansion amplitude with a huge increase towards 
the later times. Expansion is possible for both negative and positive values of 
(3. 
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a vs t for cosmological V(^(t)) vs t (0«Mp) for cosmological 

constant dominated Einstein gravity constant dominated Einstein gravity 



(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during expansion phase for << tial during expansion phase for (j) « Mp with 
Mp with Vo = 2x10-3m4, Bio = Mp, Bn = Vo = d.GxlO-^M^, B^ = 0.5, a = 0.18, /3 = 2 . 
1, a = 0.5, 13 = 1. 


a vs t for cosmological 
constant dominated Einstein gravity 



V(^(t)) vs t {<p»Mp}for cosmological 
constant dominated Einstein gravity 



(c) An illustration of the behavior of the scale (d) An illustration of the behavior of the poten- 
factor with time during expansion phase for (j) » tial with time during expansion phase for (j) » 


Mp with Vo = 10’ 
1, /3 = 7. 


-8M4, 


Bi 2 = 10-®Mp, Bi 3 = Mp with Vo = 8x10-3m4, B 12 = IQ-^M, 


pi 


1, /3 = -0.01. 


Figure 18. Graphical representation of the evolution of the scale factor and the poten¬ 
tial during the expansion and contraction phase for the cosmological constant dominated 
Einstein gravity for case 1. 
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• Fig. 18(d) shows the plot of the behavior of the potential with time for large 
held supergravity potential. This graph has been obtained with the help of Eqn. 
(5.137) with parameter values Vq = 8xlO“^Mp, = 10“®Mp, a = 1, /S = 
—0.01. Proper evolution of the potential which will give rise to expanding 
universe is possible only if /3 < 0 and B 12 is << 1. Larger values of a makes 
the potential fall more steeply. Larger values of Vq mildly increase the height 
of the potential. 

For Case 2: 

In Fig. 19, we have plotted the evolution of the scale factor and the potential 


a vs t for cosmological 
constant dominated Einstein garvlty 



V(0(t)) vs t (0«iWp) for cosmological 
constant dominated Einstein garvlty 



(a) An illustration of the behavior of the scale factor (b) An illustration of the behavior of the potential with 
with time during expansion phase for cp « Mp with time during expansion phase for (p << Mp with Vq = 
Vo = Bio = Mp'^, Bit = 1 , Ai = a = Bm = IQ-^Mp^, Ai = 

2.6x10-2m|, A 2 = A 3 = A 4 = A 2 = 7.9x10-2m|, A 3 = T.gxlQ-^Mp, A 4 = 

8.8x10-2, Q = 0.13, 0 = -0.4 . 8.4x10-1, 0 = -1.85. 


Figure 19. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the cosmological constant dominated Einstein gravity for 
case 2. 


during expansion for case 2. From the figures we can draw the following conclusions: 

• Fig. 19(a) shows the variation of the scale factor with time for small held 
supergravity potential for case 2 with Vq = 1.3xlO“^Mp, Biq = Bn = 

1 , Ai = 2.6 x10-2m3, A 2 = - 10 -^M 2 , As = -IQ-^Mp, A 4 = 8 . 8 x 10 - 2 , a = 

0.13, (3 = -0.4 
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• Fig. 19(a) has been obtained from Eqn. (5.150). Detail graphical analysis 
show that an expanding universe in the present scenario is possible only if the 
values of Ai, A 2 , A 3 , A 4 lies between —10“^ to 10“^. They can take any sign 
provided their magnitudes remain within the specihed range. But expansion is 
possible for /3 > 0 only if any two of the constants (Ai, A 2 , A 3 , A 4 ) are positive 
and the other two take negative values. For /3 < 0, exoansion is possible for 
any possible combination of signs of the constants. Large values of a results in 
larger amplitude of expansion, but the nature remains the same. 

• Fig. 19(b) shows the variation of the potential with time during expansion 
for case 2 obtained with the help of Eqn. (5.142) and with parameter val¬ 
ues Vo = 2 . 1 xlO- 2 Mp^ a = 10 -^ Biq = IQ-^M-^, Ai = -lO'^M^, A 2 = 
7.9xlO“^Mp, A 3 = 7.9xlO“^Mp, A 4 = 8.4x10“^, (3 = —1.85. Detail graphical 
analysis have shown that keeping the magnitude of the parameters within the 
above mentioned range, potential evolves with the required nature to cause 
expansion only if /3 < 0 and a < 1 . 

• Higher and positive parameter values increases the height of the potential. 

For Case 3: 

In Fig. 20, we have plotted the evolution of the scale factor and the potential 
during expansion for case 3. From the hgures we can draw the following conclusions; 

• Fig. 20(a) shows the variation of the scale factor with time for small held su¬ 
pergravity potential for Vq = 1.5xlO“^Mp, B 22 = A = 4xlO“^Mp, H 23 = 

1, a = 0.96, 13 = 1. 

• Fig. 20(a) has been obtained from Eqn. (5.161). Detail graphical analysis 
show that an expanding universe in the present scenario is possible only if the 
value of A lies within the range —10“^ to 10“^. If we decreas the value of Vq 
to 10 “’^Mp(say), then expansion is possible only for values of B 22 « 
Expansion is possible for both positive and negative values of /3, however small 
positive values of f3 are more favourable for expansion. 

• Figs. 20(b) and 20(c) show the variation of the supergravity potential for an 
expanding universe for case 3 with 0 << Mp for two sets of parameter values, 
Do = SxlO-^M^, a = 0.8, B 22 = M-\ f3 = 0.15, H 23 = 1, A = -O.OxlO-^M^ 
and Vo = 1 . 2 xlO“^Mp, a = 10“®, B 22 = 10 “®M“^, (3 = —0.7, B 23 = 1 , A = 
4.3xlO“^Mp respectively. Figs. 20(b) and 20(c) have been obtained with the 
help of Eqn. (5.153). Graphical analysis show that within the allowed range 
of the parameter values, as discussed before, for values of A near the upper 
bound, expansion is possible only for negative values of (3 and values of other 


100 


a vs t for cosmological V(^)(t)) vs t (0«Mp) for cosmological 

constant domianted Einstein gravity constant domianted Einstein gravity 




(a) An illustration of the behavior of the poten- (b) An illustration of the behavior of the poten¬ 
tial with time during expansion phase for (p << tial with time during expansion phase for <p « 
Mp with Vo = 1.5x10-3m4, B 22 = M-"^, A = Mp with Vq = SxlQ-^M^, a = 0.8, B 22 = 
4x10-3m4, B 23 = 1, a = 0.96, 13 = 1. M-^, /3 = 0.15, B 23 = 1, A = -S.GxlQ-^M^. 


> 


V(0(t)) vs t {<p«Mp) for cosmological 
constant domianted Einstein gravity 



(c) An illustration of the behavior of the potential with time during ex¬ 
pansion phase for p « Mp with Vq = 1.2xlO~^Mp, a = 10“®, B 22 = 
10-®Mp-2, j3 = -0.7, B 23 = 1, A = 4.3x10-3m4. 


Figure 20. Graphical representation of the evolution of the scale factor and the potential 
during the expansion phase for the cosmological constant dominated Einstein gravity for 
case 3. 
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parameters {a, Vq, -B 22 ) much less than unity (~ 10“®). For positive value of 
expansion is possible if the value of Vq > and that of a is close to 

unity. In this case, in order to get an expanding universe, A must have small 
negative values. Larger the magnitude of jS, steeper the fall of the potential. 


6 Hysteresis from Loop Quantum gravity (LQG) model 


Loop Quantum gravity (LQG) is one of the candidate theories of quantum gravity 
which successfully resolves the big bang singularity. It is primarily based on the 
quantum geometry effects of loop quantum gravity. 

For the spatially flat case when the curvature parameter k = 0, the modihed 
Friedmann equations in this model are given by: 


H + H^ 





( 6 . 1 ) 

( 6 , 2 ) 


which is exactly same as that we get for RSII brane world cosmology (see Appendix). 
Since this analysis has already been done in detail by [11], we have not shown it here 
again. However the results of this analysis have been quoted in the Appendix for 
reference. Hence the rest of the conclusions remain same as that obtained for RSII 
brane world model. 

Modihed Friedmann equations for non hat case i.e. for A: ^ 0 in this model 
[76, 77] are given by: 


= -] = 


^ {P-Pi) (— {P2-P)] , 


3M2 


H + H^ = - = 


6 M 2 


(p + 3p) + 


Pc 

2 


3M2 Vpc 


— + kx] (p+\p 


+ 


kx f P 


— + kx 


\p, 

where we introduce two new constants pi and p 2 are dehned as: 

-3kxM^ 




P2 : = Pc(l - kx) ■ 


= -kxPa 


(6.3) 


(6.4) 


(6.5) 

( 6 . 6 ) 


Here the critical density pc, the model parameters A, y and C, can be expressed within 
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LQG setup as [76, 77]: 


(6.7) 

(6.8) 


Pc = 




7^ A ’ 

A = Ay/Sn'ylp = Ji^p, 


sin^/i - (l + 7 ^) /i^, 


X = 


2 -2 
-7 , 


for A: = +1 


for k = — 1 . 


( = sin^ ^ ~ 2 


(6.9) 


( 6 . 10 ) 


where p is the triad (which without any loss of generality will be chosen with positive 
orientation) and the Barbero-Immirzi parameter 7 can be hxed by computing the 
black hole entropy in LQG. Now let us concentrate on the classical limit, A —)■ 0, 
from LQG setup one can write: 


, 2-2 
X ^ -7 /i , 

3kM^ 

Pi -) 

p 


(pc- - p) 
Pc 


-»1. 


( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 


Thus we recover the Friedmann equation in classical GR in the limit A —)■ 0. In the 
next subsections we will discuss the detailed cosmological consequences as well the 
phenomena of cosmological hysteresis from this model. 


6.1 Condition for bounce 

Bounce occurs when the following condition holds good: 

P = P 2 , (6.14) 

where the universe reaches its minimum radius a^m and maximum density [77] is 
achieved: 

Pmax = P 2 lamin^ Pc- (6.15) 

Therefore condition for bounce is given by: 

Ph = Pc- (6.16) 

The mass at bounce (neglecting the constant factor) is given by: 

Mb = ptal = Peal- (6.17) 
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Therefore an infinitesimal variation in mass at bounce can be expressed as: 


5M ^ 


(6.18) 


By setting the following constraint condition; 


SM = -6W = 


pdV, 


(6.19) 


we get the expression for change in amplitude of the scale factor at each successive 
cycle as: 

n 3 _ - ^ ('0 20 ) 


di^Oimin) _ ^ pdV 


1 

pc 


Thus, the amplitude of the scale factor increases iff: 


pdV < 0 


( 6 , 21 ) 


as the critical density pc > 0 always. We also observe that now the increase in 
scale factor depends on the critical density Pc but is independent of the curvature 
parameter k. Hence the result holds good for k = ±1 i.e. for both open and closed 
universe. 


6.2 Condition for acceleration 


From Eq. (6.4), at bounce the condition for acceleration is given by: 


d 

a 


6 M 2 


{pc + Spfe) + (1 + kx) 


6 M 2 

kx 

72 A 


(1 + kx) 


3 

Pc + ^Pb 

2 Ck 

72 A 


, 1 

kx + ^ 


( 6 . 22 ) 


Thus it implies that whether the condition for acceleration violates the energy condi¬ 
tion, now depends upon the values of different parameters of the LQG model present 
in the above expression. 

For different values of the curvature parameter k we get the following constraint 
conditions for acceleration at bounce as: 


2xM\l + x) 2CM^ Pc(3 + 4x) 
72 A 1 + 2x + 72 A 3(1 + 2x) 


for k = +1 


Pb> < 


2xM^ (1 - x) 2CM^ p,(3 - 4x) 

^ 72 A l-2x 72 A 3(1-2x) 


(6.23) 


In Figs. 
LQG model. 


21 , we have shown the phenomena of bounce and acceleration in the 
We can draw the following conclusions from the above hgures: 
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H vs p for LQG model 
-Fork=1 
-Fork=-1 
-ForkzO 



(a) An illustration of the bouncing condition for 
a universe withfc = l,pc = lOM^ for red curve 
and fc = 1, w = 1/3, pc = lOM^ for blue curve 
and fc = —1, w = l/3, Pc = lOM"* for green 
curve. 


avspforLQGniodel(k=-1) 



a vs p for LQG model (k=0) 



(b) An illustration of the acceleration condition 
at the time of bounce for a universe with an equa¬ 
tion of state w = 0, k = 0, Pc = lOA/^. 


a VS p for LQG model (k=1) 



H(Mp) H(M,) 

(c) An illustration of the acceleration condition (d) An illustration of the acceleration condition 
at the time of bounce for a universe with an equa- at the time of bounce for a universe with an equa¬ 
tion of state w = 0, k = —1, pc = lOM^ tion of state w = 0, k = 1, Pc = lOA/^. 

Figure 21. Graphical representation of the phenomena of bounce and acceleration for 
LQG model. 
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• In Fig. 21(a), we have plotted the r.h.s of Eq. (6.3) using the relation p = 
0-3(1+’^) with /c = l,pc = lOM^ for red curve and A; = 1, w = 1/3, pc = lOM"^ 
for blue curve and k = —1, w = 1/3, pc = lOM^ for green curve. 

• We have used w = 1/3, since we require a soft equation of state for causing the 
acceleration and expansion. The case tc = 0 causes bounce equally well. The 
value of Pc chosen is arbitrary. 

• From Fig. 21(a), we can also say that the behavior of the Hubble parameter 
w.r.t. density is nearly same for any value of k. 

• From Fig. 21(a), we get the bounce at p = pb = lOM^ for all the three cases. 

• Figs. 21(b) to 21(d), we have shown the necessary condition of acceleration 
(a > 0) at the time of bounce for /c = 0, —1, 1. Here we have plotted 
the r.h.s of Eqns. (6.4) and (6.3). These plots have also been obtained for 
m; = 0, Pc = lOMh 

• Thus from Fig. 21 we can conclude that for this model, bounce is possible for 
closed, open and flat universe. 

6.3 Condition for turnaround 

Turnaround or re-collapse occurs when the following criteria is achieved: 

p = pi, (6.24) 


when the universe reaches its maximum radius amax and minimum density can be 
written as [77]; 

3A:M2 

Pi = Pmin ~ ^ -• (6.25) 

( 1 ^ 

^max 

Therefore, in place of Eq. (6.20), we get: 

Sttmax = ( 6 - 26 ) 

Unlike for the bounce case, the condition for an increase in expansion amplitude at 
turnaround now depends on the curvature parameter only (apart from the sign of 
the work done). Substituting the values of curvature parameter Eq. (6.26) can be 
recast in the following form: 





for k = +1 


6 a 


max 


< 


(6.27) 


V 



for k = —1. 
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Let us explicitly mention the two possible physical outcomes from Eq. (6.29) appear¬ 
ing in the present context: 

• For k = -|-1, in order to get an increase in the amplitude of the scale factor 
after each successive cycle, we need f pdV < 0. 

• For k = —1, in order to get an increase in the amplitude of the scale factor 
after each successive cycle, we need f pdV > 0. 

Rest of all the conclusions remain same as that we obtained for the case of bounce 
from LQG setup. 


6.4 Condition for deceleration 

Using the dehnition of pi = —kxPc from Eq. (6.4), at turnaround time scale the 
condition for acceleration is given by: 


pi + 3pt > 


6M^x/c 


(6.28) 


Further substituting the expression for pi for different values of the curvature pa¬ 
rameter k in non-flat case we get the following constraint conditions for acceleration 
at bounce as: 






for k = +1 


Pt= < 


(6.29) 






for k = — 1 . 


Just like the case of acceleration within LQG setup, the condition for deceleration at 
turnaround also depends on the LQG model parameters. Additionally it is important 
to note that whether this constraint violates the energy condition or not, solely 
governed by the numerical values of the constants. 

In Fig. 22, we have shown the phenomena of turnaround and deceleration in 
LQG model. We can draw the following conclusions from the above hgures: 


• In Fig. 22(a), we have plotted the r.h.s of Eqns. (6.3), using the relation 

p = with fc = —1, in = l, Pc= lOM^ for blue curve and fc = 1, w = 

1, Pc = lOM^ for green curve. For this case we have used w = 1, since we 
require a stiff equation of state for causing contraction and deceleration. 

• From Fig. 22(a), we get the turnaround at p = pt = 0 for k = 0, —1 and at 
p = Pt = respectively. 
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H vs p for LQG model 
-ForkzO 
-Fork=-1 
Fork=1 



P(in<) 


(a) An illustration of the turnaround condition 
for a universe with k = 0, Pc = lOM'^ for red 
curve, k = —1, re = 1, Pc = lOM"^ for blue curve 
and k = 1, w = 1, Pc = lOM^ for green curve . 



a VS H for LQG model (k=0) 




(b) An illustration of the deceleration condition 
at turnaround for a universe with /c = 0, Pc = 
IOMA 


a VS H for LQG model (k:1) 





(c) An illustration of the deceleration condition (d) An illustration of the deceleration condition 
at turnaround for a universe with an equation of at turnaround for a universe with an equation of 
state w = 1, k = —l,pc = lOM"'. state w = 1, k = 1, pc = 


Figure 22. Graphical representation of the phenomena of turnaround and deceleration 
for LQG model. 
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• Figs. 22(b), 22(c) and 22(d) show the necessary condition of deceleration (a < 
0) at the time of tnrnaround. Here we have plotted the r.h.s of Eqns. (6.4)and 
(6.3). This plot has also been obtained for the same parameter values as the 
earlier graph. 

• Thus Fig. 22 shows graphically that the phenomenon of turnaround is possible 
for the DGP model having w = 1. 

6.5 Evaluation of work done in one cycle 

The expression for the total work done is same as that given by Eq. (4.52) in DGP 
model. But in order to get an expression similar to Eq. (4.63), we see that just like 
in DGP model, here also the Friedmann equations given by Eq. (6.3) and Eq. (6.4), 
which are highly complicated. Hence to get a analytical solution, just like in the 
case of DGP model, here also we can use the early and late time approximations to 
Eq. (6.3). 

At early time, p » pi, hence we can neglect pi in Eq. (6.3). Then we get: 

2 


3M2 


^ p{ — i.p2-p) 

\Pc 


Solving the above equation for p and using the following relation [76]; 

Pi+ P2= Pc 

we get the following two-fold solution for the energy density: 

Pc± JpI- 4:i3M‘^p^ + P 1 P 2 ) 

P= -^^-• 


(6.30) 


(6.31) 


(6.32) 


But we know that at the time of bounce {H = 0), the solution of this equation is: 


P — P2 — Pc 


(6.33) 


which we approximately get by setting if = 0 in Eq. (6.32) and taking the positive 
signature in the solution. Therefore the physically acceptable solution for the density 
p is given by: 

Pc+ \/pc -4(3M2pc(f)%PiP2) 


Pphys — P — 


(6.34) 


Using the energy conservation equation (continuity equation); 

p-|-3ii(p-|-p) = 0 (6.35) 

and Eq. (6.30), we get the following form of the second Friedmann equation written 

2 1 


as: 


a 

a 


2M2p, 


■(2p-p2)(p + p), 


(6.36) 
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from which we get the following expression for the pressure p in early times as: 



2MV, 

(2p-p2) Va/ 



(6.37) 


Further using Eq. (6.34) into Eq. (6.37), we get the following expression for the 
second and third integral of Eq. (4.53) as: 



V / 


a^adt. (6.38) 


At late times, p << p 2 , hence the following condition holds good: 


{P2 - P) 
Pc 


1 . 


(6.39) 


Hence using pi = 3kM‘^ Eq. (6.3) reduces to standard Friedmann equation 
given by: 

Using the standard second Friedmann equation for acceleration and the above equa¬ 
tion we get the following expression for pressure p as: 



U(0) = 




(6.41) 


Therefore the expressions for the first and last integral in Eq. (4.53) are given by 




pdV 


j 3M^ (2dda -|- -|- kd) dt 


(6.42) 


Therefore the complete expression of the work done in one single expansion-contraction 


- no - 
















cycle is obtained by substituting Eq. (6.38) and Eq. (6.42) into Eq. (4.53) i.e. 


/ 


pdV 



(6.43) 


From the above expressions we can conclude that, at early times, work done depends 
on the LQG model parameters like the critical density Pc, but depends only on the 
curvature parameter k at late times. 


6.6 Semi-analytical analysis for cosmological potentials 

As we have done for the other models, here also we will consider the cases of three 
different potentials as mentioned earlier. In this subsection we will denote Planck 
mass by Mp. 

6.6.1 Case I: Hilltop potential 

Since we need to substitute the expression for the scale factor into Eq. (4.64), we 
need to find separate expressions for the scale factor for both early and late times 
during expansion and contraction phases of the Universe. 


A. Expansion 
i) Early time 


At early time within the interval amin < a < a' we will do the analysis for the 
case when the curvature parameter k ^ 0, because as we have already seen, that the 
curvature term is necessary for causing the turnaround at late times. To serve this 
purpose we can use the approximate Friedmann equation given by Eq. (6.30), where 
p is now given by the hilltop potential. Then substituting the resulting expression 
for the Hubble parameter H from LQG setup into Eq. (4.67), we get an integral 
equation of the following form; 



(6.44) 
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Now it is important to note that the exact analytical solution of the left hand side of 
the above integral equation is given in the Appendix for RSII model. But for LQG 
setup to get analytical result, we consider the following field redefinition; 


0 

Mp 



(6.45) 


where now we will solve for A instead of 0. Simplified analytical expression for A was 
possible only if we consider the small field limiting case i.e. (p/Mp « 1, which has 
been elaborately discussed below. 

For this case we can expand the exponentials upto linear order and apply the 
following constraint conditions within LQG setup as: 


P/3^ 

(1 + (3) 
Vo{l + /3{l+p\)) 

Pc 


« 1 , 
« 1 . 


(6.46) 

(6.47) 


Since we are in the small A limit, these conditions are possible provided we choose 
the values of the parameters of the model accordingly. The solution for A is now 
given by: 


A 


+ ^o) 

Q 


(6.48) 


where we introduce two constants Q and Eq given by: 


Q 

Eo 



1/2 

(1 + /5)'/' 


QAj + 


VoU 

3M2 





(6.49) 

(6.50) 


Here A* is the value at bounce which corresponds to the time t = ti. 

Further substituting the expression for A back into the approximated version of 
the Friedmann equation, we get the following expression for scale factor as: 


a{t) = El exp 
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where we introduce a new constant Ei given by: 


El = Oj exp 
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3M2 
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(6.52) 


ii) Late time: 


At late times within the interval a' < a < amax, the Friedmann equation is given 
by Eq. (6.40), which is the standard Friedmann equation appearing in the context 
of classical GR. Thus, in order to generate the condition for turnaround, we need to 
take the curvature parameter k ^ 0 for this specihc case. Thus, using the standard 
classical GR version of the Friedmann acceleration equation and considering that 
we have been solving for the cases when the contribution from the canonical kinetic 
term is smaller compared to the potential i.e. p ~ —p ~ V (0) for expansion and the 
contribution from the canonical kinetic term is larger compared to the potential i.e. 
p K. p K. 0^/2 for contraction, we get a differential equation for the scale factor of 
the form: 


a 

a 



(6.53) 


and solving this equation we get the two-fold expressions for the scale factor of the 
following form: 


a{t) = exp 


(y/^t + E') ± e 2^*+®2 _ 2k 


(6.54) 


and the form of E 2 and E 2 can be evaluated from the appropriate boundary condi¬ 
tions. 

Substituting the above expression back into the Friedmann equation, we get the 
expression for the potential as 


1/(0) = ^ 2(1 + 2fce(®2+2v^b) + g(-2E'-2e(^2+2V^hC+4A:-2v/E^t)_ 


(6.55) 


B. Contraction 

From Eq. (4.68), we see that for the contraction phase, under the approximation 
which we have assumed, the analysis becomes independent of the choice of the po¬ 
tential. Hence, the hnal results for the LQG model holds good for any form of the 
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potential. 


i) Early time: 


At eaarly times within the interval amin < a < a', following the same analysis as we 
have already done for the case of expansion, with the new expression for density, the 
integral equation is now given by: 


dcj) 


2 \ 1/2 

_ 1 V 

v^Mp V 2p, 


= — dt. 


3</>2 


Solving the above integrals, we get the following expression for 0 as: 


2Ml 


3(t + + g 


where we introduce a new constant dehned as: 


1 /2M^ 

Es = -Ut^ . 

Here 0* is the value of the derivative of the scalar held at bounce. 
Integrating the above equation we get the solution for 0 as 


where 


(6.56) 


(6.57) 


(6.58) 


0(^) = 2y|^tan \VQ{E3 + t)^jMp) + E'^ (6.59) 

E'^ = (l)i- 2^^^tcin~^{V6{E3 + ti)^/^c/Mp (6.60) 


Further using the above expression in the Friedmann equation, we get the ex¬ 
pression for scale factor as: 


a(t) = E 4 ^ exp 


(-£, + i) ln[M| + 3(E, - i)VJ V motA 


6 Mr, 


{E3-t)2pc 


"py M2+3(£;3-t)2pc 
where we introduce a new constant E 4 , dehned as: 


E4, = ttiCXp 


(6.61) 


(-E, + U) + 3(£3 - ti?Pc]J MiJlt-t.) 


3 ^i) pc 


6Mr 


(Es-tippc 


^P\J M^+3{E3-UPp, 

where Oj is the value of the scale factor at the time of bounce fj. 


(6.62) 
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ii) Late time: 


At late times within the interval a' < a < Qmax, the Friedmann equation given 
by Eq. (6.40), which is the standard Friedmann equation as appearing in the context 
of classical GR. Thus, in order to generate the condition for turnaround, we need to 
take k ^ 0 for this case. Thus, using the standard Friedmann acceleration equation 
and considering that we have been solving for the cases when the p ~ —p ~ V(</>) 
for expansion and p k. p k. 0^/2 for contraction, we get a differential equation for 
the scale factor of the form; 


2 a^ 2k 

a, H-1-— 0. 

a a 


(6.63) 


Now solving the above differential equation for a(t), we get the following expression 
for the scale factor as: 


a{t) = InverseFunction 


-k 


1 

+ ProductLog[^- 


{t + E,) 


(6.64) 


where Eg, Ej and Ki are the arbitrary integration constants whose values are deter¬ 
mined from the appropriate choice of boundary conditions 

C. Expression for work done 

The expression for work done is given by the sum of Eq. (6.42) and Eq. (6.38). 
Evaluation of the integrals using the corresponding scale factor for small field limit¬ 
ing situation as have been calculated in this section gives 



fFi + W2, 


(6.65) 


^ Here InverseFunction (/)(a:) represents the inverse of the function f{x). ProductLog or Lambert 
w function or omega function gives the principal solution for w in z = we'", for any complex number 
z. By implicit differentiation, one can show that all branches of w satisfy the differential equation: 
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where Wi and W 2 is defined as: 


Wi = k InverseFunction 
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(6.67) 


where di...d 2 o are constants that depends on the parameters of the model whose ex¬ 
plicit forms have been given in the appendix. 


6.6.2 Case II: Natural potential 
A. Expansion 

i) Early time: 

^Here Ei(z) represents the exponential integral function, where Ei(z) = f ^-^dz and Erfi repre¬ 
sents the imaginary error function given by, Erfi( 0 ) = 
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At early time, within the interval Omm < a < a', we will do the analysis for the 
case when kT^O, because as we have seen, the curvature term is necessary for causing 
the turnaround at late times. We can use the approximate Friedmann equation given 
by Eq. (6.30), where p is now given by the natural potential. Then substituting the 
resulting expression of H into Eq. (4.67), we get an integral equation given by 


/■ / jy Y'" (hW!(|)y^ 

J sinih 



dt (6.68) 


The exact solutions of the above integrals are given in the Appendix. In order to 
simplify the expressions, we solve the above integrals for two limiting cases: 


a) 0// « 1: 


For this case we take small argument approximations of the trigonometric functions 
after which we get the solution for 0 as 


4>{t) 

f 


Egexp 



(6.69) 


where 



(6.70) 


Here 0* is the value of 0 at the initial time of bounce fj. 

Substituting the expression for 0 back into the Friedmann equation, we get the 
expression for scale factor as 


a{t) = Egexp 




where 




Here a* is the value of the scale factor at bounce. 


(6.71) 


(6.72) 


b) 0// » 1: 


For this case, considering the fact that for large argument the cosine function is 
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small and the sine function can be approximated to unity, the solution for 0 in this 
case simplifies to 
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(6.74) 


Here is the value of the scalar field at the time of bounce. 

Substituting back the expression for 0 into the Friedmann equation and inte¬ 
grating, we get the expression for the scale factor as 


a(t) = Eiiexp 
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2 . ^ 2(3i?io/^-I-tjV'o)pc 1 
/ sm(—^==-)+4fip, 

3./^P(V,-2pP 


ii) Late time: 


At late times, within the interval a' < a < amax the analysis is independent of 
the choice of the potential. Consequently the results shown in this subsection hold 
good for natural potential. 

B. Contraction 

As has already been mentioned that while performing the analysis for expansion. 
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the conclusions for contraction phase is independent of any potential, hence remains 
same. 


C. Expression for work done 

The expression for work done is given by the sum of Eq. (6.42) and Eq. (6.38). 
Evaluation of the integrals using the corresponding scale factor analytically has been 
possible for small held limiting situation, which as have been calculated in this section 
as shown below: 


pdV = Di + D 2 , 


(6.77) 


where Di and D 2 is dehned as: 
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(6.79) 


where di...d 24 are constants that depends on the parameters of the model whose 
explicit forms have been given in the appendix. Hence from our analysis we can see 
that the phenomenon of hysteresis is achieved due to non-zero work done for natural 


119 

































potential in LGG setup. 


6.6.3 Case III: Coleman-Weinberg potential 
A. Expansion 

i) Early time: 


At early time within the interval, Omm < a < to do the analysis we use non¬ 
vanishing curvature parameter 7 ^ 0 , because as we have seen, the curvature term 
is necessary for causing the turnaround at late times. We can use the approximated 
version of the Friedmann equation as given by Eq. (6.30), where the energy density 
p is now characterized by the Coleman-Weinberg potential. Then substituting the 
resulting expression of H into Eq. (4.67), we get an integral equation given by 




1 + 
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—_ Yo_ f pf 

~ 3M2 J 


(6.80) 

To compute the left hand side of the above integral equation, we consider the follow¬ 
ing field redefinition; 

as we had done for the case of hilltop potential. Here now we will solve for the 
transformed or redefined field A instead of 0. Simplified analytical expressions for 
the above integral was possible only in the small field limiting case, 0 /Mp << 1 , 
which has been elaborately discussed below. 

For this case we can expand the exponentials upto linear order after which we 
get the following solution for A as: 


A = 


~ ^ 12 ) (4a + /3) 



1/2 ’ 


(6.82) 


where E 12 is the arbitrary integration constant given in terms of the model parameters 
as: 


Ei2 — 


«.-A.(2 + a) + 


3M2 


1/2 


(6.83) 


Here A* is the initial value at the time of bounce. 
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While obtaining the above expression we have assumed the following three con¬ 
straint conditions: 


Ho(l + (a-f/3A)(l + 4A) 

<< i. 

(6.84) 

Pc 

{a + /3A)(1 -h 4A) << 1, 

(6.85) 

{4a + /3)^^ ■ 

(6.86) 


Since we are in the small A limit, these conditions are satished if we also choose the 
parameters of the model appropriately. 

Further substituting the expression for A back into the Friedmann equation, we 
get the following expression for scale factor as: 


a(t) = Eisexp 
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where Ei^, is the arbitrary integration constant given in terms of the model parameters 
as: 


Ei3 = aiexp 
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where we introduce two constants E' and E" dehned as: 
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(6.87) 


fi''- ( 1 + 1 ) (i-|)(i+“>)+(!)+ 4 “+'5) (1+1 


( 2 a + I) (1 - 1^(1 + a)) . 


( 6 . 88 ) 


Here a* is the value of the scale factor at the time of bounce. 

ii) Late time: 


At late time, within the interval a' < a < amaxi the result is independent of any 
specific choice of potential. We have already seen that for late times for the hilltop 
potential did not depend on the form of the potential, hence the results shown in 
that section holds good for Coleman-Weinberg potential. 

B. Contraction 
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As has already been mentioned that while performing the analysis for expansion, 
the conclusions for contraction phase is independent of any choice of potential, hence 
the conclusion remains the same as mentioned earlier for other potentials within 
LQG setup. 

C. Expression for work done 

The expression for work done is given by the sum of Eq. (6.42) and Eq. (6.38). 
Since the evolution of the scale factor with time is same as that for the hilltop poten¬ 
tial, the hnal expression for the work done also remains the same as that obtained 
for the hilltop potential. But now the constants depend on the parameters of this 
model under consideration. Hence we can see that the phenomenon of hysteresis is 
possible for Coleman-Weinberg potential in LQG setup. 


6.7 Graphical Analysis 
6.7.1 Case I: Hilltop potential 

All the graphs in this section and in the following sections have been plotted in units 
of Mp =1, Ho = 1, c = 1, where Mp is the Planck mass, Hq is the present value of 
the Hubble parameter and c is the speed of light. 

In Fig. 23 we have shown the evolution of the scale factor and the potential 
for early and late time expansion phase. From the above plots, we can draw the 
following conclusions: 

• Fig. 23(a), shows the plot of the scale factor in the small held limit for hilltop 

potential given by Eq. (6.51) with the parameter values Vq = 1.8xlO“^Mp, p = 
3, (3 = 1, Eo = lOMp, Pc = 0.86Mp^ = 1. 

• Detail graphical analysis show that expansion is possible for both positive and 
negative values of (3. However, for large positive values of (3 (such as > 2) 
expansion is possible only if the values of Vq and Eq are large. For (3 < —1, 
expansion is not possible. And for (3 close to —1, expansion is obtained only if 
Vo and Eq is 0(10“®). Increase in p only results in an increase in the amplitude, 
keeping the pattern same. The plot is almost independent of the value of Pc 

• Fig. 23(b) shows the plot of the behavior of the potential with time for 
small held hilltop potential. This graph has been obtained with the help of 
Eqn. (6.48) with parameter values Vq = 1.8xlO“®Mp, p = 3, (3 = 1, Eq = 
0.38Mp, Pc = 0.86Mp. The evolution of the potential is not much ahected 
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a vs t for LQG model (early time) 


V(0(t)) vs t {^«Mp) for LQG model 




t(inM/) t(inMp-') 


(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during the early expansion phase tial during early expansion phase for (j) << Mp 
for (j) « Mp with Vq = 1.8xlO“^Mp, p = 3, /? = with Vq = 1.8xlO“^Mp, p = 3, /? = 1, Eq = 
1, Eo = lOMp, Pc = 0.86M4, El = 1. 0.38Mp, pc = OMM^ . 


a vs t for LQG model (late time) 


V(M))vstforLQG model 



(c) An illustration of the behavior of the scale (d) An illustration of the behavior of the po- 
factor with time during late time expansion phase tential during late time expansion phase with 

fc = -1.0. Aa = 0.1M2, =-30, yfc = 1.0 . 


with E 2 = Mp 


E'2 = 0 . 6 , 


Figure 23. Graphical representation of the evolntion of the scale factor and the potential 
dnring the expansion phase for LQG model. 
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by the value of p and pc (only the amplitude changes). /3 < 0 do not give the 
required evolution of the potential for causing expansion. Expansion is possible 
for larger positive values of (3 only if Vq is close to Larger the value 

of Vo; more linear the fall of the potential becomes. Expansion is possible only 
if the value of Eq < IMp. 

• Fig. 23(c) shows the plot of the scale factor in the late time expansion phase 
for hilltop potential given by Eqn. (6.54). This plot has been obtained for 
E 2 = Mp, E 2 = 0.6, k = —1.0. Since the solution is true for any potential, 
Fig. 23(c) is also true for all the three potentials. Detail graphical analysis 
have shown that the pattern of the graph remains same for k = 1, only its 
amplitude decreases. Expansion is possible only for E 2 < 1. Larger values 
of E 2 increases the amplitude of expansion. Expansion is possible for both 
negative and positive values of the integration constants. 

• Fig. 23(d) shows the evolution of the potential in the late time expansion phase 
for hilltop potential. Fig. 23(d) has been obtained by with the help of Eqn. 
(6.55) with the parameter values E 2 = O.IM^, E 2 = —30, k = 1.0. Expansion 
is possible for negative values and slightly positive values of E 2 only and for 
getting the correct form of the potential with positive E 2 , we need to make E 2 
much less than unity. Very large values of E 2 do not give the required nature 
of potential. Larger the magnitude of E 2 , larger is the height of the potential. 
In this case, we get a potential decreasing with time, which is required for 
expansion, only for k = 1 and not for k = —1. 

In Fig. 24 we have shown the evolution of the scale factor and the potential for early 
and late time contraction phase. From the above plots, we can draw the following 
conclusions: 

• Fig. 24(c) shows the plot of the scale factor in the early time contraction 

phase for hilltop potential given by Eqn. (6.61). This plot has been obtained 
for E 3 = E 4 = 10^, pc = 0.86Mp. Detail graphical analysis have 

shown that as we increase the value of F4 and F3, the amplitude of expansion 
increases but the nature of the graph remains the same. Both the amplitude 
and nature is almost independent of any change in the value of pc- 

• Fig. 24(d) shows the evolution of the potential in the early time contraction 
phase for hilltop potential. Fig. 24(d) has been obtained by with the help 
of Eqn. (6.59 with the parameter values Vq = 6.3xlO“^Mp, p = 5, f3 = 
0.35, Eg = 10.2Mp, E 3 = 5M“^, pc = 0.86Mp = 1.4. Detail graphical analysis 
show that in order to get the correct nature of the potential which will result 
in contraction, we need /3 > 0. The other parameters do not affect the nature 
of the graph, but only increases its amplitude. 
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(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the po- 
factor with time during the late contraction phase tential during late contraction phase with Vq = 
with A 4 = 100M-\ r = 6 . A 3 = 152M-\ A 4 = 730Mp, P = 

0.52, r = 1.62, p=l . 


a vs t for LQG model (early time) V(^(t)) vs t (^«Mp)for LOG model (early time) 



(c) An illustration of the behavior of the scale fac- (d) An illustration of the behavior of the po- 
tor with time during early time contraction phase tential during early time contraction phase with 
with A 3 = 10-®Mp-i, A 4 = 10^ Pc = 0MM^. Vo = O.SxlQ-^M^, p = 5 , /3 = 0.35, = 

10.2Mp, A3 = 5Mp-i, Pc = 0.86M4 . 

Figure 24. Graphical representation of the evolntion of the scale factor and the potential 
dnring the contraction phase for LQG model. 
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• The solution for the late time contraction phase given by Eqn. (6.64), con¬ 
tains inverse functions. Hence graphically solutions can be obtained only after 
solving the equation numerically. 

6.7.2 Case II: Natural potential 

In Fig. 25 we have shown the evolution of the scale factor and the potential for early 
expansion phase. From the above plots, we can draw the following conclusions: 

• Fig. 25(a), shows the plot of the scale factor in the small field limit for natural 
potential given by Eq. (6.71) with the parameter values Vq = 4.1xlO“^Mp, Eg = 
10 ^ during the early phase of expansion. 

• Higher values of Vq and Eg only increase the amplitude of expansion, keeping 
the nature of the plot unchanged. The change in the amplitude and nature of 
the plot is almost independent of the value of pc- 

• Fig. 25(b) shows the plot of the behavior of the potential with time for 
small field natural potential during the early phase of expansion. This graph 
has been obtained with the help of Eq. (6.69) with parameter values Vq = 
1.3xlO“^Mp, / = 20Mp, Eg = 33, Pc = 0.86Mp. Large values of Vq and / 
increases the height of the potential. Only for certain range of values of Eg, we 
get the required evolution of the potential just like in case of the behavior of 
natural potential in DGP model. The evolution is almost independent of the 
value of Pc- 

• The late time expansion plots remain same as given in Figs. 23(c) and 23(d). 
The amplitudes can be adjusted accordingly by changing the magnitudes of 
different parameters. 

• Fig. 25(c), shows the plot of the scale factor in the large field limit for natural 
potential given by Eq. (6.75) with the parameter values Vq = 1.5xlO“^Mp, / = 
IMp, Eio = 50Mp, Ell = 1, Pc = 0.86Mp during the early phase of expansion. 

• Higher values of the pc, Vg and / only increase the amplitude of expansion, 
keeping the nature of the plot unchanged. 

• Very large values of Vg (close to makes the expansion uneven. Expansion 
is possible for both negative and positive values of Eig. Increasing the value 
of Eig can both increase and decrease the amplitude and makes the expansion 
more non linear and linear alternately. 

• Fig. 25(d) shows the plot of the behavior of the potential with time for large 
field natural potential. This graph has been obtained with the help of Eqn. 
(6.73) with parameter values Vg = 4.4xlO“^Mp, / = IMp, Eig = 6.4Mp, pc = 
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avstforLQG model 


V(^(t)) vs t (^«f) for LOG model 




t(in Mp'^) t(in Mp'^) 

(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during the early expansion phase tial during early expansion phase for (p « f with 
for (/) « / with Vb = 4.1 x10-3m4, Eg = lO^. Vb = l.SxlQ-iM^, / = 20Mp, Ag = 33, = 

0.86M4 . 

a vs t for LOG model (early time) V(0(t)) vs t (^»f)for LOG model 




t(inMp‘^) t(inMp'^) 

(c) An illustration of the behavior of the scale (d) An illustration of the behavior of the poten- 
factor with time during early time expansion tial during early time expansion phase for ^ >> / 
phase for p » f with Vb = 1.5xlO“^Mp, / = with Vq = 4.4xlO“^Mp, / = IMp, Eio = 
IMp, Eio = 50Mp, All = 1, Pc = OMM^. QAMp, p^ = OMM^ . 

Figure 25. Graphical representation of the evolution of the scale factor and the potential 
during the early time expansion phase for the LQG model. 
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0.86Mp. Detail analysis show that larger values of / gives smoother evolution 
of the potential. There are certain ranges of values for Vq and Eiq, appearing 
after certain intervals, which gives rise to the correct nature of potential. Larger 
values of pc may decrease or increase the height of the potential. 

• If we compare the amplitudes of the scale factor in Fig. 24(c) with Fig. 25(a) 
or 25(c), we hnd that after one cycle of expansion and contraction, we can 
get a net increase in amplitude of the scale factor provided the parameters are 
chosen properly. 

6.7.3 Case III: Coleman-Weinberg potential 

a vs t for LQG model (early time) V(^(t)) vs t (0«Mp)for LOG model 




tpnMp-’) t(inMp-') 

(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the po- 
factor with time during the early expansion phase tential during early expansion phase for (f> << 
for (j) « Mp with Vq = 3xlO“^Mp, a = Mp with Vq = 2.7xlO“^Mp, /3 = 0.4, E 12 = 
0.08, /3 = 0.7, E' = 18.4, = 0.86, A 13 = 1. 0.61Mp, a = 0.02, = 0.86 . 

Figure 26. Graphical representation of the evolntion of the scale factor and the potential 
dnring the early expansion phase for the LQG model. 

In Fig. 26 we have shown the evolution of the scale factor and the potential for 
early small field expansion phase for supergravity potential. From the above plots, 
we can draw the following conclusions; 

• Fig. 26(a), shows the plot of the scale factor in the small field limit for su¬ 
pergravity potential given by Eqn. (6.87) with the parameter values Vq = 
3xlO“^Mp, a = 0.08, = 0.7, E' = 18.4, pc = 0.86, E 13 = 1 during the early 

phase of expansion. 
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• Expansion is possible for both negative and positive valnes of f3. However, for 
negative valnes of f3, expansion is possible only if a > 0.5 and E' > 1. Larger 
valnes of Vq and E' increases the amplitnde of expansion. The expansion is 
almost independent of the valne of pc- 

• Fig. 26(b) shows the plot of the behavior of the potential with time for small 
held snpergravity potential dnring the early phase of expansion. This graph 
has been obtained with the help of Eqn. (6.82) with parameter valnes Vq = 
2.7xlO“^Mp, (5 = 0.4, Ei 2 = 0.61Mp, a = 0.02, pc = 0.86. Correct form of 
the potential is obtained for both negative and positive valnes of (3. However, 
for negative valnes of /3, E 12 mnst be mnch less than nnity. Larger the valne of 
Vo; smaller the valne of E 12 allowed. For large positive /?, proper evolntion of 
the potential is possible only for very small valne of Vq. pc has very negligible 
effect on the evolntion of the potential. 

• If we compare the amplitndes of the scale factor in Fig. 26(a) with Fig. 24(c), 
we hnd that after one cycle of expansion and contraction, we can get a net 
increase in amplitnde of the scale factor provided the parameters are chosen 
properly. 


7 Hysteresis from Einstein-Gauss-Bonnet brane world gravity 
model 


Einstein Ganss-Bonnet gravity is based on the theory of higher (Eve) dimensional 
brane world scenario where the simplest RS brane world action for single brane 
model (RSH) is modified by the addition of a Ganss Bonnet term {Cgb) which is 
a combination of the Ricci scalar R, the Ricci tensor Rab and the Riemann tensor 
Rabcd and is given by: 

Cgb = R" - ^RabR^^ + RabcdR^^^^- (7.1) 


Gonseqnently the action in five dimensional space-time takes the following form 

I H'-‘ X ./ — 7i R — \ rv { — A R . ^ R^^ A- R . ^ 




R-2A + a[R^ - ARabR^^ + RabcdR^^^^ ' 

+ j {Ctt"' - a) , 


(7.2) 


where qab is the metric in the 5D bnlk and 


g,, = d,X^d,X^~gAB (7.3) 

®One can generalize this action in arbitrary dimension D > 5. It is important to note that at 
D = 4 the contribution from Gauss-Bonnet term is a topological invariant and does not contribute 
to the field equations. For rest of the analysis we fix the space-time dimension to be Z? = 5 for the 
sake of simplicity. 
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is the induced metric on the brane with X^[x^) being the coordinates of an event on 
the brane labeled by 

In the present context the modihed Friedmann equation is given by [78]; 


K 




h{a) 


+ eH‘ 


4q! f 3k — ehia) ^ ^ 
1 + — - 


-[ 2 


(7.4) 


where a is the single brane tension, a is the Gauss-Bonnet coupling, A is the 5- 
D cosmological constant, e = -fl, — 1 for space-like or time-like extra dimension 
respectively. Here the the function h{a) is given by the following expression: 


h{a) = ek + 






(7.5) 


where /r is a constant. 

Since the above equations are highly complicated, following the analysis of [78], 
we can write the above equation as: 

C{p + aY={A± H^) {B + H^)\ (7.6) 


where ± corresponds to e: = +1 and e = —1 respectively. Here A,B and C are 
dehned as: 


A : 

k h{a) 

“ a2 ^ a2 ^ 

) 



(7.7) 

B : 

3k 3 

2 a2 8a 

eh{a) 

3 3k 

eA 

(7.8) 

2 a2 

2 ’ 

C : 

36 \Sa) 

> 0. 



(7.9) 


7.1 Condition for bounce 
A. Space-like extra dimension with e = 1 


Let p = pfe is the corresponding energy density at which bounce occurs i.e. the 
Hubble parameter H = 0 condition achieved. Hence substituting the Hubble param¬ 
eter if = 0 in Eq. (7.6), we get: 


Pb = 


^/AB 

~C^ 


a 


(7,10) 


where C = y/C is also a constant. A, B are now hxed at the bounce. 
Total mass at the bounce is given by 


'^/AB 


Mb = Pbol 


C 


— a 


(7.11) 











and infinitisimal change in mass at bounce can be expressed as: 


6 Mk 


C 


fB5Aal 
V 2A 


+ 6Bal + 3Bal6ab 


(7.12) 


where inhnitisimal change in A, B i.e. <5^4, 6 B can be expressed using Eq. (7.9) with 
£ = +1 as: 


6 A 

SB 


—2k ^ 

- - — Od “ 1 “ -- 7 - 

c? (1 - (A - A) 4a) 


2 k yua ® 

- -da — - -- ITT - -da. 

“ (1 - (-4 - Ji) 4a) 


(7.13) 

(7.14) 


Further substituting Eq.(7.13) and Eq. (7.14) into Eq. (7.12) and using the energy 
conservation equation (or continuity equation) we get the following expression for 
the change of the amplitude of the scale factor after one cycle as: 


Sa . =_ 

^^mzn /o 9 


X') 


(7.15) 


where we introduce a new model dependent parameter X' dehned as: 


X' = 


7)4 

C" 


-2k + 


-2 



3Bat 


. (7.16) 


Thus we see that the condition for an increase in the amplitude of the scale factor 
depends on A through X', which in turn depends on the curvature parameter k and 
the model parameters like /i, a etc. 

For k = +1, 0, —1, Eq. (7.15) can be rewritten as: 


f pdV 

(3aai,„ - V') 


for k = +1 


damin ^ 


f pdV 




for k = 0 


(7.17) 


f pdV 


for k = —1. 
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where J' and Z' defined as: 


X' = 


V — ZA 
^ c 



^- 2 + 




1- A-- 


^— 
'min 


Aa 


W' = # 


-2 \ 


-2 \ 


\ _B_/ \ _|_ Q 

'l-4Aa) J 2A Wl-4Aa) J 


mm 


y! _ ZA 


2 + 




1-1 A+^^ )4a 


min 



-2 + 




+ 3Ba%, 


, for k 


for k 


|4o 

min 


(7,18) 

and now the expressions for A and B for different values of curvature parameter 
k = +1, 0, —1 are given by: 


, for k 


and 


/i(a^ 


A = M ^min) 


h{ar, 


(1 

\ ^min 


for k = +1 


for /c = 0 


, for k = —1. 


(7.19) 


2 ai 
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h{ar^ 


2 ai 


3 3 

ScK 2 a'L. 


^ 3 /i(ci^ 

So; 2al 


3 
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A 

J’ 


2 a 


3 h(ar 
So; 2a^ 


So; 2al 


A 

-, for k = +1 

2 


for k = 0 


A 

—, for k = —1. 
2 ’ 


(7.20) 


For the various values of the curvature parameter k = +1,0, —1 one can point out 
the following characteristics from our analysis: 
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= +1 


= 0 


= - 1 . 



• For k = +1 case if we consider the situation where 3cra^j„ > F', i.e. if the 
brane tension to be large, then an increase in the amplitude of the scale factor 
is possible for a positive sign of the work done i.e. § pdV > 0. And if the brane 
tension is small enough then a negative sign of the integral gives expansion 
with increasing amplitude. 


• If we follow the analysis of [78], it shows pictorially that bounce does not occur 
when the extra dimension is space-like along with fc = 0. We can also do similar 
pictorial study for the case when k ^ 0 and see whether this statement remains 
true for /c 7 ^ 0. 


• We also see that for k = —1 case depending on the relative signs of A, B and 
O', Samin can be positive for both § pdV > or < 0 cases. Thus just like for 
k = -|-1, if 3craF„ > Z\ i.e. if we consider the brane tension to be large, then 
an increase in the amplitude of the scale factor is possible for a positive sign 
of the work done i.e. f pdV > 0. Additionally if the brane tension is small 
enough then a negative sign of the integral i.e. f pdV < 0 gives expansion with 
increasing amplitude. 


B. Time-like extra dimension with s = —1 


In this context the condition for bounce is given by Eq. (7.10) with the expres¬ 
sion of A and B now given by Eq. (7.9) with e = —1. Hence the expression for 
inhnitesimal change in mass at bounce dMf, is again given by Eq. (7.12), where the 
expressions for 6 A and 6 B are now given by: 


dB 


2 k ^ 

H—5-0C1 


2 k ^ 
■^oa 


pa 


-5 


(l+(A+i,)ia) 


pa 


-5 


(l + (A+i)4a) 


5a, 


5a. 


(7.21) 

(7.22) 


Further substituting Eq. (7.21) and Eq. (7.22) into Eq. (7.12) and using the energy 
conservation equation (or continuity equation) we get the following expression for 
the change of the amplitude of the scale factor after one cycle as: 


5a 


min 


f pdV 


(7.23) 


133 





where the parameter X' is defined as: 


X' = 


C 


2 k + 


-2 


1 + (A+^)4a 

min 


2 k + 


B 

Ia 


-2 

k'^min 




+ 3Baz 


. (7.24) 


For k = +1, 0, —1, Eq. (7.23) can be rewritten as: 


^ Ojrn.in. ^ 


f pdV 


(3<^o:Ln - i") 


§ pdV 


- J') 


^ pdV 




for k = +1 


for A: = 0 


for k = —1. 


(7,26) 


where Y', J' and Z' defined as: 


V' — YA 
^ c 


2 + 




X' = K j' — Ya 

J — C' 
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(l+4Aa) J 2A Wl+4Aa) 


+ 3-SaLn 


7' — YA 

Zj — C! 


-2 + 
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B_ 

2A 
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+ 3Bai 


+ ^BaY 

(7.26) 

and now the expressions for A and B for differfent values of curvature parameter 
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k = +1, 0, —1 are given by: 


4a 


h{a) 


A = 


,2 ’ 


1 

4a 


' —1 + + -aA — 1 I , for /c = +1 

^min 4 


for k = 0 


'-l + 4^ + |aA-l 


for k = —1. 


(7.27) 
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for k = +1 


B = 


3 kictYriin^ 3 A. 

— ^ — + _ for A: = 0 

8 a 2 aL„ 8 a 2 ’ 


A 


for k = —1. 


(7.28) 


8 « ‘^aLn 2 ’ 

For the various values of the curvature parameter k = +1,0, —1 one can point out 
the following characteristics from our analysis: 


• Here for k = +1 case we see that depending on the relative signs of A, B and 

O', Samin can be positive for both f pdV > or < 0. Thus, if > Y', i.e. if 

we consider the brane tension to be large, then an increase in the amplitude of 
the scale factor is possible for a positive sign of the work done i.e. f pdV > 0 . 
Additionally if the brane tension is small enough then a negative sign of the 
integral i.e. f pdV < 0 gives expansion with increasing amplitude. 

• In [78], the authors show pictorially that the bouncing condition is reached if 
either 0 < H < 2A or B > 2A and a < 0 for k = 0 case. We can also do 
similar pictorial study for the case when k ^ 0 and see whether this statement 
remains true for k 0. 

• Also for k = —1 case we see that depending on the relative signs of A, B and 

O', damin can be positive for both §pdV > or < 0. Thus just like for k=+l, if 

3(Ta^j„ > Z', i.e if we consider the brane tension to be large, then an increase 
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in the amplitude of the scale factor is possible for a positive sign of the work 
done i.e. ^ pdV > 0. And if the brane tension is small enough then a negative 
sign of the integral i.e. ^ pdV < 0 gives expansion with increasing amplitude. 

7.2 Condition for acceleration 
A. Space-like extra dimension with e = 1 


In order to get the condition for acceleration, we need to use the second Fried¬ 
mann equation which we get after differentiating Eq. (7.6) w.r.t. time t for e = -fl 
and using the energy conservation equation to get rid of p. The resulting equation 
is given by: 


H + 


d 

a 


2 >C'‘^{p + p){p + a) ^ 


Z 

Y 


(7.29) 


where Z and Y are functions of A, B and their time derivatives, which are in turn 
functions of scale factor a{t) and its time derivative a respectively. Their explicit 
expressions in terms of A, B and their derivatives are given by: 


Y = + 4:AB + 2B^ + 6H^ + 8BH^, (7.30) 

AB^ AB B 

Z = —— + 2—— + Ah^ + 2 H{Ab + AB + BB + BHA- 
H H 

(7.31) 


where H is the Hubble parameter. 

The general condition for acceleration is given by the following expression: 

(p + p)(p + ,j)<A^(H^Y-Z) (7.32) 

Using Eq. (7.10) and setting = 0 in the above equation, we get the following 
condition for acceleration at bounce as: 

VAB ( Z ^ aC \ 

^TIb) 

where A,B, Z are fixed at bounce. Thus, whether this condition violates the energy 
condition or not, depends upon the values of the different parameters of the model. 

Finally, in terms of the scalar field and its potential, the condition for acceleration 
can be written as: 



^ -Z)-A\V{^) + a) (7.34) 

Thus we see that the condition for acceleration in terms of the scalar held now not 
only depends on the kinetic term, but also on the fourth power of the time derivative 
of 0. 
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B. Time-like extra dimension with e = —1 


The second Friedmann equation is again given by Eq. (7.29) i.e. 


H + H^ 


a 

a 


+ p){p + a) faV 
D \a) 


E 

D' 


where now the expression for D and E are given by; 


D = AAH^ + AAB - 2B^ - 6H^ - 8BH^, 

AB^ ABB 

E = —— + 2 —— + AH^ + 2 H{Ab + AB-BB- BH^. 
H H 


(7.35) 


(7.36) 

(7.37) 


Rest all the conditions for acceleration in terms of the pressure and the scalar held re¬ 
mains the same with expressions for D and E now given by Eq. (7.36) and Eq. (7.37) 
respectively. 


In Figs. 27, we have shown the phenomena of bounce and acceleration in the 
Einstein Gauss Bonnet gravity model. While performing this analysis, we have con¬ 
sidered that p is very small and hence can be neglected. We can draw the following 
conclusions from the above hgures: 

• Figs. 27(a) and 27(b), have been plotted with the help of Eqns. (7.4), (7.5), 

(7.6) and (7.9) using the relation p = with k = 0, A = lOM^, B = 

40Mp, cr = — C = 100. Bounce occurs for both negative and positive 

brane tension. Detail analysis of the scenarios under which bouncing is possible 
has been discussed extensively in [78]. 

• From Fig. 27(a), we get the bounce at p = pb = 12.78Mp. 

• In Fig. 27(b), we have shown the bouncing condition for k = 1, £ = 1, w = 
1/3, C = 100, A = 10M2, B = 40M2, a = -lO-^M^, C = 100, A = 
6 xlO“^Mp. The condition remains valid for the other combinations like (k = 
—1, £ = —1), (k = 1, e = —1) and (k = —1, £ = 1). But for the last case, 
bouncing occurs at a very small value of p, while for the other cases, it occurs 
at some larger value ofp, which depends upon the choice of our parameters. 

• For the cases (k = 1, s = ±1), bounce is possible only for negative brane 
tension. But for (k = —1, £ = ±1), it is possible for both positive and negative 
brane tension. Bounce is also possible for ta = 0. From Fig. 27(b), we get the 
bounce at p = pb = 1.9Mp. 
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H vs p for Einstein Gauss 



10 15 20 25 30 35 40 


P(inMp^) 


H vs p for Einstein Gauss 
Bonnet gravity model (k=1 ,e=1) 



(a) An illustration of the bouncing condition for a (b) An illustration of the bouncing condition 
universe withfc = 0, A = lOMp, B = A0Mp, cr = for a universe withfc = 1, e = 1, w = 
C= 100. 1/3, C = 100, A = 10M2, B = dOM^, a = 

-lO-^Ml, C = 100, A = 6x10-3m2. 


a vs Ph for Eisntein Gauss 
Bonnet gravity model (k=0) 



(c) An illustration of the acceleration condition at the time of bounce for 
a universe with an equation of state w = 1/3, k = —0, A = —lOM^, B = 
10M2, C' = 10, fj = 1M4 


Figure 27. Graphical representation of the phenomena of bounce and acceleration for 
Einstein Gauss Bonnet gravity model. 


138 














• Fig. 27(c), we have shown the necessary condition of acceleration (a > 0) at the 
time of bounce for w = 1/3, /c = 0, A = —lOM^, B = lOM^, C = 10, a = 
IMp. This plot has been obtained with the help of Eqn. (7.29). But since Eqn. 
(7.29) is highly complicated, in order to know the parameter space for which 
we get acceleration at bounce, we have instead set = 0 and p = in Eqn. 
(7.29) and analysed for what values of the parameters do we get acceleration 
at bounce. From detail graphical analysis, we have found that this possible if 
the condition AAB + 2B‘^ < 0 is satisfied for {p + a)> 0, and vice versa. This is 
possible for both positive and negative values of A, B, a provided the above 
conditions are satisfied. Similar analysis can be done for k = 1,-1 also. The 
results remains nearly same. 

• Thus from Fig. 27 we can conclude that for this model, bounce is possible for 
closed, open and flat universe. 


7.3 Condition for turnaround 
A. Space-like extra dimension with e = 1 


Following the same line of treatment as for bounce we get the condition for turnaround 
for £ = +1 as: 


Pt 


^/AB 

C 


— a 


(7.38) 


where A and B are fixed at turnaround. 

Following the same analysis as for bounce we get the following expression for 
change in the amplitude of the scale factor after each successive cycle is given by: 
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(7.39) 
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for k = —1. 


where where the expressions for E, J, Z remains same as in the bounce section 
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for y', J', Z', only with ttmin replaced by ttmax he. 
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(7.40) 

and now the expressions for A and B for differfent values of curvature parameter 
k = + 1 , 0 , —1 are given by: 
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(7.42) 


Therefore the conclusions also remain same as for bounce case. 
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In [78], the authors show pictorially that for the case when fc = 0, re-collapse 
occurs when A > 0, i? > 0 and a < 0 or is small but positive. We can also do similar 
pictorial study for the case when k ^ Q and hnd similar conditions on A, B and a 
which gives re-collapse. 

B. Time-like extra dimension with e = —1 


The condition for turnaround remains same as Eq. (7.38). 

Following the same analysis as for bounce we get the expression for change in 
the amplitude of the scale factor after each successive cycle as: 
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for k = —1. 


where the expressions for Y, J, Z remains same as in the bounce section for 
only with amin replaced by amax i-e. 


V — ZA 
^ c 


2 + 


—2 


\ \ ^max / 


^-|2 + 


—2 
^^max 


\ \ ^max J / 


X = I J — Va 
•J — C' 


7 _ YA 

Zj — C! 


f-^^max \ B 

(l+4Aa) J 2A 




{l+AAa) 


2 

max 


-2 + 


-2 


1+1 

^max 


4o 


2A 


-2 + 


-2 


1+1 -4-V— 14a 


-|- 3Bal 


(7.44) 


for k 


for k 


, for k 


and now the expressions for A and B for different values of curvature parameter 


- 141 - 


= +1 


= 0 


= - 1 . 


























k = + 1 , 0 , —1 are given by: 
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Therefore the conclusions also remain same as for bounce case. 

In [78], the authors have shown pictorially that re-collapse occurs if 0 < i? < 2A 
for = 0. We can also do similar pictorial study for the case when 7 ^ 0 and hud 
similar conditions on A, B and a which gives re-collapse. 
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7.4 Condition for deceleration 
A. Space-like extra dimension with e = 1 


In this subsection the analysis again remains the same giving in place of Eq. (7.32), 
the governing constraint condition is given by: 


(p + p)(p + ^)>^(H'^y-z), 

and in place of Eq. (7.33) and Eq. (7.34), we get: 
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Thus we again see that this condition now depends on the fourth power of the time 
derivative of the scalar held as mentioned earlier. 

B. Time-like extra dimension with e = —1 


In the presesnt context the analysis again remains the same as the hnal govern¬ 
ing equations are exactly same as Eq. (7.47), Eq. (7.48) and Eq. (7.49) in which only 
the signatures of Y, Z, A, B changes accordingly due to e = — 1. 

In Fig. 28, we have shown the phenomena of turnaround and deceleration in 
Einstein Gauss Bonnet gravity model. We can draw the following conclusions from 
the above hgures: 

• Fig. 28(a), has been plotted with the help of Eqns. (7.4), (7.5), (7.6) and 

(7.9) using the relation p = with fc = 0, A = 1, i? = 10, a = 

— C = 200. Detail graphical analysis show that the condition of 
turnaround is achieved only when 2l > 0, i? > 0, a < 0 or is small positive. 
For further detail analysis of this case, one may refer to [78]. 

• From Fig. 28(a), we get the turnaround at p = pt = 0.7Mp. 

• Fig. 28(b) shows turnaround using Eqns. (Eqns. (7.4), (7.5), (7.6) and (7.9)) 

with w = 1, k = 1, e = 1, a = C = 200. For this case we have used 

w = 1, since we require a stiff equation of state for causing contraction and 
deceleration. The parameter space causing turnaround remains same as that 
for bounce case. From Fig. 28(b), we get the turnaround at p = pt = l.OM^. 

• Fig. 28(c), we have shown the necessary condition of deceleration (a < 0) at 
the time of turnaround for w = 1, fc = 0, A = 10, i? = 15, a = liW^, C = 10. 
This plot has been obtained with the help of Eqn. (7.29). But since Eqn. (7.29) 
is highly complicated, in order to know the parameter space for which we get 
deceleration at turnaround, we have instead set = 0 and p = ph fa Eqn. 
(7.29) and analysed for what values of the parameters do we get deceleration 
at turnaround (similar to the analysis done for bounce). From detail graphical 
analysis, we have found that this possible if the condition AAB + 2B^ > 0 is 
satished for {p + (j)> 0, and vice versa. This is possible for both positive and 
negative values of A, B, a provided the above conditions are satished. Similar 
analysis can be done for k = 1,-1 also. The results remains nearly same. 

• Thus from Fig. 28 we can conclude that for this model, turnaround is possible 
for closed, open and hat universe. 
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H vs p for Einstein Gauss H vs p for Einstein Gauss 

Bonnet Gravity modei (k=0, £=1) Bonnet gravity modei (k=1 ,£=1) 




(a) An illustration of the turnaround condition (b) An illustration of the turnaround condition 
for a universe with k = 0, A = IM^, B = for a universe with an equation of state w = 
lOM^, cr = C = 200. 1, fc = 1, e = 1, cr = -lO-^M^, C = 200. 


a vs pt for Eisntein Gauss 
Bonnet gravity model (k=0) 



(c) An illustration of the deceleration condition at turnaround for a uni¬ 
verse with an equation of state w = 1, k = 0, A = 10, B = 15, a = 
lO-^M^, C" = 10. 


Figure 28. Graphical representation of the phenomena of turnaround and deceleration 
for Einstein Gauss Bonnet gravity model. 
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7.5 C. Evaluation of work done for one cycle 

The expression for the total work done is same as that given by Eq. (4.52) in DGP 
model. But we can also express it in terms of the scale factor by using Eq. (7.29) 
and Eq. (7.6) in order to get an expression of the pressure in terms of the scale 
factor, which we can then substitute into the integral for the work done to get: 

A. Space-like extra dimension with e = 1 


fpiv = f3a^l(A + + H^) 

+ j3i^(A + HyBfB + it. 

(7.50) 


B. Time-like extra dimension with e = —1 

jp,v = a H-) (^± 

+ j,3h^,{A - + H^-) it. 

(7,61) 

Thus the work done now depends not only on the scale factor, but also on the dif¬ 
ferent parameters of the model like the coupling constant, brane tension etc within 
the present setup. 

7.6 RSII limit 

It has been shown in [79], that in the limit, a —)■ 0, Eq. (7.4) reduces to RSII scenario. 
Then all the results shown in [11], for the case when all the other parameters are 
made zero and the extra dimension is time-like, will also hold true for the above 
model. These results have been shown once again in the appendix for the sake of 
clarity. 

7.7 Semi-analytical analysis for cosmological potentials 

In the analysis below, we will denote Planck mass by Mp. Here we will perform the 
analysis for fc = 0 case. One can repeat the analysis for k = ±1 as well. 


=F 1 dt 


=F 1 df 
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7.7.1 Case I: Hilltop potential 
A. Expansion 


Al. Space-like extradimension {e = +1) 

Using Eq. (4.67) and Eq. (7.6) for e = +1 and using the explicit form of p, which is 
specified by the hilltop potential as mentioned earlier, we get an integral equation of 
the following form; 




-{^^+1:^)'2Va + VaMr±2 


jVo {l + p (jf^) 


v^/?p(k 


p-1 


Un 


3M2 


dt, 


where 


R = A 


[C 

— -3B±2i -a 
AA \ A 


(7.52) 

(7.63) 


While arriving at the above integral, we have considered that p is small enough such 
that the term ap/a^ can be neglected and a should be small such that A is a large 
number hence the following constraint holds good; 




(7.54) 


By neglecting p, the quantities A and B becomes independent of scale factor, hence 
allowing us to attain some analytical approximate solutions. To compute the left 
hand side of the above integral equation, we will again use the following redefinition 
of the field; 

^ (7.55) 




= e 


Analytical solution was possible only for the case in the small filed limiting case i.e. 
when (p/Mp << 1 which has been discussed below. 

In this limit Expanding the exponentials upto linear order and applying two 
conditions; 










3B±4j^a + 4J^Vo 


(A+f) 


«1, 

(7.56) 

«1, 

(7.57) 
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we get the solution for A as: 


A 


Vopl3V2 

3M2 


t + Fo 



(7.58) 


where we introduce a new constant Fq, which is dehned in terms of the model pa¬ 
rameters as: 

^0 = R + - (^ + (7.59) 

where S is given by: 


5 





(7.60) 


Here A* is the value at bounce. 

Substituting the above expression in the Friedmann equation, we get the expres¬ 
sion for scale factor as 
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where the overall factor Fi can be expressed in terms of the model parameters as: 


Fi = a^exp 


R + 2-/CV,, 


-2(A + - 






p(5-l)Fo 


± ^\Jr + 2VCVq 




2(^+f) 


1/2 


\p{S - 1) 


Vbp/3\/2 

3iW2 


\ 


{R+^^o-{A+ms 


(7.62) 


Here a* is the value of the scale factor at the time of bounce. Thus we see that the 
scale factor varies exponentially with time in the small held limit during expansion. 


A2. Time-like extra dimension 


The results remain same as for space-like extra dimension with the expression for R 
now given by: 


R 




4A 


+ 3B±2 



(7.63) 


B. Contraction 


Following the same procedure as for the expansion phase, but with the expression 
for density is now given by 0^/2, we get the following solutions: 


1 + 


1 - I (-2 (H + f) + VAr) (-t + F 2 ) + 1 


vT?(-2(A+f )+FAR) 


vT?(-2(A+f )+FAR) 


0(t) = 


2tVR 


-2 f H + Varj - F 2 (^-2 +11 + Tar 1 1 


-2 (A + f) + v^) 


+ ^3, 


,(7.64) 


(7.65) 
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where we introduce two arbitrary integration constants F 2 and F 3 , which is given in 
terms of model parameters as: 


F2 

Fs 


^f-2y/R(^-2{A+^) + y/ARj tf{y/AR-B-2A) 


2A + B - ^/AR 


2tf\fR 


A + ^ \ + \fAR ) + F 2 


+ 


2A + B - y/AR ’ 

■2 ^ + \/AR^ tf 

-2 (A + I) + ^/ARj t) 


(7.66) 


(7.67) 


Further substituting back into the Friedmann equation, we get the following solution 
for the scale factor as: 


a{t) = F 4 exp 


[D^ f (-2 + D' + D'F^D" 

V ^ V ?>D'‘^D" 




(7.68) 


where we introduce a arbitrary integration constant F 4 , which is given in terms of 
model parameters as: 


F^ 


ttf exp 



(-2 + D' + D'F^D" 
^D'^D" 


D'D"tf) 



Here D' and D" is given by: 


D' 


D" 


1 

—2(y4 + —) + 'Jar, 


(7.69) 


(7.70) 

(7.71) 


and ttf is the value of the scale factor at turnaround. Here also the dependence is 
exponential but includes higher powers of time. 


C. Expression for work done 


The expression for work done is obtained by substituting the expressions for scale 
factor into the Eq. (7.50), taking care of the assumptions that we have made earlier 
in the context of Einstein-Gauss-Bonnet brane world model. 


Al. Space-like extra dimension: 


149 






















pdV = fsi-fit 

max + f2t max - fst max + fit max - fht max 


fat 

max — ht max + fit min - f2t min + ht min 

fitJnin + htmin “ /e^Ln + A^mm) 

h{Erf[i-fio + fntmax)] - Erf [{-ho + futmin)]) 

(7.72) 


where fi....fii are the constants whose functional form depends on the parameters 
as appearing in the expression for the scale factor. Their explicit forms have been 
given in the appendix. 

B. Time like extra dimensions: 

Since the solutions for the scale factor in the expansion and contraction phase remain 
same baring the constants, the hnal results for work done in this case also remain 
same. 

Thus we see that since the work done after a complete cycle is non-zero, at least 
for the approximated expressions for the scale factor which we have evaluated in this 
paper. Thus hysteresis phenomenon can explained for hilltop potential in the small 
held limit for Einstein Gauss Bonnet brane world gravity. 


7.7.2 Case II: Natural potential 
A. Expansion 

Repeating the analysis as we have done for the hilltop potential, we get the ex¬ 
pressions for the scalar held and scale factor for the two limiting cases-i) small held 
limit (p/f << 1 and ii) large held limit (p/f » 1, which we discuss elaborately in 
the following subsections. 


i) f/f « 1 

Al. Space-like extra dimension (e = -|-1) 

Taking the small argument approximations of the trigonometric functions, we get 
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the solutions for the scale factor and scalar held as 


a(t) = F 5 exp 




1 / 2 ^ 


1/2 


Ht) 

f 


= Fq exp 


Vot 


Ml 

V2 


-2 {A+^)+Vu(l±^^y^^^ 


1/2 


, (7.73) 


(7.74) 


where we introduce three arbitrary integration constants F 5 , Fq and U, which are 
given in terms of model parameters as: 


F 5 


Fe 

U 


tti exp 




1 / 2 ^ 


1/2 


/ 


exp 


VoU 


l/2\ 


^(^-2(7l+f) + v^(l±^) j 


yii? ± 2yCK), 


(7.75) 


(7.76) 

(7.77) 


with the dehnition of R is same as appearing in the previous case i.e. for hilltop 
potential. 

Thus the scale factor depends on time exponentially with only hrst power or 
linear power of time. 


A2. Time-like extra dimension (e = —1): 

The conclusion remain same as was for the hilltop potential. 


ii) 4>/f » 1 

Al. Space-like extra dimension {e = +1): 

Using the large argument approximations of the trigonometric functions and ap¬ 
plying the conditions: 



W » 



(7.78) 

(7.79) 
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i.e. W is a large quantity where it is defined as: 


W = 


1 

2 


Vu - 



(7.80) 


and U has already been defined before, we hence get the solution of the scalar field 
as: 


4>{t) 

f 


2 tan ^ 


exp 




(7.81) 


where we introduce three arbitrary integration constant Fy, which is given in terms 
of model parameters as: 


F7 = f\n 


W 

— tan 



Voti 

3/2 


(7.82) 


Substituting back into the approximate Friedmann equation we get the following 
expression for the scale factor as: 


a(t) = Fg exp 




V 


wVgln(l + e^) 4 fC 


2Vo 
3/2 W 




where the functions G and Fg are given by: 
2 F 7 2Vot 


G = 


W 3PW’ 


Fg = tti exp 


-vw[(i-±j§]u 




+ 4 Ic 


2Vo 
3/2 TV 


+ IF V 


(7.83) 


(7.84) 

.(7.85) 


A2. Time-like extra dimension {e = —1): 

The conclusion remain same as was for the previous case. 

B. Contraction 


Al. Space-like extra dimension {e = +1): 

Since the contraction phase is independent of any potential, the conclusions remain 
same as for hilltop potential 
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A2. Time-like extra dimensions {e = —1): 


In this case also the conclusions remain same as for hilltop potential 

C. Expression for work done 

The expression for work done is obtained by substituting the expressions for scale 
factor into the Eq. (7.50), taking care of the assumptions that we have made earlier 
in the present context. 

Al. Space-like extra dimension {e = +1): 



mm 


f^^min T fd^min T f7tmin 

min — max ^ 


(7.86) 


where / 1..../13 are the constants whose form depends on the parameters in the ex¬ 
pression for the scale factor. Explicit forms of the constants are given in the appendix. 

A2. Time-like extra dimensions (e = —1): 

Since the solutions for the scale factor in the expansion and contraction phase remains 
same baring the constants, the results for work done also remains same. 

Thus we see that since the work done after a complete cycle is non-zero, at least 
for the approximate expressions for the scale factor which we have evaluated in this 
subsection earlier. Thus hysteresis phenomenon is possible for the natural potential 
in the small field limit for Einstein Gauss Bonnet gravity. 

7.7.3 Case III: Coleman-Weinberg potential 
A. Expansion 

Al. Space-like extra dimension (e = +1): 

The original integral equation being complicated, we follow the same procedure as 
has been done before i.e. use the field redefinition; 


0 /Mp = e\ 


(7.87) 
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Analytical solutions were possible only for small field limit 0/Mp << 1 which we 
have discussed below. 

In this limit expanding the exponentials upto linear order and applying the 
following constraint conditions: 


2v^K)(4a +/3)A 
{y/U+ 2y/CV^a) 
2y/CVo{Aa + l3) 
v^(yA^-2(A + f)) 
4/3A 

(4q; + /3) 


«1, 

(7.88) 

«1, 

(7.89) 

«1, 

(7.90) 


we get the solution for redehned or transformed held A as: 


A 


(4q; + /3) 


( v^t 
\ 3M2 



yA^-2 (A + f) 


1/2 




(7.91) 


where we introduce two new constants H' and Fg dehned in terms of the model 
parameters as: 


H' = VU + 2^/CVoa, 

F - 
-cg — 






{Aa + (A) 


1/2 3M2' 


(7.92) 

(7.93) 


Here A* is the value at bounce. 

The solution for the scale factor turns out to be: 


a{t) = Fio exp 


- 2 ^ 


1/2 


+ F(^I{4:Cx + /3) > t 


(4a + /5)t^3^1 


(7.94) 

where we introduce two new constants I and Fio dehned in terms of the model 
parameters as: 


I = 


Fio = tti exp 


V2Vo{Aa + /3) 


-2{A+ §) 


( VaIt 




2{ A+ — 

' 2 


1/2 


+ Fgl(4a /3) f U 


(7.95) 

(4a + /3)t2j^ 


(7.96) 
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A2. Time-like extra dimension (e = —1): 


The conclusions remain same as for e = 1 case, only the expression for different 
parameters change with e = 1 replaced by £ = — 1. 

B. Contraction 

Al. Space-like extra dimension {e = +1): 

Conclusions for hilltop potential, holds true for this case also. 

A2. Time-like extra dimension (e = —1): 

Conclusions remain same as for the case of hilltop potential. 

C. Expression for work done 

Al. Space-like extra dimension (e = +1): 

The general expression for the scale factor i.e its time dependence is same for super¬ 
gravity potential as for hilltop potential in the small held limit, hence the conclusions 
for work done remains valid for this case also, thus giving rise to the phenomenon of 
hysteresis. 

A2. Time-like extra dimension {e = —1): 

The conclusion is same as above. 


7.8 Graphical Analysis 
7.8.1 Case I: Hilltop potential 

All the graphs in this section and in the following sections have been plotted in units 
of Mp = 1, Ho = 1, c = 1, where Mp is the Planck mass, Hq is the present value of 
the Hubble parameter and c is the speed of light. 

Fig. 29 shows the evolution of scale factor and potential during expansion and 
contraction phase for Einstein Gauss Bonnet gravity model. We can draw the fol¬ 
lowing conclusions: 

• Fig. 29(a), shows the plot of the scale factor in the small held limit for hilltop 
potential given by Eqn. (7.61) with Vq = p = 3, Fi = 1, Fq = 
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a vs t for Einstein Gauss 
Bonnet gravity modei 


V(^(t)) vs t (^«Mp)for Einstein Gauss 
Bonnet gravity modei 



(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during expansion phase for << tial during expansion phase for (j) « Mp with 
Mp with Vo = p = 3, Fi = 1, Fo = Vo = 2.7xlO-^M^, p = 4, Fq = 54Mp, A = 

O.lMp, A = 0.1M2, F = 0.1M2, C = 10, S = O.IM^, B = 0AM^,C = 500, i? = 30M^, S = 
0.1, i? = 0.1M2. 0.3, /3 = 0.05. 


a VS t for Einstein Gauss 
Bonnet gravity modei 



(c) An illustration of the behavior of the scale 
factor with time during contraction phase with 
A = lOM^, B = 20M2, Fa = 150Mp, F 4 = 1. 


V(^(t))vstfor Einstein Gauss 
Bonnet gravity modei 



(d) An illustration of the behavior of the po¬ 
tential during contraction phase with Vg = 
10-6A/4, p=l, F2 = 24Mp \ Fa = IMp, A = 
2M^, B = -2M^, R = 10-3M2, /3 = -0.9 . 


Figure 29. Graphical representation of the evolution of the scale factor and the potential 
during the expansion and contraction phase for Einstein Gauss Bonnet gravity model. 
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O.lMp, A = 0.1M2, B = 0.1M2, C = 10, ^ = 0.1, R = O.IM^. 

• Detail graphical analysis show that the amplitude of expansion increases with 
increase in the value of C, decrease in A and B. Expansion is not possible for 
A < 0. For R, S < 0, expansion is possible for large values of other parameters. 
For B < 0, expansion is possible if Vq and A are small, R and F are negative. 
In this case, S can be both negative and positive. Again for expansion to occur, 
the parameter space {R, F, S < 0), {A, B > 0) and {A,B,R >0), (S', F < 0) 
not allowed. But if we look at the expression of Fq, then R, A cannot take 
negative values. This will make the expression imaginary, which is unphysical. 
For F < 0, S' > 0 but small, expansion is possible if p lies between 1 and 3. 

• Fig. 29(b) shows the plot of the behavior of the potential with time for small 
field hilltop potential. This graph has been obtained with the help of Eqn. 
(7.58) with parameter values Vq = 2.7xlO“^Mp, p = 4, Fq = 54Mp, A = 
O.lMp, B = 0AMp,C = 500, R = SOM^, S = 0.3 /9 = 0.05. Correct form of 
potential is possible for both negative and positive values of /3. In this case, 
the condition 

\[r^+2^^q > 2{A + B/2) (7.97) 

should always be satisfied for getting real outputs which we can conclude from 
Eqn. (7.58). Larger values of S makes the potential flatter for longer time. 
Larger values of other parameters make the potential fall more linearly and 
decreases its slope. Both negative and positive values of B are allowed provided 
the above condition is satisfied. S' < 0 not allowed because we do not get the 
correct nature of the potential. 

• In Fig. 29(c), we have plotted Eqn. (7.68), which results in the contraction 
phase of the universe provided we choose the value of the constants accordingly. 
This plot has been obtained for A = lOM^, B = 20Mp, F 3 = 150Mp, F 4 = 1. 
Output is possible provided 

VAR/2>2{A + B/2), (7.98) 

which we can also conclude from Eqn. (7.68). Contraction is possible provided 

B,R>A (7.99) 

and 

A,R,B>0. (7.100) 

Higher values of the parameters increases the amplitude of expansion. If we 

compare Fig. 29(a) with Fig. 29(c), we find that there occurs a net increase in 
the amplitude of the scale factor after one expansion-contraction cycle. 
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• Fig. 29(d) shows the plot of the potential during the contraction phase given 
by Eqn. (7.65). This plot has been obtained for Vq = p = 1, F 2 = 

24M“^, F 3 = IMp, A = 2Mp, B = —2Mp, R = {3 = —0.9. Larger 

values of the integration constants F 2 and F 3 makes the rise more nonlinear. 
For /3 > 0 rise of the potential is possible only for very large values of F 2 . 
Larger values of other parameters increases the amplitude of the potential. 

7.8.2 Case II: Natural potential 

All the graphs in this section and in the following sections have been plotted in units 
of Mp =1, iLo = 1) c = 1, where Mp is the Planck mass, Hq is the present value of 
the Hubble parameter and c is the speed of light. 


a vs t for Einstein Gauss 
Bonnet gravity modei 



V(^(t)) vs t (^«f)for Einstein Gauss 
Bonnet gravity model 





(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the poten- 
factor with time during expansion phase for tial during expansion phase for (j) « f with 
(j) « f with Vo = SxlO-^Mj), Fo = I, A = Vo = SxlO-^M^, Fq = I, A = 

C = 1, R = 40.6M2. 


0.8M2, B = A.2M^ 


= SxlO-^M^ 
0.1M2,C= 10 


i? = 40.6M2, / = 0.1Mp 


O.lMph B = 


Figure 30. Graphical representation of the evolntion of the scale factor and the potential 
dnring the expansion for Einstein Ganss Bonnet gravity model. 


Fig. 30 shows the evolution of scale factor and potential during expansion phase 
for Einstein Gauss Bonnet gravity model. We can draw the following conclusions: 

• Fig. 30(a), shows the plot of the scale factor in the small held limit for natural 
potential given by Eqn. (7.74) with Vq = 8 xlO“^Mp, F 5 = 1, A = O.SM^, B = 
4.2M2, C = 1, R = 40.6M2. 
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• Detail graphical analysis show that the amplitude of expansion increases with 
increase in the value of R. From Eqn. (7.74), we can conclude that i?, A < 0 
not possible. Expansion is possible for B < Q. But large negative values of B 
expansion is not possible. The nature of the graph is almost independent of 
the value of C and Vq. 

• Fig. 30(b) shows the plot of the behavior of the potential with time for small 
field natural potential. This graph has been obtained with the help of Eqn. 
(7.74) with parameter values Vq = 8xlO“^Mp, Fq = 1, A = O.lMp, B = 
0.1Mp,C = 10, R = 40.6Mp, / = O.lMp. The conclusions regarding the 
allowed parameter space for the expansion of the scale factor holds true for 
this case also. Very large values of Vq and A are not allowed since they cause 
uneven oscillations in the potential. 

• For large field case given by Eqns. (7.83) and (7.81), the behavior is similar, 
hence have not been shown here explicitly. 

• If we compare Fig. 30(a) with Fig. 29(c), we find that there occurs a net 
increase in the amplitude of the scale factor after one expansion-contraction 
cycle. 

7.8.3 Case III: Coleman-Weinberg potential 

All the graphs in this section and in the following sections have been plotted in units 
of Mp =1, Hq = 1, c = 1, where Mp is the Planck mass, Hq is the present value of 
the Hubble parameter and c is the speed of light. 

Fig. 31 shows the evolution of scale factor and potential during expansion phase 
for Einstein Gauss Bonnet gravity model. We can draw the following conclusions: 

• Fig. 31(a), shows the plot of the scale factor in the small field limit for super¬ 
gravity potential given by Eqn. (7.94) with Vq = 10“^Mp, Fg = IMp, Fig = 
1, A = O.IM^, B = O.IM^, C = 1, R = AAMl, f3 = 1.5, a = 0.33. 

• Detail graphical analysis show that the amplitude of expansion increases with 
increase in the value of R. But the expansion is almost independent of a, f3, Fg. 
Negative values of B gives large expansion. For larger values of R, Fg cannot 
take large values. Expansion occurs for both negative and positive values of 
(3. Larger values of Vg makes the graph more linear. For A < 0, expansion is 
not possible. However, B can take negative values, but both B and R taking 
negative values not allowed. 

• Fig. 31(b) shows the plot of the behavior of the potential with time for small 
field supergravity potential. This graph has been obtained with the help of 
Eqn. (7.91) with parameter values Vg = 1.2xlO“'^Mp, Fg = O.OMp, A = 
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a vs t for Einstein Gauss 
Bonnet gravity modei 


V(^(t)) vs t (^«Mp)for Einstein Gauss 
Bonnet gravity model 
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100 
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(a) An illustration of the behavior of the scale (b) An illustration of the behavior of the po- 
factor with time during expansion phase for ^ < < tential during expansion phase for (j) « Mp 
Mp with Vo = 10-®M4, Fq = IMp, Fio = with Vq = Fg = 6.6Mp, A = 

1, A = 0.1M2, B = 0.1M2, C = I, R = 0.1M2, B = 0.1M^,C = 10, i? = 20M^, a = 
4.4M2, /3 = 1.5, a = 0.33. 0.38, j3 = 0.85 . 

Figure 31. Graphical representation of the evolntion of the scale factor and the potential 
dnring the expansion for Einstein Ganss Bonnet gravity model. 

O.lMp, B = 0.1Mp,C = 10, R = 20Mp, a = 0.38, (3 = 0.85. The conclusions 
regarding the allowed parameter space for the expansion of the scale factor 
holds true for this case also. Very large negative values of (3 are not allowed. 
Expansion is possible only for small negative and positive values of (3. Other 
parameter values only increases the amplitude, keeping the nature of the graph 
unchanged. The plot is almost independent of the value of C and R. 

• If we compare Fig. 31(a) with Fig. 29(c), we find that there occurs a net 
increase in the amplitude of the scale factor after one expansion-contraction 
cycle. 

8 Conclusion 

In this paper, we have further explored the role of hysteresis in making cyclic models 
of universe as an alternative proposal to inflationary paradigm. The idea, originally 
proposed by the authors in refs. [10, 11], have been studied by us for a wide variety of 
cosmological models, especially, higher dimensional gravity setup. The basic analysis 
for getting the conditions which lead to an increase in expansion maximum, remains 
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the same as discussed in refs. [10, 11]. The most interesting outcome of this analysis 
is the dependence of the bouncing and turnaround conditions on the various model 
parameters. Through this analysis we hnd that the phenomenon of hysteresis is very 
robust and is eternal for all of these models. We study essentially those effective held 
theoretic models which can give rise to both the conditions for bounce and turnaround 
and at the same time also satishes the observed features of the present universe. 
Like in the previous analysis as performed in the ref. [11], here also we hnd that 
inhationary conditions are not necessary for causing cosmological hysteresis. Hence 
the models analyzed in this paper can also be treated as an alternative prescription 
to inhation. Most of the results indicate that the value of the scale factor maximum 
and minimum after each cycle, depends not only on the signature of the hysteresis 
loop integral but also on the relative amplitudes of the model parameters, or in 
other words, we can hne tune the model parameters and get an amplitude increase 
after each cycle even if the signature of the hysteresis loop integral is positive. While 
doing the analysis we have not constrained the potential by any particular form. The 
potential can have any general form like a power series, oscillatory, etc, but with well 
dehned minimum/minima which is essential for generating the required randomness 
or mixing of the held in the phase space {(j), 0) so that its value during contraction 
and expansion are uncorrelated. 

We have also tried to hnd the nature of dependence of the scale factor on time by 
directly solving the equation of motion for the scalar held under diherent scenarios 
and varying forms of the cosmological potential. Exact solutions of the equation of 
motion for the scalar held being highly complicated, we applied certain valid and 
possible approximations while doing the analysis by redehning the held in terms of 
other dynamical parameter, in order to get an approximate analytical expression for 
the dependence of the scale factor on time. This in turn helps us to directly guess 
the nature of the behavior of the scale factor as well as the scalar held during each 
cycle for diherent models. We have also derived an explicit expression for work done 
per cycle, hence have shown that its value comes out to be non zero for the cases 
which we have studied explicitly, provided we choose the parameters of the models 
accordingly. It shows that the phenomenon of hysteresis and a cyclic universe with 
an ever increasing scale factor, can be produced by a wide variety of models and 
potentials. This will also help us to constrain the parameters of the models and the 
potentials, such that we get the required results compatible with observations. 

The future prospects of our work are mentioned below; 

• Through this analysis we have seen the beautiful correlation between purely 
thermodynamical principle and relativistic models and how the former can be 
used for extracting interesting results from the later. But the models that we 
have considered are the variants of minimally coupled gravity frameworks. It 
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would therefore be interesting to investigate what new features arises once we 
relax this constraint. Also we can check whether modihed gravity models like 
for variants of f{R) gravity, two brane-world model in presence of the Einsitein- 
Hillbert term and the Einstein-Hillbert-Gauss-Bonnet gravity setup succeeds 
in generating a cyclic universe with increasing amplitude of expansion. 

• One would also like to ask what other observational signatures we can get from 
such models which can be tested using CMB. We have also not yet verihed 
whether these models produces a universe with the right amount of anisotropy 
and inhomogeneity present in the universe. In future we plan to connect these 
analyses with CMB observations, by rigorous study of the cosmological pertur¬ 
bation theory in various orders of metric fluctuations and computation of two 
point correlations to get the expressions for scalar and tensor power spectrum 
in this context. Hence we extend the study of this paper to compute the pri¬ 
mordial non-Gaussianity in CMB from three and four point correlations. We 
also plan to derive the explicit expression for various modihed consistency re¬ 
lations between the non-Gaussian as well as other cosmological parameters in 
the present context. 

• Also, our present analysis have been performed for three different potentials. 
It would be interesting to check whether we can formulate any general form of 
a potential, which can be used to repeat the analysis of this paper for all the 
models. 

• We also carry forward our analysis in the development of density inhomo¬ 
geneities, which is the prime component to form large scale structures at late 
times. Also the specihc role of cosmological hysteresis in the study of cosmolog¬ 
ical perturbations i.e. for interacting/decoupled dark matter and dark energy 
have not been explored at all earlier. We have some future plan to do some 
computations from this setup. 

• Further using the reconstruction techniques we want to study the generic fea¬ 
tures of scalar held potentials in the framework of cosmological hysteresis. 

Acknowledgments 

SG would like to thank Department of Theoretical Physics, Tata Institute of Funda¬ 
mental Research, Mumbai for providing me Visiting (Post-Doctoral) Research Fel¬ 
lowship. SG take this opportunity to thank sincerely to Prof. Soumitra SenGupta, 
Prof. Sayan Kar, Prof. Sandip P. Trivedi, Prof. Shiraz Minwalla, Dr. Subhabrata 
Majumdar and Dr. Supratik Pal for their constant support and inspiration. SG 
take this opportunity to thank all the active members and the regular participants 


- 162 - 



of weekly student discussion meet “COSMOMEET” from Department of Theoretical 
Physics and Department of Astronomy and Astrophysics, Tata Institute of Funda¬ 
mental Research for suggesting various crucial issues at the initial stage of the work, 
which hnally helped us to improve the quality of the work. SC also thanks Indian 
Association for the Cultivation of Science (lACS), Kolkata and Physics and Applied 
Mathematics Unit (PAMU), Indian Statistical Institute (ISI), Kolkata for extending 
hospitality during the work. Additionally SC take this opportunity to thank the 
organizers of 28th lAGRG Meeting, 2015, Raman Research Institute (RRI) and In¬ 
dian Institute of Science (IISG) for providing the local hospitality during the work. 
SB sincerely thanks Prof. T. P. Singh for constant support and inspiration. Last 
but not the least, we would all like to acknowledge our debt to the people of India 
for their generous and steady support for research in natural sciences, especially for 
theoretical physics. 

9 Appendix 

9.1 Hysteresis from RSII brane world model 

Below we have quoted the results for RSII brane world cosmology with time like 
extra dimension necessary for causing the bounce, which was studied in detail by the 
authors of [11]. 

The Friedmann equations in this model are given by: 

a 

a j 0 Pc J 

W + W" = - = -^|(/5 + 3p)-^(2p + 3p)| . (9.2) 

a 3 [ Pc J 

9.1.1 Condition for bounce 

At bounce, by setting the Hubble parameter H = 0, we get: 

Pb = Pc (9.3) 

The mass at bounce (neglecting the constant factor) is given by: 

Mb = pbol- (9-4) 

Therefore, the inhnitesimal change in the mass content at bounce is given by: 

5M = 5{pcabY ■ (9.5) 

From energy conservation, we get: 

5M + 5W = 0, (9.6) 


mu 


-p 




(9.1) 
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where 6 W is the work done during each expansion-contraction cycle which is given 
by f pdV which includes contribution from the area of the hysteresis loop. 

Further setting 


6 M = -6W = 



(9.7) 


we get the expression for change in amplitude of the scale factor at each successive 
cycle as: 




3 _ 







(9.8) 


Thus we see that the change in minimum of the expansion factor depends upon the 
work done and the density at bounce. 


9.1.2 Condition for acceleration 


From Eq. (9.2), condition for acceleration is given by: 

47rG(pc+ Pc) > 0 (9.9) 

i.e. bounce can be obtained without violating the strong energy condition. 

9.1.3 Condition for turnaround 


Though the authors of [11] have suggested different ways of getting the condition for 
turnaround, but in there analysis they have used the presence of density which take 
negative values at late times to be the cause for turnaround. The authors of [11] 
have assumed that the late time behavior of the universe is governed by: 

a"' 

where k = and A > 0, n < 2. Here H = A<0, n = 0 corresponds to a negative 
cosmological constant, whereas A = 1, n = 2 describes a universe which is spatially 
closed. Setting the Hubble parameter H = 0, we get the condition for turnaround 
as: 

H 

= Ti 
at 

where pt is the density and at the expansion factor at turnaround. They are connected 
via the following equation: 


(9.11) 


pt = 


3Mi 

2 



(9.12) 


where pt and at are the density and scale factor at turnaround respectively. 

The change in amplitude of the scale factor after each successive cycle is given 
by: 



which is responsible for turnaround. Here two extreme physical situation correspond 
to: 
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• (i) the negative cosmological constant {n = 0) for which 

= - jf pdV , (9.14) 

• (ii) the spatially closed universe (n = 2) for which 

^Omax = - (^tax = 

9.1.4 Condition for deceleration 

Turnaround occurs only if; 

p + 3p>0 (9.16) 

Therefore, turnaround can be obtained without violating the energy condition in 
RSII setup. 

9.2 Exact expressions for 0 integrals 
DGP model-Hilltop potential 


i) The exact solution to the left hand side of the integral in Eq. (4.70) for early 
times is given by 


Xi = 


(-4(i + 0 (1 + (<i>lMtyi) 

4(-2 + p)/3p 

{(j)/M^f-P Hypergeometric2Fl -{(j)/M4)P/3 


A{-2 + p)j3p 

ii) The exact solution to the left hand side of Eq. (4.81) for late times is given by 


(9.17) 


I, = 


(-m - 


M-2+p) (Dg+ 


1/2 


WM^r-{Wlpp + 


Da 


-Hypergeometric2Fl 


2Vo/3 


Dg 


4(-2 + p) {Dg + 


1/2 


(9.18) 


where 


Dg = ^ 


2 En 


3M| 


1 

H— 


(9.19) 
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DGP model- Natural potential 


i) Exact solution to the integral on left hand side of Eq. (4.110) for early times 
is given by 


X. = 


m 


1/2 r 


1 + COS I j 


1 1/2 r 


-InIcosI jy I I +liilsinl - 


sec 


2 / 

(9.20) 

ii) Exact solution to the integral on left hand side of Eq. (4.115) for late times is 
given by 




where 


^ + Zlpcos(2 


£) 




n - , 1 




Dp + 
(9.21) 


(9.22) 


2Vo 

zmJ 


Cosmological constant model- Hilltop potential 


Exact solution to the integral on left hand side of Eq. (4.110) for expansion is given 
by 


X. = 






4(-2 +p) (Dk + 


1/2 


{(j)/Mi)‘^ I 1 + -Hypergeometric2Fl 






Vo/3 

3M? 


2 ’ p ’ p 


Da 


4(-2 +p) (^Da + 


1/2 


(9.23) 


where 


Do = ^ + - 
^ 3M| 3 ■ 


(9.24) 


Cosmological constant model- Natural potential 
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Exact solution to the integral on left hand side of Eq. (5.96) for expansion is given 
by 


Xfi — \ Di 


where 




3M| 


0,-1 
tan 


^+ Df cos (j 


Df -hh 

/ zmI 




^ + D^cos(f 


3M| 


Df + -bt, 

/ ^ 3M? 


(9.25) 


Df = 


Vc 


3M2 


A 

4^3- 


(9.26) 


LQG model- Hilltop potential 


Exact solution to the integral on left hand side of Eq. (6.44) for early time expansion 
is given by: 


l7 = - 


Vo 

(0/Mp)2-yi - ^(1 + (0/Mp)P/3)AppellFl 


(9.27) 

-1 + -14’-( 




1 I Vq I Vo(<i>/Mp)Pl3 \ 1/2 

-2 + P){ 1 Pp 


-l + ('^) 

^ Pc ' 


LQC model- Natural potential 


Exact solution to the integral on left hand side of Eq. (6.68) for early time expansion 
is given by 


Xg 





2En 


Pc 


cos(0/2/)+ 



^l-2|icos(^/2/)' 

/l _Yo_Vo 
V Pc pc W ^ J 


(9.28) 
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9.3 Expressions of the constants in work done analysis 


9.3.1 Dvali-Gabadadze-Porrati (DGP) brane world model 


Hilltop potential 


/ 


0-6 — 




\3Mi) 

2(1 + /3)V2 


/ 


07 — 


04 = 


/3p 


Vo 

3M5 


1/2 \ 


2(1 + /3)V2 


Vo 

3M? 


2 I ^0 

^ ' 3M| 


1/2 


(1+/?)^/" , J_ 

/3p 2r, 


I# .3 /^(2+ (!-*» 


o 


05 

Cli 


2 V 07 


2y/aj 


1/2 

02 = O;^ 

3/2 / 

03 = 0/ jr ^ 




Vo 

_3M|_ 

Vo A ^ 

3MJ y /3p 2rc 


\ 

/ 


(9.29) 


Natural potential 
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9.3.2 Cosmological constant dominated Einstein gravity model 


Hilltop potential 


h 

h 

h 

bb 

be 

h 

bi 

b2 

bs 




V3AP 


_Vo_o2 2 


2 (^(1 + /^) + t ) 


Bn 


\3M2j 


/3Y 


4 (^(1 +/^) + t ) 

2// 




Cn 


1/2 


/3l/2^(^) 

^ 2 ^3M2^ 

f,3(h+l2)^/{2h) 

Tu 

^i\/ 3 /(2/3) 

VW2 


(9.31) 


Natural Potential 


a) 0 << / 


Ys 

b9 



Vn 

3M2 


Vn 

3M2 


+ 



1/2 




(9.32) 
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h) (f)» f 


bio 

bi3 

h 

h 

bi2 

bi5 


A 

\3M? 3 


P 

Vn 

3M2 


1/2 


2( -X0_ I A'] 

^V3M2 ^ 3^ 


3M2 


Vn 

3M2 


\3M2 ^ 3) 

1 , hii = h/bio, &14 = Bgi^- 

hh ~ hhbio + ^ 10^6 


-3 + A), 


bioh 


- 1 


lo^e 


(9.33) 


9.3.3 Loop Quantum Gravity model 


Hilltop potential 


Qi = 


q2 = 


Vo 

3M2 


Vb 

3M2 



<^2 — 
(is = 

do 

(iy = 
dg = 


\IMI 

(1 + M/)/V6 

1 + M2 


^3g3pc/2 

Vr+^' 


(9.34) 
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dll — 

di2 = 

dig = 

di8 = 

d20 = 

di3 = 
di4 = 
di5 = 
die = 

di7 = 


SEjeP^il + + pj 

Ak{l + Mp)V^ 

1 

2\/E2 

E'2 

1^ Efpc 

V 4 

M 

4g2 

(h 

2 \ q2 

\/^Q2 

2 

12e5f/(2<?2)^ 


d[j = V3El{12M^ - p,) 


Natural potential 


(^22 — 


K) 

3M2 


Pc 


1/2^ 


1/2 


C?21 — 


6Eqk 


d22 
E^ 

d23 = ^ Pc(12Mp — Pc) 


2d. 


22 


P _ ^gPc 
ao4 — 


2d 


22 


9.3.4 Einstein Gauss Bonnet Gravity brane world model 
Hilltop potential 


1/2 


9i = 


92 = 


' E + 4:JjVo-{A + B/2y/^ 


S 


p{S-l)Eo 


-9i 


93 = 


\^lR + AiJ^Vg-{A + B/2))S 


-9i 


\i^lR + A^I^Vg-{A + B/2))S 


(9.35) 


(9.36) 


(9.37) 


171 



94: = 



g^ = -2 + D' + D'F^D" 

9, = D'D" 

97 = D‘^D" 

/i = 210^1(28^4^^^^ + 15^4^5(3^4^5 - ^ 90)9197 + 9^6(3^4^5 - ^ 949^96 + 91)97 + 9^4^? 
/2 = hQAglgl{lAgiglgl + bgli^g^g^ - ?> 9 &) 9&97 + 8(85(45(5 - 255)5'? 

/s = 4205?5e(14545556 + 854(65455 - 255)5657 + 8545?) 

fi = -1205?5^ (28545^56 + I8545557 - 85457) 
h = 8155|56(45?56 + 57) 

fe = 28 O 5 I 5556 
/r = 2854^56 


/8 

/g 


252 O 5557 

3 ( 91 + 32 )^/( 253 ) 

773 _ 

a/^ 


v/l + A(l + B)Fl 1^- 


■1 + 


aO 


v/rTA(l + S) 


/lO — S'! ^ 8 /( 253 ) 

/ii = VW2 (9.88) 


Natural potential 


„ 3 (a/1 + ^4 + a/ 1 + Ai? — ctC'^) 

/12 = -^5-- 


(9.89) 
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